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Molecular  hydrogen  is  interesting  to  investigate  in  a  porous  material.  It  is 
strongly  quantum-mechanical  with  respect  to  its  translational  and  rotational 
degrees  of  freedom  as  well  as  its  ortho-para  conversion.  Zeolite  (13X)  was 
selected  as  porous  material  for  its  almost  monodisperse  pore  diameter  of  13  A. 

The  studies  mainly  focussed  on  Nuclear  Magnetic  Resonance  data.  The 
hydrogen-zeolite  system  was  probed  in  coherent  pulse  and  continuous  wave 
(cw)  mode  at  268  MHz. 

The  pulse  data  revealed  a  strong  peak  in  the  transverse  relaxation  time  T2 

at  about  10  K.  A  similar  peaking  was  observed  for  the  nuclear  spin  echoes. 
This  behavior  was  interpreted  in  terms  of  supercooling  and  implies  that  hydrogen 
in  zeolite  solidifies  about  5  degrees  below  its  bulk  melting  point.  Two  features 
are  especially  remarkable.  Firstly,  the  T2  was  found  to  be  nonhysteretic  and 

emphasizes  a  purely  energetically  driven  process.  Secondly,  the  supercooling 
transition  temperatures  exhibit  a  linear  parahydrogen  concentration  dependence 
of  -0.10±0.01  K/%.  This  implies  that  parahydrogen  can  be  supercooled  more 
easily  than  orthohydrogen. 
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The  continuous  wave  data  show  a  glass-type  lineshape  of  130  kHz  width 
at  4.2  K  with  a  =  0.45  ±0.3.  It  is  a  combination  of  broad  wings  with  a  narrow 
component  in  the  center.  This  behavior  is  displayed  upon  warming  from  4.2  K 
when  the  narrow  line  disappears  and  the  broad  contribution  remains.  No  sharp 
transition  at  10  K  is  observed.  The  line  broadens  by  a  factor  of  1.5  due  to 
orientational  ordering  when  the  sample  is  cooled  to  0.8  K.  This  latter 
observation  indicates  that  the  pulse  data  peaks  at  10  K  do  not  originate  from 
orientational  ordering. 

The  ortho-para  conversion  was  extracted  from  the  time  dependence  of 
the  lineshape.  The  orthohydrogen  contents  was  followed  over  1600  hours  which 
yielded  two  different  conversion  constants.  For  the  first  500  hours  the 
conversion  is  characterized  by  /c,  =0.425 ±0.006  %/h  which  is  only  a  fourth  of 

the  bulk  value  and  is  explained  by  the  reduced  number  of  nearest  neighbors  in 
the  constrained  geometry.  For  later  times  the  conversion  accelerates  to 
k2  =2.21  it 0.075  %/h  and  is  explained  by  clustering  effects.  The  relative  shape 

of  the  line  changes  over  time  as  the  narrow  component  of  the  line  decays  faster. 

Isochoric  pressure  versus  temperature  data  were  recorded  while  the 
hydrogen  sample  was  taken  out.  The  nature  of  a  two-component  system  was 
confirmed.  The  adsorption  energy  of  hydrogen  on  the  zeolite  surface  was 
determined  to  be  270  ±50  K. 

Hydrogen-deuteride  was  also  adsorbed  into  zeolite.  The  interpretation  of 
the  pulse  data  peaks  in  terms  of  supercooling  was  supported  by  the  fact  that  HD 
showed  a  similar  behavior  without  orientational  degrees  of  freedom. 

Similar  but  more  extreme  results  were  obtained  from  a  second  zeolite 
(5A)  with  even  smaller  pore  size  (9  A)  in  agreement  with  the  previous  data. 
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CHAPTER  1 
INTRODUCTION 


The  system  to  be  discussed  in  this  dissertation  is  hydrogen  confined  to 
zeolite,  a  porous  material.  Experimental  investigations  were  carried  out  by 
means  of  Nuclear  Magnetic  Resonance  at  low  temperatures.  These  four 
characteristics  -  molecular  hydrogen,  zeolite,  low  temperatures  and  NMR  - 
determine  the  scope  of  this  thesis.  Every  one  of  these  will  receive  special 
attention  in  the  following  chapters.  In  addition,  the  effects  that  are  expected  and 
observed  will  be  discussed.  The  introduction  will  lead  into  the  subject  by 
explaining  some  general  principles  that  are  important  for  the  understanding  of 
the  results  and  provide  the  background  for  the  remainder  of  the  thesis.  The 
overview  is  also  meant  to  motivate  the  performed  experiment. 

For  the  understanding  of  many  effects  it  is  necessary  to  be  familiar  with 
the  special  properties  of  hydrogen.  Hydrogen  is  the  simplest  existing  diatomic 
molecule  and  can  therefore  be  used  as  a  model  molecule  to  understand 
fundamental  interactions  and  phenomena.  When  it  is  cooled  below  13.8  K  it  is 
also  the  simplest  molecular  solid.  The  beauty  of  hydrogen  is  that  the  involved 
interactions  are  understood  in  terms  of  first  principles. 

1.1.    Quantum  Nature  of  Hydrogen 

Hydrogen  is  special  in  its  properties.  It  is  a  very  interesting  substance  as 
it  is  explicitly  quantum-mechanical  in  its  nature.  There  are  three  realizations  of 
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this  fact:  the  existence  of  two  species,  the  large  zero-point  motion  (ZPM)  and 
the  quantum  rotor  property. 
1.1.1.  Ortho-  and  Parahvdroaen 

One  of  the  most  fascinating  properties  of  molecular  hydrogen  is  the 
explicitly  quantum  mechanical  classification  into  two  species.  Hydrogen  occurs 
naturally  as  a  diatomic  molecule.  Its  nuclear  constituents  are  spin-1/2  protons 
with  the  electronic  spin  contributions  paired  in  a  symmetric  1X*  molecular  ground 
state.  The  wave  function  can  be  written  as  a  product  of  nuclear,  rotational  and 
vibrational  wave  functions.  The  vibrational  groundstate  is  symmetric  and  does 
not  influence  the  considerations.  The  quantum-mechanical  requirement  of  a 
totally  antisymmetric  wavefunction  for  two  indistinguishable  fermions  can 
therefore  be  realized  in  two  distinct  ways,  leading  to  two  distinguishable 
molecular  species:  ortho-  and  parahydrogen.  Orthohydrogen  has  a  symmetric 
nuclear  spin  wavefunction  (/  =  1)  and  an  antisymmetric  orbital  wavefunction 
(J  odd).  The  nomenclature  is  such  that  the  species  with  the  highest  spin 
quantum  number  is  called  "ortho."  Parahydrogen  has  an  antisymmetric  nuclear 
wavefunction  (/  =  0)  and  a  symmetric  orbital  wavefunction  (J  even).  It  is 
important  to  notice  that  both  species  of  hydrogen  are  bosons. 

Heteronuclear  molecules  such  as  hydrogen  deuteride  (HD)  do  not  display 
distinct  species  because  of  the  distinguishable  nature  of  their  constituents. 

Because  of  the  small  moment-of-inertia  /  for  hydrogen,  the  separation  of 

*2 

the  rotational  states  is  very  large:  Ej  =  Bj  J(J  +1)  with  Bj  =  —  =  85.4  K.  The 

quantum  number  J  is  therefore  a  meaningful  entity  at  low  densities.  The 
intermolecular  interactions  which  tend  to  mix  different  rotational  states  are  very 
weak,  and  the  ground-state  is  the  J  =  0  parahydrogen  state.  This  still  holds 
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when  the  surface-molecule  interaction  of  the  restricted  geometry  is  included,  as 

discussed  in  Chapter  2. 

1.1.1.1.        Quadrupolar  interaction 

Parahydrogen  is  spherically  symmetric  with  no  electric  or  magnetic 
moments.  The  lack  of  nuclear  spin  degrees  of  freedom  for  parahydrogen  makes 
it  also  undetectable  by  Nuclear  Magnetic  Resonance.  This  is  in  contrast  to 
orthohydrogen  with  spin  1  which  also  has  an  electric  quadrupole  moment.  The 
molecules  tend  to  orient  at  low  temperatures  in  order  to  minimize  the  anisotropic 
interaction 


which  occurs  when  the  molecular  axes  are  aligned  at  a  90  degree  angle,  a  T" 
configuration  between  quadrupoles.  Here  A  is  the  lattice  parameter,  R  the 
intermolecular  distance,  Fan  angular  function  of  the  angles  measured  relative  to 
the  connecting  line  between  molecules,  and  r  =  0.8  K. 

The  order  parameters  specifying  the  degrees  of  freedom  for  this 
interaction  are  components  of  the  second  rank  quadrupolar  tensor  with  respect 
to  the  local  reference  axes  (x,  y,  z) 


which  leads  to  two  additional  parameters  besides  the  three  local  reference  axes. 
The  reason  for  only  two  additional  parameters  is  the  vanishing  orbital  angular 
momentum,  which  means  that  (J,)  =  0  is  "quenched."1  This  is  the  case  when 
time-reversal  symmetry  is  not  broken  which  can  be  assumed  here,  even  though 
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a  magnetic  field  is  applied.  The  ratio  of  the  magnetic  interaction  energy  (mK)  to 

h2 

the  involved  rotational  energies  (Sj=  — =  85.4K)  is  the  reason  for  this 

approximation:  the  magnetic  energy  is  only  of  negligible  influence.  The 
argument  for  "quenching"  is  related  to  the  non-degenerate  nature  of  the 
rotational  groundstate.  A  non-degenerate  groundstate  must  be  real  otherwise 
real  and  imaginary  parts  of  the  wavefunction  would  separately  be  solutions, 
which  contradicts  non-degeneracy.  The  angular  momentum  operators  are 
imaginary  and  the  expectation  value  (J,)  with  /'  =  x,y,z  is  therefore  purely 
imaginary.  On  the  other  hand  it  must  also  be  real  due  to  hermiticity  and  the  only 
expectation  value  satisfying  both  conditions  is  zero.  At  higher  temperatures  this 
requires  more  careful  consideration  because  the  Jz  =  ±1  are  also  occupied  and 
could  contribute  to  the  rotational  Zeeman  energy.  This  is  not  the  case,  however, 
because  the  transition  frequency  amongst  these  states  is  much  higher  than  the 
dipolar  frequency  and  is  therefore  averaged  out.  All  off-diagonal  matrix 
elements  are  also  zero  by  the  properties  of  angular  momentum. 

One  additional  order  parameter  for  orthohydrogen  is  the  alignment 

(7=0^  =  ^(3^-2) 

For  Jz  =  0  the  alignment  is  <r  =  1,  which  depicts  a  prolate  ellipsoid  probability  for 
the  wavefunction.  For  the  case  of  an  oblate  ellipsoid  where  Jz  =  1  the  alignment 
is  cr  =  -1/2. 

The  other  additional  order  parameter  is  the  eccentricity  described  by 


1  =  (j2*-4) 
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and  is  zero  for  any  rotationally  symmetric  body. 
1.1.2.  Zero-Point  Motion 

Hydrogen  is  the  molecule  with  the  lowest  mass.  It  possesses  only  a 
weakly  polarizable  electron  cloud  leading  to  weak  intermolecular  forces.  These 
two  facts  combined  with  the  groundstate  energy  EQ=jho)  of  a  three- 
dimensional  harmonic  oscillator  lead  to  a  large  zero-point  motion  (ZPM).  This 
implies  that  the  molecule's  oscillation  amplitude  is  a  sizeable  fraction,  15%  for 
hydrogen,  of  the  intermolecular  lattice  spacing  even  at  zero  temperature.  This  is 
a  fascinating  feature  and  is  the  cause  for  interesting  properties.  It  also  renders 
calculations  more  difficult  because  a  harmonic  approximation  about  the 
intermolecular  potential  minimum  is  no  more  justified.  Also  correlations  between 
molecules  are  important  to  consider  when  an  understanding  for  the  dynamics  is 
developed. 

In  order  to  gain  insight  into  the  ZPM,  a  harmonic  approximation  is 
presented,  with  all  the  caution  in  place.  The  standard  interaction  potential  used 
is  the  Lennard-Jones  potential 


where  o  is  the  hard  core  radius  and  e  the  potential  minimum.  The  first  part  is  a 
mathematically  convenient  form  of  a  repulsive  interaction  and  is  a  consequence 
of  the  Pauli  exclusion  principle  when  electronic  clouds  start  to  overlap.  The 
second  term  is  the  attractive  Van  der  Waals  potential  due  to  induced  dipole- 
dipole  interaction.  No  permanent  dipoles  are  involved  but  one  dipole  is  created 
by  thermal  fluctuation  fields  and  another  dipole  is  induced  in  a  neighboring 
molecule. 
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A  more  sophisticated  potential  was  empirically  found  by  Silvera  and 
Goldman2 


with 


V(r)  =  exp(a-fir-yr2)  +  f(r)\  £ 

[/=6.8,10  '  'J 


2  ' 

=  exp- 

1.28 Rm  1 

r 

=  1 


for  r  <  1.28/=?„ 


for  r>1.28R. 


where  the  constants  are  a  =  1.71 3,  /?  =  1.5671, /  =  0.00993  and  C6  =-12.14, 
C8  =-215.2,  C10  =-4813.9  and  C9  =143.1  all  in  atomic  units  and  Rm  =3.41  A  at 
the  well  minimum  for  which  the  ninth  order  term  is  excluded. 

The  simple  Lennard-Jones  potential  is  expanded  up  to  second  order  in 
the  displacement  u  relative  to  the  radius  of  minimum  potential  energy  which 
leads  to  a  harmonic  oscillator  with  Hamiltonian 


H  = 


■+—ku2 


2m  2 


with  k  =  ZV" 


where  Z  is  the  number  of  nearest  neighbors.  The  groundstate  energy  of  a 
harmonic  oscillator  is  well-known  to  be 


5=%  =  ^ 
2       2  V  m 


The  corresponding  wavefunction  is  of  Gaussian  form 


with  width  W.  The  width  describes  the  spread  of  the  wavefunction  in  space  and 
the  localization  of  particles  on  the  lattice.  The  expectation  values  for  the  kinetic 
and  potential  energies  over  this  groundstate  wavefunction  are 


The  results  imply  two  things:  firstly,  the  shallow  potential  with  small  curvature 
reduces  the  potential  energy  and  secondly,  the  small  mass  increases  the  kinetic 
energy  raising  its  relative  importance  even  more. 

A  calculation  of  the  minimum  energy  with  respect  to  the  width  leads  to 


in  terms  of  the  relevant  parameters.  It  is  especially  large  for  small  m  and  small 
V",  which  is  the  case  for  quantum  crystals,  such  as  hydrogen. 

Hydrogen  has  a  large  kinetic  energy  relative  to  its  small  potential  energy 
and  would  be  too  localized  for  a  bound  state  so  that  it  reduces  its  kinetic  energy 
by  spreading  out  to  a  higher  potential  energy.  This  process  of  spreading  the 
wave  function  leads  to  a  new  arrangement  of  the  lattice  parameter.  A  guideline 

is  r0-2(u2)V2  where  r0  is  the  distance  of  minimal  potential  and  u  is  the 
displacement  about  a  lattice  site.  If  this  quantity  is  smaller  than  the  hard-core 


3  n2 

4mW2 
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radius  the  particles  spend  appreciable  time  in  each  others'  hard  cores  which  is 
an  unphysical  situation.  The  strong  repulsive  interactions  cost  energy  and  a  new 
minimum  energy  must  be  found.  The  solution  to  the  dilemma  is  an  expansion  of 
the  lattice  from  r0  to  A.  This  procedure  reduces  <T>  and  increases  <PE>,  which 
reestablishes  the  situation  of  classical  solids  and  is  a  compromise  for  a  certain 

A-2(u2V/2 

width  W.  Empirically  it  is  found  that  — - —  arranges  itself  around  unity. 

(7 

1.1.2.1.       De  Boer  parameter 

It  is  interesting  to  define  a  measure  for  the  quantum  nature  of  a  molecule. 
This  can  be  approached  by  rewriting  the  Schroedinger  equation  in  terms  of 
reduced  quantities  by  scaling  the  distance  and  the  potential  with  the  Lennard- 
Jones  parameters 

r  V 
r'=-      and      V  =- 
a  £ 


which  results  in 


2nri£oz 


(V*)2+V 


The  prefactor  of  the  reduced  kinetic  energy  term  expresses  the  relative 
dominance  of  the  kinetic  energy  compared  to  the  potential  energy.  This  is  the 
relevant  quantity  for  quantum  behavior  and  is  defined  as  the  de  Boer  parameter 


v/2m£cr2 
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which  is  large  for  small  energy  and  mass.  This  is  consistent  with  the  earlier 
statements  and  calculations  about  the  nature  of  quantum  crystals.  The 
parameter  is  0.196  for  hydrogen3  which  is  much  larger  than  for  classical  crystals 
with  values  below  0.05. 

1.1.3.  Crystalline  Potential  and  Quantum  Rotor 

At  temperatures  above  the  melting  point  (13.8  K)  orthohydrogen  oscillates 
rapidly  between  its  orbital  magnetic  states  M.  This  generates  an  effective 
spherical  symmetry  as  it  is  perceived  by  its  neighbors.  In  this  respect  it  is  similar 
to  parahydrogen  at  high  temperatures.  At  lower  temperatures  where  hydrogen 
solidifies  this  is  no  longer  true.  The  transition  rate  between  M  states  decreases 
and  the  non-spherical  character  of  orthohydrogen  manifests  itself.  The  average 
potential,  the  crystal  field,  at  the  site  of  an  orthohydrogen  molecule  is  of  lower 
symmetry.  This  is  a  self-enforcing,  cooperative  process  as  more  orthohydrogen 
molecules  with  lower  symmetry  reduce  the  free  rotation  of  others. 

The  crystalline  potential  V(0,<D)  is  a  function  of  the  molecular  orientation. 
As  already  emphasized,  J  is  a  good  quantum  number  even  for  the  case  of 
substrate  interaction.  This  renders  it  only  necessary  to  consider  the  three  pure 
states  \JM)  which  are  spherical  harmonics  of  second  order.  If  the  crystalline 
potential  is  expanded  in  terms  of  spherical  harmonics  it  is  clear  by  orthogonality 
that  for  the  expectation  value  only  second  order  harmonics  of  the  potential 
contribute.  The  potential 

V(S,  co)  =  AZX\  +  AXX\  +  AyX\ 
is  therefore  a  quadratic  function  of  the  direction  cosines 


Xx  =sin<5cos<u 
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Ay  =sin<5sintu 
Xz  =  cos  8 


of  the  molecular  axis.  The  eigenfunctions  to  this  problem  are  for  /'=  x,y,z 


At  this  point  the  quantum  rotor  property  of  hydrogen  becomes  imminent. 
By  this  statement  is  meant  that  orthohydrogen  cannot  be  perceived  as  a  static 
entity.  Due  to  its  large  ZPM  it  is  in  rotation  even  at  low  temperatures.  The  J  =  1 
state  is  therefore  not  a  rigid  ellipsoid  but  its  motion  averages  the  shape  into  a 
sphere.  As  will  be  seen,  the  angle  dependence  of  the  dipolar  interaction  is 
governed  by  (3cos20-1)  which  is  proportional  to  the  spherical  harmonic  of 
second  order.  It  is,  however,  of  little  significance  to  define  an  angle  0  between 
the  connection  line  of  the  protons  with  the  magnetic  field.  Instead  it  is  necessary 
to  define  the  angles  /?,  <p  between  one  of  the  crystal  field  axes,  e.g.  z,  and  the 
magnetic  field  and  average  the  angular  dependence  over  all  space  isotropically. 

The  angle  0  between  the  magnetic  field  and  the  proton  axis  can  be 
rewritten  in  terms  of  the  angles  /?,  q> 


cos 0  =  Xx  sin/?cos  <p  +  Xy  sin /Jsin  q>  +  Xz  cos 0 
=  cos  5  cos  p + sin  <5sin/?cos(  <p  -  co) 


and  the  average  of  the  dipolar  interaction  along  one  crystal  axis  Xz  is  taken  over 
all  space 
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(3  cos2  0  -  i)2  =  j-  J  (3  cos2  0-1)  cos2  <5sin  SdadS 
=  |(3cos2/?-1) 

The  result  is  the  famous  reduction  factor  of  2/5.  It  expresses  that  in  the 
case  of  a  quantum  rotor  the  dipolar  interaction  is  rescaled  to  40  %,  leading  to 
narrower  lines  in  the  quantum  case  than  expected  classically. 

1.2.    Open  Questions  for  Confinement 

This  short  overview  about  interesting  properties  of  hydrogen  gives  a  flavor 
for  the  type  of  experiment  that  can  be  performed  with  hydrogen  in  a  porous 
material.  The  substrate  represents  a  geometric  confinement  and  influences  the 
three-dimensional  properties  of  hydrogen.  The  questions  concerning  the 
changes  that  arise  from  such  confinement  are  various.  Some  of  them  are  as 
follows: 

1)  What  impact  does  the  surface  potential  have  on  hydrogen  molecules  being 
closer  to  the  wall  than  others?  Will  different  environments  result  in  different 
subsystems? 

2)  Is  the  ZPM  in  a  more  localized  environment  dominant  enough  to  introduce 
fluctuations  that  wash  out  some  of  the  expected  effects? 

3)  Does  orientational  ordering  of  quadrupoles  occur  and  if,  at  what 
temperature?  This  may  lead  to  new  results  when  the  surface  potential  is 
assisting  in  the  dynamical  slowing  down  of  rotational  motion. 
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4)  Is  the  melting  transition  influenced  by  increased  fluctuations  due  to  lower 
dimensionality  and  a  pore  size  that  reduces  the  size  of  the  critical  solid  nucleus 
in  the  formation  of  a  solid  as  described  by  nucleation  theory?  What  is  the  impact 
of  grain  boundaries  forming  along  rough  substrate  walls? 

5)  Is  a  glass  phase  possible  as  a  result  of  inhibited  and  disturbed  solidification? 

6)  How  is  the  ortho-para  conversion  influenced  by  the  surface  potential,  the 
dynamical  restriction  due  to  the  geometry  and  the  reduced  number  of  nearest 
neighbors? 

7)  At  last  a  very  challenging  idea  is  brought  forward:  would  it  be  possible  to 
reach  a  superfluid  state  in  hydrogen  if  the  solidification  was  suppressed  to 
temperatures  where  superfluidity  is  conceivable  and  theoretically  predicted? 

This  catalogue  of  questions  is  well  worth  a  thorough  experimental 
investigation  that  might  give  answers  to  these  open  problems. 

Chapter  2  will  expand  on  the  concepts  and  ideas  just  mentioned.  This  will 
allow  for  an  understanding  and  appreciation  of  the  importance  of  the  previous 
questions.  The  conclusion  chapter  7  will  try  to  answer  them  based  on  the 
obtained  data.  Chapter  3  focuses  on  some  special  topics  in  Nuclear  Magnetic 
Resonance  as  they  pertain  to  the  experiment.  Chapter  4  describes  the 
experimental  "hardware"  and  gives  an  introduction  to  the  technical  specialties  of 
the  experiment.  The  experimental  data  are  presented  in  chapter  5  while  chapter 
6  contains  the  computer  simulation  of  the  NMR  lineshape.  Chapter  7  closes  the 
dissertation  with  a  conclusion  that  restructures  the  experimental  data  as  they 
pertain  to  physical  effects  rather  than  a  mere  listing  of  facts  and  explanations. 


CHAPTER  2 
THEORY 


2.1.    Surface  Interactions 


The  properties  of  a  molecular  liquid  and  solid  are  expected  to  be  modified 
in  contact  with  a  surface.  The  main  reason  is  the  interaction  potential  between 
the  molecules  and  the  surface  which  has  to  be  incorporated  in  a  theoretical 
treatment.  The  impacts  of  the  surface  are  various.  The  potential  energy  arising 
from  an  attractive  potential  such  as  the  Van  der  Waals  potential  restricts  the 
molecular  motion,  tending  to  localize  the  particles  in  potential  wells.  For  particles 
adsorbed  on  a  surface  the  number  of  nearest  neighbors  is  reduced,  which  is  a 
direct  consequence  of  the  geometrical  constraints.  This  in  turn  reduces  the 
dominance  of  the  intermolecular  interaction  therefore  increasing  the  relative 
importance  of  the  external  surface  potential.  The  two-dimensional  character  on 
the  surface  implies  a  reduced  number  of  accessible  modes  compared  to  three 
dimensions. 

These  facts  should  lead  to  a  modified  phonon  spectrum,  reduced 
diffusivity,  a  distinct  first  monolayer  behavior,  induced  dipole  moments,  reduced 
entropy4  and  less  intermolecular  interaction. 

As  a  consequence,  the  system  could  display  altered  macroscopic  effects. 
These  macroscopic  manifestations  of  the  surface-molecule  interaction  could  be 
modifications  in  viscosity,  orientational  ordering  and  solidification  behavior,  just 
to  mention  a  few. 
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Being  familiar  with  these  facts,  one  must  keep  in  mind  the  dominant  zero- 
point  motion  of  hydrogen  which  is  counteracting  the  localizing  surface  potential 
and  might  reduce  the  effective  surface  potential. 
2.1.1.  Theory  of  Surface-Molecule  Interaction 

Many  experiments  and  theoretical  calculations  have  been  performed  for 
hydrogen,  methane  and  helium  on  vycor,  zeolite  and  two-dimensional 
substrates,  favorably  grafoil. 

A  comparison  between  the  properties  of  zeolite  and  grafoil  is  in  place  as 
the  calculation  to  be  discussed  was  performed  for  grafoil,  but  the  experiment  in 
this  dissertation  was  carried  out  using  zeolite.  The  attraction  potential  of  grafoil 
is  about  a  factor  of  two  larger5  than  the  interaction  strength  of  most  other  solids, 
including  zeolite.  In  addition,  the  surface  is  two-dimensional  in  comparison  with 
a  curved  geometry  in  zeolite.  An  additional  complication  for  zeolite  is  the 
tortuosity  of  the  restricted  topology.  The  question  of  dimensionality  arises  in  this 
context.  The  main  idea  is  nevertheless  similar  and  reduces  to  investigating  a 
system  that  interacts  and  is  partially  bound  to  the  surface  of  a  substrate. 

The  following  calculations  will  lead  to  the  ground  state  energy,  the  two- 
dimensional  phonon  spectrum6  and  the  purity  of  the  zero  rotational  state7  of 
hydrogen  on  a  substrate. 

The  starting  point  for  any  calculation  of  this  type  must  be  the  Hamiltonian 
for  such  a  system.  It  consists  of  the  translational  and  rotational  kinetic  energies, 
the  molecule-substrate  interaction  energy  and  the  intermolecular  interaction 
energy,  neglecting  weak  anisotropic  intermolecular  interactions: 


H  =  HT+HR 


15 


where  M  is  the  mass  of  the  molecule,  L  the  angular  momentum  operator  and  v 
the  isotropic  intermolecular  interaction  which  depends  only  on  the  relative 
distance  between  molecular  centers  of  mass.  The  anisotropic  substrate- 
molecule  interaction  VA  depends  on  the  distance  itself  and  coupies  the  rotational 
and  translational  parts  of  the  Hamiltonian.  For  most  calculations,  however,  the 
wavefunction  is  approximated  by  an  independent  product  of  spatial  and 
rotational  contributions.  This  allows  for  writing  two  independent  eigenequations: 

HT<D  =  ET0 

for  the  translational  part  and 

=  f„¥ 

for  the  effective  rotational  part,  where  the  argument  of  the  anisotropic  substrate 
potential  is  now  the  position  of  the  center  of  mass  averaged  over  the 
translational  motion. 

The  translational  part  is  treated  to  calculate  the  phonon  spectrum.  It  was 
emphasized6  that  the  calculation  must  incorporate  the  third  dimension  of  the 
adsorbed  layer.  The  translational  wave  function  is  therefore  separated  into  an 
in-plane  and  a  vertically  directed  Gaussian  part.  The  ground  state  energy  is  a 
superposition  of  three  energy  contributions 
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E  E2     +  +  E^ 

The  energy  Ez  is  the  kinetic  energy  associated  with  the  motion  in  z- 
direction  plus  the  potential  energy  of  the  laterally  averaged  substrate  potential  in 
z-direction  and  is  evaluated  using  the  z-component  of  the  wavefunction. 

The  E^  term  is  correspondingly  related  to  the  kinetic  energy  in  the  x-y 
plane  plus  the  intermolecular  interaction  energy,  averaged  over  the  z-component 
of  the  wavefunction  and  evaluated  with  the  in-plane  wavefunction.  The  form  of 
the  intermolecular  interaction  is  the  earlier  described  Silvera-Goldman  potential 
omitting  the  C9term. 

The  third  contribution,  E^,  also  involves  both  components  of  the 
wavefunction.  It  represents  the  vertical  potential  energy  due  to  the  periodic 
variation  of  the  adsorption  potential  averaged  over  in-plane  displacements. 

The  solution  to  this  problem  is  determined  variational^  by  a  minimization 
of  the  total  energy  as  a  function  of  the  parameters  in  the  two  Gaussian 
wavefunctions.  This  leads  to  a  two-dimensional  self-consistent  phonon  theory 
because  each  wave  function  contribution  appears  in  two  equations.  This  makes 
the  process  iterative  and  several  iterations  are  necessary  to  obtain  a  self- 
consistent  system  of  equations.  The  process  starts  with  the  determination  of  the 
z-contribution  of  the  wave  function  because  it  appears  solely  in  the  Ez  equation. 
It  has  been  shown  that  the  self-consistent  phonon  approximation,  including 
anharmonicities  from  the  zero-point  motion,  leads  to  a  ground  state  energy 
which  is  by  about  15%  higher.8 

The  ground  state  is  almost  entirely  determined  by  the  Ez  contribution  and 
the  total  energy  is  about  500  K  for  hydrogen  on  graphite.  The  phonon  spectrum 
is  characterized  by  a  phonon  gap  of  46.6  K,  a  maximum  in  the  transverse 
phonon  mode  of  64.9  K  and  a  maximum  in  the  longitudinal  mode  of  83.8  K.  This 
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suggests  that  it  is  more  advantageous  for  the  system  to  rotate  away  from  the 
symmetry  direction.  The  phonon  .density  has  a  large  width  of  48  K  because  of 
the  substantial  zero-point  motion. 

These  calculations  are  backed  by  agreement  with  experimental  data 
which  are  within  10%.  The  conclusion  is  that  the  third  dimension  must  be 
incorporated  in  calculations  dealing  with  molecule-surface  interactions.  It  is 
intuitively  also  clear  that  the  restricted  geometry  limits  the  long  wavelength 
phonons  which  couple  to  the  substrate  phonon  modes  and  are  the  predominant 
excitation  at  low  temperatures. 

Novaco6  initially  performed  the  calculation  for  a  commensurate  phase  but 
extended  it  later  to  the  incommensurate  phase.  The  result  is  that  the  energy  for 
the  incommensurate  phase  lies  only  a  few  Kelvin  (3%)  above  the  commensurate 
phase  groundstate.  This  is  of  importance  for  zeolite  where  no  commensurate 
phase  can  be  expected  due  to  the  surface  roughness  and  curvature. 

Returning  to  the  rotational  part  of  the  total  Hamiltonian,  another  valuable 
piece  of  information  can  be  extracted.7  This  is  related  to  the  mixing  of  rotational 
states  under  the  influence  of  an  additional  nonspherical  potential.  The  ordinary 
rotational  ground  state  is  known  to  be  J  =  0  and  could  be  subject  to  mixing  with 
higher  order  states. 

The  sum  in  the  Hamiltonian  extends  over  single  particle  terms  which  can 
be  solved  individually  with  the  rotational  wave  function  expanded  in  a 
superposition  of  free  rotor  states  \JM).  The  single  particle  energy  contributions 
are  then  summed  up  to  yield  the  total  energy. 

The  use  of  spatial  wave  functions  such  as  in  the  ground  state  calculation 
leading  to  the  phonon  spectrum  (splitting  into  in-plane  and  out-of-plane 
contributions)  already  presumes  a  J  =  0  state.  If  the  results  turned  out  to  be 
different  from  this  assumption  a  new  Ansatz  would  become  necessary.  This  is 
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not  the  case,  however,  since  hydrogen  exists  to  99.9%  in  the  J  =  0  state, 
confirming  the  general  assumption  of  pure  states. 

Furthermore,  as  the  calculation  shows,  hydrogen  molecules  adsorbed  on 
a  substrate  behave  like  slightly  hindered  three-dimensional  rotors.  It  was  also 
confirmed  by  analysis  of  neutron  scattering  data  that  hydrogen  on  surfaces,  such 
as  grafoil,  behaves  like  a  three-dimensional  rotor.9  Further  experimental  support 
for  three-dimensional  behavior  is  superfluid  helium  in  vycor.10  The  three- 
dimensional  network  displays  a  bulk-like  transition  independent  of  fractional 
filling. 

For  rotationally  excited  states  an  interesting  picture  was  deduced  from 
perturbation  theory11  and  numerical  studies.12  Rotationally  excited  states  couple 
with  the  vibrational  states  associated  with  motion  perpendicular  to  the  surface. 
This  coupling  results  in  a  motion  where  the  molecules  line  up  with  the  normal  of 
the  substrate  at  the  maximum  amplitude  of  the  zero-point  motion  and  lie  parallel 
at  the  center.  This  implies  a  mixing  of  the  M  =  0  and  the  M  =  ±1  states  inside 
the  J  =  1  manifold. 

Recent  experimental  data13  suggest  strong  mixing  of  rotational  states  for 
hydrogen  deuteride  (HD)  on  vycor.  Rotational  transitions  of  as  high  as 
SJ  =  3  and  SJ  =  5  were  observed  which  violate  the  bulk  selection  rules.  This 
may  be  a  consequence  of  the  mismatch  of  center  of  mass  and  center  of  rotation 
coordinates  in  HD. 

The  important  result  is  that  the  standard  classification  of  hydrogen  in  ortho 
and  para  species  in  their  pure  forms  still  holds.  The  same  should  also  be  true 
for  zeolite  with  its  weaker  interaction  potential  and  more  three-dimensional 
character.  This  is  valuable  information  for  the  ortho-para  conversion. 
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2.1.2.  Effects  of  Adsorption 

Several  interesting  results  from  other  experiments  and  calculations  wiil  be 
quoted  because  of  their  relevance  to  the  experiment  to  be  discussed  in  this 
dissertation. 

Related  to  the  phonon  spectrum  are  heat  capacity  measurements  of 
helium  in  vycor  indicating  a  decreased  Debye  temperature14  in  two  dimensions. 
This  is  of  importance  for  the  creation  of  phonons  carrying  away  the  released 
energy  in  the  ortho-para  conversion  process. 

The  behavior  of  the  first  monolayer  on  the  substrate  has  been  observed 
to  be  different  from  subsequent  layers.  Superfluidity  in  helium15  decreases  and 
disappears  as  the  coverage  approaches  one  layer.16  This  seems  to  suggest  that 
superfluidity  only  occurs  on  top  of  a  more  localized  first  layer  smoothing  out  the 
surface. 

Experiment17-18  and  theory19  seem  to  agree  upon  the  influence  of  the 
zero-point  motion  on  the  separation  between  hydrogen  molecules  and  the 
adsorbant.  The  expectation  is  that  the  large  compressibility  leads  to  a  reduced 
nearest  neighbor  separation.  This  is,  however,  not  the  case  because  the  ZPM, 
enhanced  by  the  restriction  from  the  pore  environment,  counteracts  the  attractive 
influence  of  the  surface,  hardly  allowing  the  molecule  to  get  close  to  the  surface. 
This  is  illustrated  by  3He  being  further  away  from  the  surface  than  4He  reflecting 
their  relative  ZPM.20  Methane,21  however,  is  forced  closer  to  the  wall  at  low 
temperatures  due  to  smaller  ZPM. 

Induction  of  dipoles  was  observed1722  from  strongly  enhanced  infrared- 
absorption  spectra  of  hydrogen  in  vycor.  The  conclusion  of  induced  dipoles 
follows  from  the  independence  of  the  infrared-absorption  signal  over  time.  This 
means  that  ortho-  and  parahydrogen  give  rise  to  the  same  signal  which  renders 
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the  electric  quadrupole  moment  of  the  ortho  species  irrelevant  and  shows  the 
strength  of  the  surface  interaction. 

A  special  effect  on  the  diffusion  rate  is  the  purely  geometric  restriction.23 
The  parallel-pore  model24  was  proposed  to  describe  a  monodisperse  pore-size 
distribution.  In  this  case  the  diffusion  Ds  is  related  to  the  free  bulk  diffusion  DQ  in 
an  empirical  relationship  with  the  porosity  O  and  the  tortuosity  factor  52S 


The  tortuosity  is  a  measure  of  the  geometric  complexity  of  the  substrate.  A 
numerical  value  was  found25  to  be  5=  5.6  for  hydrogen  in  vycor  of  about  10  A 
radius.  This  is  close  to  the  dimensions  of  zeolite  and  can  give  an  idea  of  the 
purely  geometric  influence  on  the  diffusion  rate.  Using  $  =  50  %  for  the  porosity 
of  zeolite  the  tortuosity  effect  reduces  the  diffusion  by  roughly  a  factor  of  10. 
This  implies  strongly  inhibited  diffusion  for  the  confined  system. 

This  survey  of  the  variety  of  influences  caused  by  the  surface-molecule 
interaction  illustrates  the  strong  motivation  for  the  intensive  investigation  of  this 
type  of  system. 

2.2.    The  Liquid-Solid  Transition 

There  have  been  many  attempts  to  explain  and  describe  the  liquid-solid 
transition.  The  early  phenomenological  theories  relate  critical  behavior  of  solid 
state  properties  to  the  onset  of  melting.  Such  theories  are  characterized  by  the 
increase  in  vibrational  amplitudes,26  the  disappearence  of  bulk  moduli27  and  the 
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proliferation  of  defects  (dislocations)28-29  upon  occurrence  of  the  liquid  state. 
These  are  models  for  the  mechanical,  rather  than  the  thermodynamic,  instability 
of  the  crystal.  The  study  of  a  thermodynamic  criterion,  the  minimum  chemical 
potential,  has  been  pursued  by  molecular  dynamics  studies.30  They  predict 
melting  which  takes  place  in  a  layer-wise  fashion  starting  from  the  surface. 
Caution  towards  the  results  of  such  calculations  is  necessary,  considering  the 
limited  number  of  particles  involved. 

2.2.1.  Elastic  Instability  Theory 

The  original  Born  criterion  of  vanishing  bulk  moduli  was  revised,31 
remedying  the  erroneous  second  order  and  homogeneous  nature  of  the 
transition.  The  difference  between  the  two  theories  hinges  on  the  experimental 
observation  that  one  of  the  bulk  moduli  does  not  completely  disappear  at  the 
transition  temperature  but  only  inside  the  melt.  This  implies  a  two-phase  model 
where  the  dilatation  on  the  liquid  side  is  8'm  which  is  different  from  S*m  on  the 
solid  side  of  the  transition.  The  shear  modulus  vanishes  or  nearly  vanishes 
when  8  =  8'm.  At  the  transition  the  entropy  2RT  is  released  upon  entering  the 
liquid  state.  This  reduces  the  free  energy  and  is  the  reason  why  there  is  another 
dilatation  in  the  solid  state  which  is  associated  with  the  same  free  energy.  This 
point  is  the  true  transition  where  the  isothermal  and  discontinuous 
transformation,  involving  a  latent  heat,  occurs.  The  released  entropy  also 
represents  the  nucleation  barrier  which  renders  this  a  heterogeneous  theory,  as 
the  solidification  occurs  at  strained  sites  and  surfaces  where  the  free  energy  is 
higher. 

2.2.2.  Dislocation  Theory 

A  more  microscopic  view  of  crystallization  was  taken  by  Cotterill29  et  al.  by 
using  a  model  incorporating  dislocations.  The  decisive  assumption  is  that 
dislocations  also  exist  in  the  liquid.  The  dislocation  density  is  extremely  high  and 
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there  is  no  crystalline  material  around  the  mutually  contiguous  dislocation  cores. 
This  makes  it  imprecise  to  describe  the  liquid  as  a  crystal  saturated  with 
dislocations  cores.  Instead,  the  liquid  can  be  referred  to  as  a  pseudodislocation 
saturation.  Viscosity  and  other  liquid  properties  are  explainable  by  the  rapid 
motion  of  pseudodislocations.  Crystallization  occurs  upon  cooling  by  elimination 
of  pseudodislocations  without  recreation  when  the  temperature  is  below  melting. 
The  elimination  of  pseudodislocations  is  driven  by  an  increase  in  the  free  energy, 
as  the  liquid  state  exists  below  its  equilibrium  temperature.  This  phenomenon  is 
called  supercooling  and  indicates  a  metastable  state.  Nucleation  theory  (as 
discussed  later  in  this  chapter)  applies  to  this  system  in  a  standard  way.  If  the 
viscosity  of  the  melt  is  high  or  the  system  is  cooled  rapidly  into  the  solid 
temperature  region  the  removal  of  dislocations  is  inhibited.  The  dislocation 
annihilation  rate  is  then  slower  than  the  solidification  process.  This  solid  is 
microscopically  similar  to  a  liquid  incorporating  a  high  concentration  of 
dislocations,  which  are  immobile  and  incapable  of  diffusion  and  interaction.  A 
polycrystal  with  a  grain  size  of  near-atomic  dimension  develops.  Such  a  solid  is 
called  a  "glass"  and  follows  in  a  straightforward  fashion  from  this  model. 

A  few  properties  of  glasses  shall  be  mentioned  in  this  context.  The 
ground  state  of  a  regular  solid  is  a  single  crystal  which  minimizes  the  dislocation 
energy.  A  glass,  on  the  other  hand,  is  one  realization  of  a  metastable  state 
frozen  in  at  some  arbitrary  configuration.  The  reason  for  its  existence  can  be  its 
frustration  due  to  the  incompatibility  of  the  interactions  with  the  lattice.  The 
specific  volume  of  a  glass  is  larger  than  that  of  a  crystal  which  is  an  immediate 
consequence  of  the  abundance  of  frozen  dislocations  and  the  increased 
disorder.  A  glass  temperature  exists  below  the  bulk  melting  point  from  where  on 
the  mobility  of  dislocations  is  zero,  the  elimination  of  defects  impossible  and  the 
glass  phase  persistent. 
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2.2.3.  Two-Dimensional  Melting  Models 

A  two-dimensional  dislocation  theory  is  especially  interesting  in  the  light  of 
the  two-dimensional  liquid-solid  transition.  Mermin,32  following  Peierls's33  earlier 
arguments,  showed  rigorously  that  conventional  crystalline  long-range  order  is 
excluded  in  two  dimensions  for  power-law  potentials.  This  does  not  pertain  to 
orientational  ordering,  however.  Kosterlitz  and  Thouless34-35  proposed  another 
model  for  the  existence  of  long-range  order.  This  long-range  order  is  of  a 
topological  type.  It  originates  from  the  dislocation  theory  of  melting.36  The  idea 
is  similar  to  the  already  discussed  model  of  Cotterill29  et  al.  The  liquid  phase 
close  to  the  melting  point  has  a  local  structure  similar  to  the  solid  but  dislocations 
exist  which  are  mobile  and  induce  viscous  flow.  The  solid  state  is  rigid  because 
the  dislocations  are  immobile.  They  pair  up  in  the  solid  state  thereby  reducing 
their  energy  from  a  logarithmically  diverging  to  a  finite  contribution.  The  phase 
transition  from  the  solid  to  the  liquid  state  is  characterized  by  the  dissociation  of 
pairs  into  single  dislocations.  The  critical  temperature  is  calculated  to  be 

=  //a2(1+T) 
°  *nkB 

where  ^  is  the  shear  modulus  (in  two  dimensions)  and  r  Poisson's  ratio.  The 
connection  with  the  revised  Born  model31  described  earlier  in  this  context  is  of 
interest.  From  this  theory  the  shear  modulus  cannot  be  zero  in  the  transition 
region,  otherwise  the  transition  to  a  solid  would  occur  at  zero  temperature. 

Starting  from  the  dislocation  model,  Halperin  and  Nelson37  have  worked 
out  a  more  detailed  theory  of  two-dimensional  melting.  They  assume  the 
unbinding  of  dislocation  pairs  which  are  referred  to  as  dipoles  because  they  pair 
up  with  two  different  polarities.   This  yields  an  exponential  decrease  of  the 
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translational  order  parameter.  Instead  of  forming  an  isotropic  liquid  at  this  point, 
calculations  show  a  six- fold  anisotropy  with  the  bond-angle  correlation  decaying 
algebraically.  At  a  higher  temperature  finally,  the  dissociation  of  these 
disclination  dipole  pairs  drives  another  transition  to  an  isotropic  liquid.  The 
intermediate  liquid-crystal  state  is  possibly  bypassed  or  occurs  simultaneously  in 
the  presence  of  a  substrate  due  to  its  orientational  bias.  This  would  transform 
the  two  continuous  phase  transitions  into  one  first-order  transition. 

This  idea  of  a  virtual  region,  where  two  transitions  are  hidden  in  the 
experimentally  observable  first-order  transition,  was  adopted  by  Tallon.38  In  a 
fashion  similar  to  the  revised  Born  theory,  an  entropy  term  is  included.  This 

AV 

could  be  a  volume-independent  term39-40  of  AS  =  R\n2  as  — a  -»  0  which  is  the 

case  in  a  porous  material  where  the  available  space  is  constant  as  a  geometric 
restriction  and  the  molecules  are  trapped  in  the  pores.  This  entropy  term  is 
possibly  associated  with  the  transition  from  a  uniform  distribution  of  the  crystal 
cell  size  to  a  random  one.  Another  possibility  is  the  communal  entropy  arising 
from  the  fact  that  particles  in  the  liquid  state  have  access  to  the  whole  fluid 
instead  of  being  localized.  This  contribution  is  significant  for  large  diffusivities 
but  small  for  a  tortuous  geometry  such  as  zeolite.  The  transition  becomes  first 
order  by  going  discontinuously  from  below  the  dissociation  to  above  the 
disclination  of  dislocation  dipoles  for  two  and  three  dimensions. 

This  leads  to  the  conclusion  that  melting  in  two  and  three  dimensions  is 
likely  to  be  very  similar.38   The  difference  would  be  the  varying  width  and 
separation  of  the  two  virtual  transitions,  both  being  larger  in  two  dimensions. 
2.2.4.  Molecular  Dynamics  Studies 

Studies  from  a  different  starting  point  are  performed  using  the  molecular 
dynamics  approach.    In  this  case  thermodynamic  quantities  such  as  the 
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chemical  potential  are  investigated.  Of  special  interest  is  Broughton30  et  al.'s 
conclusion  that  melting  occurs  in  a  layer-wise  fashion  starting  from  the  outermost 
crystalline  layer  and  progressing  into  the  bulk  region.  The  mechanism  involves 
the  promotion  of  particles  from  the  solid  to  the  liquid  layer  with  an  accompanied 
increase  in  mobility  in  the  remaining  solid  layer.  The  empirical  rule  has  been 
established  that  each  layer  melts  when  it  reaches  the  bulk  ratio  of  kinetic  to 
potential  energy  at  the  melting  point.  Using  this  criterion,  the  outermost  layer 
melts  at  about  72%  of  the  bulk  melting  temperature,  the  second  layer  at  90%, 
while  already  the  third  layer  is  close  to  the  bulk  melting  point. 
2.2.5.  Implications  for  the  Investigated  System 

The  previous  remarks  are  of  fundamental  importance  for  a  system  such 
as  hydrogen  in  zeolite  where  the  purely  geometrical  surface-to-volume  ratio  is  as 
much  as  1/4  and  about  80%  of  all  hydrogens  have  contact  with  the  substrate.  It 
is  therefore  difficult  to  decide  upon  the  dimensionality  of  the  adsorbed  system. 
The  prediction  of  a  similar  melting  process  for  two  and  three  dimensions 
eliminates  some  of  these  problems. 

Implications  for  the  nature  of  the  transition  are  also  given  by  the  existence 
of  hysteresis  upon  heating  and  cooling.  The  absence  of  hysteresis  would  lead  to 
the  conclusion  that  the  transition  is  continuous  and  the  melting  domain  not  a 
heterophase  region. 

Furthermore,  it  has  been  argued41  that  in  the  case  of  grain  boundaries  the 
transition  should  always  be  first  order.  As  will  be  seen  later,  grain  boundaries 
may  have  substantial  impact  on  the  solidification  process. 

The  above  theories  and  projections  are  a  rich  testing  ground  for 
experimental  investigations.  A  variety  of  work  is  necessary  in  order  to 
understand  the  complicated  process  of  melting  and  solidification.  The  outcome 
will  certainly  be  severely  affected  by  the  impacts  of  dimensionality  such  as  the 
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reduced  number  of  nearest  neighbors  and  the  substrate  character  with  an 
interaction  potential  that  leads  to  limited  diffusion. 

2.3.    The  Superfluid 

Superfluidity  for  helium  is  a  well-known  phenomenon  in  the  low- 
temperature  community.  The  important  issue  is  that  the  freezing  point  of  helium 
is  below  the  transition  to  the  superfluid.  This  is  the  case  because  helium  does 
not  solidify  at  any  temperature  without  pressure.  Therefore  the  superfluid  state 
is  feasible  and  transforms  by  Bose  condensation  from  a  normal  liquid  at 
7A=2.17K. 

Superfluidity  is  expected  to  occur  for  any  system  of  bosons.  Interest 
exists  therefore  in  a  potential  superfluid  hydrogen  phase.  Calculations  have 
shown  for  an  ideal  Bose  gas  that  the  transition  temperature  to  a  superfluid42  is 
characterized  by 

*2   /  „\2/3  (  11  V5/3 
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where  M  is  the  mass  of  the  molecule,  Mp  the  proton  mass,  n  the  concentration, 
p  the  density  and  g  the  degeneracy  for  a  single-particle  state  which  is  unity  for 
parahydrogen.  This  results  in  a  predicted  transition  temperature  of  =  3  K  for 
helium  compared  to  the  experimental  value  of  Tl=2.17K.  It  is  in  good 
agreement  and  can  therefore  be  applied  to  parahydrogen  with  a  reasonable 
chance  for  a  realistic  prediction  of  the  superfluid  transition  temperature.  The 
calculation  for  hydrogen  yields  =  6.6  K.  This  low  value  for  the  transition 
temperature  is  the  reason  why  superfluid  hydrogen  has  not  been  observed 
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experimentally:  hydrogen  freezes  at  13.8  K  before  it  reaches  the  superfluid 
transition.  It  is  expected,  however,  that  hydrogen  enters  the  superfluid  state  if 
the  liquid  is  supercooled  to  below  the  superfluid  transition  temperature.  The 
quality  of  the  above  estimate  for  the  superfluid  temperature  is  stressed  by  the 
fact  that  the  transition  temperature  increases43  with  a  larger  interaction  potential. 
The  energy  scale  of  the  Lennard-Jones  potential  between  hydrogen  molecules  is 
e=  37  K  which  is  much  larger  than  e=  10.2  K  between  helium  atoms.  The 
calculation  is  therefore  likely  to  be  a  better  estimate  for  hydrogen  than  for 
helium.  Other  authors44  are  less  confident  about  the  value  of  7"v  They  scale 
the  difference  between  TK  and  with  the  de  Boer  parameter  which  for 
hydrogen  is  only  2/3  of  the  value  for  helium.  The  de  Boer  parameter  describes 
the  degree  of  "quantumness"  which  leads  the  authors  to  believe  in  a  reduced 
experimental  transition  temperature.  On  the  other  hand,  it  seems  that  a  larger 
ZPM  would  rather  be  less  favorable  for  a  collective  process,  such  as  Bose 
condensation. 

In  summary,  superfluid  hydrogen  can  be  expected  below  6  K  but  the  liquid 
state  must  be  extended  down  to  the  superfluid  transition  temperature  which  is 
possible  when  the  liquid  is  kept  metastable. 

Properties  of  the  superfluid  have  been  projected  partly  in  analogy  to 
helium.  The  dispersion  curve  should  also  exhibit  phonon  and  roton  excitations. 
The  slope  of  the  phonon  curve  is  given  by  the  sound  velocity.  The  roton  gap  A 
is  estimated  by  scaling  the  minimum  roton  energy  of  helium  by  the  ratio  of  lattice 
spacings.  This  results  in  a  gap  A  of  around  3kTx.  This  would  lead  to  a  roton 
gap  of  20  K.  A  recent  experiment45  of  hydrogen  in  vycor  suggests  a  tentative 
explanation  of  heat  capacity  data  by  a  roton  gap  of  23  K.  This  occurs  in  addition 
with  supercooling  at  10  K  which  seems  a  very  high  temperature  in  the  light  of 
superfluidity. 
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So  far  the  discussion  only  pertained  to  parahydrogen.  For  orthohydrogen 
with  a  single-particle  degeneracy  of  g  =  9  for  spin  and  orbital  angular  momentum 
degeneracy  the  superfluid  transition  drops  to  1.5  K.  For  this  system  further 
magnetic  properties  could  exist  which  would  generate  additional  diversity  in  the 
phase  diagram. 

The  question  arises  if  supercooling  by  about  8  K  is  a  feasible  operation. 
Several  options  and  speculations  have  been  proposed.  It  has  been  suggested43 
to  investigate  surface  films  on  different  smooth  substrates  where  hydrogen  may 
not  wet  and  stay  liquid.  This  could  be  done  on  top  of  a  monolayer  of  substrate- 
coating  material  (deuterium  or  neon)  generating  a  less  dense  fluid.  The  setup 
would  lead  to  a  two-dimensional  superfluid.  It  has  been  shown  by  Widom46  that 
superfluidity,  excluding  intermolecular  interactions,  is  indeed  possible  in  two-  and 
even  one-dimensional  systems.  The  condition  for  the  existence  is  the  presence 
of  an  external  disturbance,  for  example  gravity  or  rotation,  which  breaks 
translational  invariance.  This  does  therefore  not  contradict  Hohenberg's47 
argument  which  rules  out  lower  dimensional  Bose  condensation  in 
homogeneous  systems,  i.e.  without  symmetry  breaking. 

As  a  confirmation  of  this  theory  the  possibility  of  surface  superfluidity48  in 
parahydrogen  has  been  put  forward  as  a  potential  explanation  of  a  gravity- 
induced  flow  below  3  K  observed  by  Alekseeva  and  Krupski.49  The  calculated 
flow  rates  are  in  agreement  with  experiment  and  result  from  vortex 
considerations  analogous  to  the  Kosterlitz-Thouless  theory. 

Inclusion  of  impurities  is  another  proposed  method  of  achieving 
supercooling  as  suggested  by  Geilikman50  who  extended  the  Lindeman  criterion 
for  melting  to  quantum  crystals  where  the  melting  temperature  Tm  is  much 
smaller  than  the  Debye  temperature  0D.    He  stresses  that  the  zero-point 
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amplitude  is  inversely  related  to  Tm  which  is  in  turn  proportional  to  (M@D)  . 
Impurities  reduce  the  Debye  temperature  if  they  are  heavy,  weakly  bound  to  the 
matrix  and  have  a  high  solubility.  By  these  means  Tm  is  reduced. 

An  interesting  possibility  at  low  temperatures  is  the  creation  of  quantum 
excitations  from  light  impurities,  so  called  impuritons.  Atomic  hydrogen  can  be 
used  which  also  obeys  Bose  statistics.  If  the  impuriton-impuriton  scattering  is 
repulsive  then  the  energy  spectrum  displays  an  acoustic  branch  which  satisfies 
the  superfluidity  condition.51  The  calculated  superfluid  transition  temperature 
is  raised  to  12  K. 

Another  alternative  is  dynamic  supercooling  which  relies  on  a  rapid 
decrease  in  temperature  and  conserves  the  liquid  state  in  a  metastable  form. 
Many  experiments  of  this  kind  have  been  performed,  being  described  in  the  next 
section.  An  explanation  using  nucleation  theory  is  possible  to  demonstrate  this 
effect  of  delayed  and  inhibited  nucleation  to  the  crystalline  state.  The  theory  will 
also  be  applicable  to  the  understanding  of  supercooling  in  small  pores  which  is 
performed  in  the  experiment  to  be  discussed  in  this  dissertation.  The  procedure 
can  substantially  decrease  the  melting  point,  by  as  much  as  5.5  degrees  in  the 
present  experiment.  Unfortunately,  the  lambda  transition  and  the  Fermi 
temperature  are  also  depressed  in  small  geometries  as  experimentally  shown  for 
4He52  and  3He.M  This  makes  it  in  turn  more  difficult  to  reach  superfluidity. 

2.4.  Supercooling 

2.4.1  Nucleation  Theory 

Nucleation  theory54  is  one  possibility  for  understanding  the  preference  of 
the  liquid  versus  the  solid  state.  It  is  a  dynamic  description  why  the  process  of 
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solidification  is  inhibited  and  therefore  the  metastable  liquid  state  preferred.  It 
builds  upon  the  formation  of  nucleation  cores  around  which  solidification  occurs 
and  describes  the  velocity  of  growth,  once  the  nucleation  cores  have  formed. 

The  starting  point  is  the  free  energy  difference  5F  between  liquid  and 
solid  state  required  to  form  a  solid  sphere  of  radius  R 

5F{R)  =  4nR2aLS-±nR2ns{fL-fs) 

where  aLS  is  the  surface  energy55  between  solid  and  liquid,  ns  the  number 
density  in  the  solid  phase  and  fL,  /sthe  free  energies  per  molecule  in  the  liquid 
and  solid  phase.  The  surface  and  volume  energy  contributions  to  the  free 
energy  are  competing.  There  exists  a  maximum  energy 

0     3  n\{fL-fLf 

which  occurs  at  a  radius 

■  a  -  2ocls 
°  ns(fL-fs) 

Figure  1  shows  a  plot  of  the  free  energy  difference  versus  the  solid  sphere  size 
for  different  temperatures.  For  R  >  R0  the  solid  sphere  grows  without  limit 
because  the  cubic  term  in  R  takes  over  and  decreases  the  free  energy.  If  a 
thermal  fluctuation  can  produce  a  solid  sphere  with  a  radius  greater  than  R0,  the 
supercooled  liquid  makes  a  transition  into  the  solid  state.  At  the  normal  bulk 
transition  temperature  of  13.8  K,  fL  and  fs  are  equal.  Therefore,  R0,  the  critical 
radius,  and  the  energy  maximum  for  SF,  which  is  the  solidification  barrier,  both 
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Figure  1.  Nucleation  Theory 
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diverge  implying  that  for  any  lower  temperature  a  large  but  finite  size  solid 
sphere  suffices  to  trigger  the  transition  to  the  solid.  At  this  point  the  liquid  state 
still  prevails,  however.  As  the  temperature  decreases  further,  the  barrier  and  the 
radius  also  decrease.  This  implies  that  a  finite,  thermally  induced  solid  sphere  is 
enough  to  induce  solidification  because,  once  the  nucleation  core  reaches  the 
critical  size  and  the  top  of  the  barrier,  it  is  energetically  favorable  to  grow  without 
bounds  and  release  energy.  The  maximum  radius  for  the  solid  sphere  is  the 
above  defined  radius  RQ.  This  is  in  contrast  to  the  behavior  in  a  porous  material 
where  the  critical  solid  sphere  is  restricted  to  the  pore  size.  The  slightly  different 
condition  will  be  described  in  the  next  section.  For  low  temperatures  the  critical 
radius  and  the  solidification  barrier,  become  small,  indicating  facilitated 
solidification. 

The  probability  for  solidification  to  occur  per  unit  volume  is  described  by 

TT=T°exp(-5F0/kT) 

with  an  attempt  frequency  per  unit  volume 

r?=(nLkT/h)exp{-®LS/kT) 

where  nL  is  the  number  of  molecules  per  unit  volume  in  the  liquid  phase  and  <DLS 
the  activation  energy  for  diffusion  of  atoms  across  the  liquid-solid  interface.  The 
self-diffusion  activation  energy  is  45  K  in  liquid  hydrogen  and  can  be  taken  as  an 
approximation  for  Ols. 

An  understanding  for  the  shape  of  the  solidification  probability  can  be 
gained  in  the  following  fashion.  The  nucleation  rate  increases  rapidly  with  a 
decrease  in  temperature  due  to  the  increased  free  energy  difference  between 
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solid  and  liquid  at  these  supercooled  temperatures.  This  implies  a  decreased 
energy  barrier.  At  low  temperatures,  however,  the  energy  barrier  becomes 
constant  but  the  lower  temperature  reduces  the  probability  of  overcoming  this 
barrier  by  means  of  thermal  fluctuations. 

Nucleation  by  quantum  tunneling  through  the  energy  barrier  has  also 
been  considered.  This  could  be  of  importance  for  quantum  crystals  like 
hydrogen  and  introduces  an  additional  mechanism  which  is  especially  active  at 
low  temperatures,  therefore  increasing  the  nucleation  rate. 

A  very  important  parameter  for  calculating  realistic  numbers  in  this  context 
is  aLS.  The  larger  this  value,  the  more  efficiently  a  liquid  can  be  supercooled. 
Unfortunately,  this  parameter  is  not  known  for  any  material  with  high  accuracy 
but  has  strong  impact  on  the  outcome  because  it  is  contained  in  the  exponent  to 
the  third  power. 

Another  supercooling  model56  which  follows  classical  nucleation  theory  in 
three  dimensions  relies  on  a  scaling  argument.  It  suggests  that  supercooling  of 
hydrogen  takes  place  at  a  reduced  temperature  of  0.8  times  the  bulk  melting 
temperature.  This  leads  to  a  higher  supercooling  temperature  of  1 1  K.  Another 
statement  is  made  from  scaling.  The  involved  degrees  of  freedom  can  be 
calculated  as  a  ratio  of  the  self-diffusion  activation  energy  and  the  bulk  melting 
point.  This  ratio  is  three  and  implies  participation  of  three  translational  and  three 
vibrational  degrees  of  freedom  to  the  process  of  diffusion  across  the  phase 
boundary. 

2.4.1.1.        Implications  for  experiment 

Nucleation  theory  can  be  used  to  devise  experiments  which  could  be  able 
to  achieve  supercooling. 
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The  first  set  of  experiments  exploits  the  dynamical  aspect  of  the  theory.  If 
hydrogen  is  cooled  rapidly  through  the  maximum  nucleation  rate  into  the 
relatively  stable  regime  at  low  temperatures,  substantial  supercooling,  maybe 
down  to  zero  degrees,  and  superfluidity  can  be  achieved. 

Experiments  of  this  kind  have  been  carried  out  at  Brown  University,57-58 
using  liquid  drops.  They  employed  a  helium  pressure  to  levitate  the  injected 
droplets  and  have  indeed  observed  supercooling  by  several  degrees.  Container 
walls  and  impurities  act  as  nucleation  cores  and  had  to  be  avoided. 

The  second  possibility  is  using  confining  geometries  into  which  hydrogen 
is  adsorbed.  This  introduces  a  cutoff  radius  Rc  into  the  system  which  is  related 
to  the  pore  size.  Even  when  a  thermal  fluctuation  is  able  to  create  a  solid  sphere 
of  radius  R0,  from  where  on  it  would  grow  without  limit  in  the  bulk,  it  is  now 
restricted  to  the  cutoff  radius  and  must  stop  growing.  As  long  as  the  liquid  state 
is  energetically  favored,  which  is  the  case  for  positive  5F,  the  system  will  stay 
liquid.  The  cutoff  radius  Rc,  is  therefore  determined  for  SF  =  0  to  be 

=  _3aLi_ 
°  ns(fL-fs) 

It  is  obvious  that  a  smaller  pore  size  results  in  stronger  supercooling,  or  better,  a 
more  depressed  melting  point.  This  process  is  not  dynamically  motivated  and 
relies  entirely  on  the  energetics. 

Experiments  following  this  reasoning  have  been  performed  mostly  using 
vycor,  a  porous  glass,  and  helium52-59-60  61  or  hydrogen.45  62  63 

There  is  an  effect  for  helium  which  is  equivalent  to  supercooling.  Helium 
already  requires  overpressure  to  achieve  bulk  solidification.  The  result  for 
helium  in  restricted  geometries  is  that  an  extra  pressure  of  several  atmospheres, 
in  excess  of  the  bulk  pressure,  is  needed  for  solidification.  This  pressure  is  a 
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function  of  pore  size  and  can  be  deduced  from  the  critical  radius  of  nucleation 
theory  treated  above: 

_    3aLS    =3aLS  Vs 
c   n,(fL-fs)     Ap  VL-V5 

This  determines  the  pressure  difference  as  a  function  of  the  pore  size.  The 
above  equation  has  been  quoted  in  the  literature52-58-59  with  a  factor  of  2  instead 
of  3  which  does  not  take  the  additional  effect  of  the  restriction  by  the  cutoff 
radius  into  account. 

The  same  considerations  can  be  made  for  hydrogen,  now  using  the  latent 
heat  term  but  neglecting  the  entropy  contribution  to  the  solidification  process. 
This  leads  to 

_  3aLS  vsTm 

where  Tm  is  the  bulk  melting  point  and  LLS  is  the  latent  heat  of  fusion  per 
molecule  which  determines  the  amount  of  supercooling  ST. 

This  model  attributes  the  extra  pressure  and  the  depression  of  the  melting 
point  to  the  additional  contribution  of  surface  energy  between  solid  and  liquid. 
No  account  is  made  for  the  failure  of  the  solid  to  wet  the  surface.  The  model  in 
the  next  paragraph  will  attempt  to  explain  this  phenomenon. 
2.4.2.  Grain  Boundaries 

A  microscopic  model64  attempts  to  explain  the  inhibition  of  the  formation 
of  solid  on  the  wetted  surface  and  the  following  rapid  nucleation.  This  model 
originated  from  the  experimental  observation65  that  the  solid-wall  interfacial 
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tension55  asw,  also  termed  surface  energy,  was  found  to  be  greater  than  the 
liquid-wall  tension  orLW  for  helium.  This  result  was  obtained  by  measuring  a 
contact  angle  0  between  solid  and  liquid  on  the  surface  of  a  copper  container 
that  was  larger  than  90  degrees  and  then  compared  with  the  equation 

aLScos0  =  aLW  -asw 

This  was  confirmed  on  a  glass  surface  and  disagrees  wuh  the  theory  of 
Landau  and  Saam.66  Their  theory  states  that  the  dense,  absorbed  film  on  a 
substrate  in  contact  with  the  liquid  is  simply  ordinary  solid  which  is  at  an 
increased  pressure  as  a  result  of  the  strong  Van  der  Waals  attraction  of  the 
substrate. 

The  increased  solid-wall  tension  is  surprising,  as  virtually  all  uniform 
substrates  favor  the  solid  state,  independent  of  their  attraction  potential.  The 
reason  is  that  the  attractive  forces  interact  more  strongly  with  the  denser  solid 
than  with  the  liquid.  The  enhancement  of  the  solid-wall  tension  therefore 
demands  a  microscopic  explanation. 

For  the  explanation  a  microscopically  rough  substrate  with  Van  der  Waals 
interaction  is  considered.  The  adsorbed  material  forms  a  monolayer  solid  on  the 
substrate  walls  due  to  the  attractive  potential.  The  nature  of  this  solid  is  crucial 
for  the  understanding  of  the  reduced  affinity  for  the  solid  state  in  the  bulk  pore. 
The  solid  conforms  to  the  rough  surface  wall  and  is  therefore  irregular  and 
microcrystalline.  Any  growth  of  a  bulk  solid  is  highly  polycrystalline.  This  would 
give  rise  to  a  large  solid  contact  surface  area  between  the  individual  solid 
crystals  and  create  many  grain  boundaries.  They  in  turn  have  a  high  grain 
boundary  energy  which  contributes  to  the  effective  interfacial  energy  between 
solid  and  liquid  and  disfavors  the  solid  state.  The  pore  environment  rs  especially 
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prone  for  such  a  behavior  with  its  large  surface-to-volume  ratio  where  the 
interference  between  crystallites  is  of  the  highest  degree.  The  film  on  the 
substrate  should  appear  as  a  foreign  wall,  not  belonging  to  the  rest  of  the  liquid 
system.  This  is  called  class-ll  growth  and  describes  the  imperfect  wetting  of  a 
film  by  the  condensed  bulk  of  the  same  material.  This  also  reconciles  the 
observation  of  a  large  contact  angle  on  a  substrate  surface  and  a  monolayer  of 
solid  directly  adjacent  to  the  wall. 

The  increased  pressure  for  solidification  of  helium  has  been  calculated 
using  the  increased  chemical  potential  associated  with  the  grain  boundary 
energy.  The  excess  average  energy  per  unit  area  for  large-angle  boundaries 
is67  E  =  0.05/za,  where  a  is  the  lattice  spacing  and  n  the  shear  modulus.  This  is 
multiplied  by  the  volume-to-surface  ratio  /  =  V/A  of  the  porous  material  and 
equated  with  the  pressure  term  of  the  energy  which  leads  to 

where  a  is  a  constant  depending  on  the  roughness  of  the  surface.  It  is  about 
equal  to  the  number  of  particles  in  contact  with  the  walls  of  the  pore  if  the 
surface  nucleates  a  crystal  at  every  particle  site.  The  roughness  parameter  has 
been  given  between  10  by  Shimoda60  et  al.  and  40  by  Dash.64  The  equation 
leads  to  the  approximately  right  overpressure  for  the  solidification. 

In  the  case  of  hydrogen  one  can  again  substitute  the  pressure  energy 
term  versus  the  latent  heat  and  arrives  at 
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using  the  same  notation  all  along. 

There  is  a  discussion  in  place  about  the  two  different  models  proposed  in 
the  discussion  of  supercooling.  Both  models  predict  a  reasonable  value  for  the 
overpressure  required  to  freeze  helium.  There  is  a  clear  distinction,  however, 
between  the  two  physical  interpretations.  In  the  grain  boundary  model  a  thin 
layer  of  adsorbate  wets  the  surface  and  solidifies  on  it.  This  is  not  the  case  for 
the  solid-liquid  nucleation  model  where  a  density  fluctuation  within  the  non- 
wetting  fluid  can  become  stable  and  solidify  under  the  appropriate 
thermodynamic  conditions. 

2.5.    Ortho-Para  Conversion  in  Molecular  Hydrogen 

One  of  the  most  fascinating  properties  of  molecular  hydrogen  is  the 
explicitly  quantum-mechanical  classification  into  two  species.  The  quantum- 
mechanical  requirement  of  a  totally  antisymmetric  wavefunction  for  two  fermions 
can  be  realized  in  two  distinct  ways,  leading  to  two  distinguishable  molecular 
species:  ortho-  and  parahydrogen.  Orthohydrogen  has  a  symmetric  nuclear 
spin  wavefunction  (/=  1)  and  an  antisymmetric  orbital  wavefunction  (J  odd).  The 
opposite  is  true  for  parahydrogen. 

The  ortho-para  ratio  is  determined  by  the  availability  of  states  given  by  the 
partition  function.  In  the  high-temperature  limit  the  even  and  odd  sums  over  J  in 
the  partition  function  yield  the  same  result  so  that  the  spin  parts  lead  to  a  3:1 
ratio  of  ortho-  to  parahydrogen  which  is  the  relative  weight  of  the  nuclear  spin 
degeneracy  for  the  ortho  and  para  species.  Because  of  the  large  separation  in 
the  rotational  energy  levels,  the  equilibrium  concentration  of  orthohydrogen 
drops  from  75%  to  minute  amounts  at  temperatures  below  about  10  K.  At  low 
temperatures  essentially  only  the  parahydrogen  groundstate  and  the  first  excited 
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state  (ortho,  J  =  1)  are  occupied  as  a  result  of  the  large  excitation  energy  in 
comparison  to  the  thermal  energy.  The  conversion  process  from  ortho  to  para 
molecules  requires,  however,  the  simultaneous  breaking  of  the  spin  and  orbital 
symmetry,  and  the  equilibrium  ratio  is  therefore  established  only  very  slowly.  For 
the  isolated  molecule  the  conversion  is  absolutely  forbidden.  In  the  condensed 
phase  the  rate  is  determined  by  the  intermolecular  magnetic  dipole-dipole 
interaction  and  magnetic  field  gradients  caused  by  magnetic  impurities. 
Because  of  this,  ortho-  and  parahydrogen  can  be  treated  as  two  different 
molecular  species. 

The  first  theoretical  calculation  of  the  ortho-para  conversion  in  bulk  solid 
hydrogen  was  performed  by  Motizuki  and  Nagamiya68  and  more  recently 
revisited  by  Berlinsky  and  Hardy,69  and  Berlinsky.70  The  two  main  interactions 
responsible  for  the  conversion  are:  1)  the  nuclear  spin  dipole-dipole  interaction 
between  two  ortho  molecules  and  2)  the  interaction  of  the  proton  spin  of  one 
molecule  with  the  rotational  magnetic  moment  of  the  other.  It  is  important  to 
note  that  parahydrogen  is  not  part  of  this  process  as  it  does  not  possess  spin  or 
orbital  magnetic  moments.  The  interaction  Hamiltonian  consists  of  two  parts 

Hint  =  Hss  +  Hra 

where  "s"  stands  for  the  spin  and  V  for  the  rotational  contribution.  The  spin-spin 
Hamiltonian  has  four  identical  contributions,  expressing  the  combinatoric  of  the 
two  pairs  of  proton  spins  involved.  The  rotational  interaction  involves  only  two 
contributions  of  the  orbital  momentum  of  one  molecule  with  the  two  spins  on  the 
other  molecule.  The  individual  Hamiltonians  are  of  the  standard  dipole-dipole 
type  (being  described  in  chapter  3.1.)  which  can  either  be  written  in  Cartesian  or 
in  spherical  coordinates  using  spherical  harmonics  with  the  convention  of 
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Rose.71  The  total  wave  function  is  a  product  of  the  orbital  and  spin  part. 
Orthohydrogen  has  3  degenerate  wavefunctions  for  the  orbit,  0Jm,  and  the  spin, 
reflecting  the  quantum  numbers  J  =  1  and  /  =  1  and  their  magnetic  states. 
This  leads  to  nine  degenerate  states.  Parahydrogen  only  has  one  state  O^^. 
This  implies  that  there  are  81  initial  states,  /,  and  9  final  states,  f,  for  a  conversion 
of  two  ortho  molecules  into  one  ortho-  and  one  parahydrogen  molecule.  The 
conversion  rate  flcan  be  calculated  using  Fermi's  golden  rule 

fl  =  ^XW",J/>l2  S(E,-Ef) 
n  if 

where  P,  is  the  probability  of  being  in  one  of  the  initial  states  /.  The  energy 
released  in  the  process  of  conversion  is 

AE  =  E,-Ef=BJl{Jl  +  y-BJf{J,  +  '\)  =  2B  =  '\7J\  K 

and  must  be  taken  up  by  the  lattice.  The  bulk  Debye  temperature1 1-72-73  is  120  K 
and  indicates  that  one  phonon  is  not  enough  to  absorb  the  released  energy. 
Instead,  two  or  three  phonons  are  necessary. 

Avoiding  most  of  the  tedious  and  involved  calculation  of  matrix  elements, 
the  results  for  the  two-phonon  process  are  discussed.  The  rate  R  is  a 
combination  of  the  dominant  spin-spin  and  the  much  weaker  spin-orbital 
contribution: 

4/is2 

where  //2  =  10/i*.  Finally,  RM  is  found  to  be 
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with  a  complicated  geometric  factor  77  =  0.7331,  depending  on  the  Debye 
temperature,  and  the  density  p.  The  interesting  feature  is  the  strong 
dependence  on  the  Debye  temperature.  Any  changes  in  this  quantity  due  to 
confinement  would  have  great  impact  on  the  conversion  rate.  Three-phonon 
processes  are  estimated  to  be  only  a  few  percent  of  the  main  contribution. 

One-phonon  processes  are  not  allowed  under  bulk  conditions,  as 
discussed  earlier.  However,  under  conditions  which  change  the  Debye 
temperature,  like  a  restricted  geometry  or  a  two-dimensional  system,  this  could 
become  important.  Berlinsky  calculated  the  one-phonon  process  and  expressed 
the  rate  as  a  function  of  density  and  complicated  functions  that  have  to  be 
extracted  from  experimental  data.  He  found  that  high-density  hydrogen  mainly 
converts  by  emission  of  single  phonons.  The  rate  for  one-phonon  processes  is 
proportional  to  Q^R'10  compared  to  G^3/?"12  for  two  phonons.74 

The  calculations  lead  to  conversion  rates  of  1.94%/h  by  Motizuki68  et  al. 
and  1.67%/h  by  Berlinsky69  et  al.  which  have  to  be  compared  with  measured 
rates  of  1.75%/h7576  and  1.82%/h.77  This  is  in  excellent  agreement  although 
several  approximations  are  introduced  into  the  theory.  The  theory  is  using  a 
Debye  approximation  with  the  energy  proportional  to  the  absolute  value  of  the 
wave  vector.  This  is  a  poor  approximation  for  the  high  frequency  part  of  the 
phonon  spectrum.  It  also  does  not  include  anharmonicities  due  to  the  large 
ZPM.  The  interactions  are  averaged  over  a  powder,  assuming  random 
orientations  of  crystallites.  The  conversion  process  is  bimolecular  or 
autocatalytic,  described  by 
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x  =  -kx2 

where  x  is  the  orthohydrogen  concentration  and  k  the  conversion  rate  constant. 
The  solution  to  this  differential  equation  is 

X(t)=—?°  — 

w  UxQkt 

where  x0  is  the  initial  ortho  concentration,  75%  in  the  case  of  hydrogen  at  STP. 

If  any  paramagnetic  impurities,  for  example  oxygen,  are  present  there  is 
an  additional  process.78  These  impurities  produce  inhomogeneous  static 
magnetic  fields  around  themselves  to  which  hydrogen  molecules  couple.  The 
calculation  procedure  corresponds  to  the  autocatalytic  process.  The  main 
difference  is  that  the  decay  is  monomolecular  and  described  by 

x  =  -kx 
with  an  exponential  solution 

x{t)  =  x0exp(-kt) 

In  the  case  of  oxygen,79  for  example,  the  conversion  rate  is  more  efficient  within 
a  sphere  of  radius  2.95  R0  than  in  the  autocatalytic  case.  In  this  sphere  mostly 
parahydrogen  is  found.  Measurements  regarding  the  influence  of  paramagnetic 
impurities80  have  been  carried  out  showing  fast  exponential  decay  rates.  The 
important  result  is  a  relationship  between  the  rate  constant  and  the  square  of  the 
magnetic  moment.  This  can  be  used,  in  an  analogy  with  experimental  data,  to 
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estimate  an  upper  bound  on  the  concentration  of  impurities  present  in  a  mainly 
clean  system. 
2.5.1.  Diffusion 

So  far,  diffusion  effects  have  not  been  taken  into  account.  Diffusion 
homogenizes  the  distribution  of  orthohydrogen  around  impurities  and  in 
depletion  zones.  This  establishes  a  perfect  equilibrium  distribution  for 
temperatures  above  8  K.81  As  long  as  this  is  the  case  the  conversion  is 
temperature  independent  because  the  energy  of  phonons  created  by  conversion 
is  much  larger  than  the  thermal  excitation  kT.  Thus  the  occupation  numbers  of 
thermal  phonons  are  extremely  low  and  have  no  impact  on  the  conversion  rate. 
Assuming  a  binomial  distribution  Wn  for  the  equilibrium  distribution  of  ortho 
molecules,  the  mean  number  of  orthohydrogen  molecules  around  each  ortho 
molecule  is  given  by 

M  =  f^nWn=Nx 

n=0 

where  N  is  the  number  of  nearest  neighbors.  The  ortho-para  conversion  leads 
to  a  deviation  from  equilibrium  which  can  be  balanced  by  self-diffusion  which 
tends  to  homogenize  the  distribution.  If  the  diffusion  is  slower,  the  conversion 
rate  is  determined  by  the  diffusion  bottleneck.  In  addition,  clustering  occurs  as  a 
result  of  the  anisotropic  interaction  between  orthohydrogen  molecules  which  is  of 
electric  quadrupole-quadrupole  character.  In  dilute  bulk  samples  orthohydrogen 
molecules  can  diffuse  by  quantum  tunneling  until  they  are  close  to  an  ortho 
neighbor.  The  two  molecules  can  then  form  a  bound  pair  as  a  result  of  their 
anisotropic  interactions  if  the  temperature  is  low  enough  (kgT  <VaniS0  <1.5  K). 
At  low  temperatures  the  overall  conversion  rate  can  therefore  be  accelerated 
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due  to  the  conversion  of  the  clusters.  For  high  temperatures,  thermally  activated 
diffusion  is  clearly  the  dominant  process,  leaving  the  molecules  well  separated, 
while  at  temperatures  around  1.6  K  clustering  becomes  significant. 

This  was  described  by  Schmidt82  for  bulk  solid  hydrogen.  He  found  that 
the  diffusion  is  large  enough  at  4  <  T  <  12  K  to  sustain  the  same  conversion  rate 
(k  =  1.9%/h)  while  at  1.57  K  the  quadrupole  effect  overcompensates  the  reduced 
diffusion  and  leads  to  an  increased  conversion  rate  after  200  hours.  The  rate 
increases  again  after  900  hours  in  the  bulk.  This  result  is  understood  in  terms  of 
the  enhanced  conversion  in  clusters  due  to  the  reduction  in  the  ortho-ortho 
separation  compared  to  a  homogeneous  distribution.  Oxygen  impurities  could 
be  excluded  by  experimental  control  of  oxygen  contents  which  did  not  influence 
the  conversion  substantially.  The  reason  is  that  oxygen  freezes  out  at  higher 
temperatures  and  has  low  solubility  in  hydrogen.  This  may  be  different  for  a 
substrate  where  oxygen  can  adsorb  onto  the  surface  and  act  as  a  catalyst  for 
the  conversion  mechanism. 

In  the  case  of  a  porous  material  the  majority  of  the  molecules  is  in  contact 
with  the  wall  and  only  a  few  nearest  neighbors,  depending  on  the  roughness  of 
the  adsorption  site,  while  the  remainder  is  in  a  bulk-like  surrounding  in  the  center 
of  the  pore  with  12  nearest  neighbors.  The  value  of  the  conversion  rate  is 
therefore  expected  to  be  different  for  a  restricted  geometry  due  to  interactions 
with  the  walls,  the  reduced  number  of  nearest  neighbors  and  the  changes  in  the 
phonon  spectrum. 

The  ortho-para  conversion  rate  in  an  environment  different  from  bulk 
needs  to  be  determined  accurately  in  order  to  characterize  the  hydrogen 
samples  used  in  subsequent  experiments. 


CHAPTER  3 
NUCLEAR  MAGNETIC  RESONANCE:  NMR 


3.1.    Dipolar  Hamiltonian 


Hydrogen  spins  in  a  magnetic  field  experience  the  Zeeman  interaction 
with  the  applied  external  field  as  well  as  inter-  and  intramolecular  dipole-dipole 
interactions.  The  intermolecular  forces  are  weak  compared  to  the  intramolecular 
forces  because  the  interactions  fall  off  with  the  cube  of  the  distance  so  that  they 
contribute  only  a  small  fraction 
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=  2.1% 


where  the  effect  of  the  different  spin  quantum  numbers  for  orthohydrogen  and 
protons  has  been  included.  This  small  term  can  be  neglected  in  the  discussion. 

The  orbital  nuclear  Zeeman  energy  associated  with  the  Hamiltonian 
Hj=hy2H0Jz  and  the  spin-orbit  interaction  Hso  =  hyp\- JHlnt  vanish,  i.  e.  are 
"quenched."1  The  argument  was  discussed  in  the  quadrupolar  section  of 
chapter  1  and  is  related  to  the  nondegenerate  nature  of  the  rotational 
groundstate.  The  dominant  part  of  the  Hamiltonian  is  therefore  the  nuclear  spin 
Zeeman  energy  with  a  perturbation  from  the  dipolar  Hamiltonian83 
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H  =  -yMM  + 1])  +  ^[l1  •  I2  -  3(P  •  n)(l2  •  n)] 


where  H0  is  the  magnetic  field  in  z-direction,  I  is  the  nuclear  spin  vector  with 
superscripts  1  and  2  for  the  two  spins  and  n  =  r12/r  the  unit  radial  vector 
between  spins.  The  dipole-dipole  term  can  be  rewritten  in  spherical  coordinates 
in  order  to  find  a  quantum  mechanical  expression 

HDD  =  ^-{  l1l2-3[/;cos0  +  sin©(/;;cos<I>  +  /;sincD)]x 

[/2cos0  +  sin0(/2cosO  +  /Jsin<I>)] } 

The  angle  0  is  taken  between  the  direction  of  the  field  and  the  intermolecular 
axis.  The  x  and  y  components  of  the  spins  can  be  combined  to  generate  raising 
and  lowering  operators 

HDD  =       {  V •  I2  - 3[l'z  cosQ  +  i sin ©(/>-•*  +  /V*)] x 
[/2cos0  +  isin0(/+2e-1*  +  /fe+i*)] } 

and  can  be  collected  into  six  terms 


HDD  =^$-(A  +  B  +  C  +  D  +  E  +  F) 


which  correspond  to  different  transitions.  The  individual  terms  are 
A  =  /;/2(1-3cos20) 

b = 1(1  -  3cos2  e)(/]/2  - 11  •  i2)  =     -  3cos2  ©)(/;/_2 + rii) 
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C  =  -f  sinGcosee-^/X2 

D  =  C  =-fsin0cos0e+i*(/]/_2  +  /V2) 

E  =  -|sin20e-2l*/X2 

F  =  E*  =-|sin20e+2i*/Vf 

The  reason  for  displaying  the  dipole-dipole  interaction  in  this  fashion  is  that  the 
different  terms  cause  distinct  transitions.  Terms  A  and  B  do  not  change  the 
overall  magnetic  quantum  number,  C  raises  it  by  one,  D  lowers  it  by  one  and  E 
and  F  increase  and  decrease  it  by  two,  respectively.  Choosing  degenerate 
states  |+-)  and  |-+)  with  the  meaning  that  spin  1  is  up  and  spin  2  is  down,  and 
vice  versa,  demonstrates  the  operation  of  each  term.  Term  A  directly  connects 
either  state  and  is  therefore  diagonal.  Term  B  flip-flops  one  state  into  the  other 
and  therefore  connects  the  superposition  of  states,  which  is  the  proper  state  in 
the  first  place.  It  is  therefore  also  diagonal.  Off-diagonal  elements  only 
contribute  by  means  of  a  second  order  correction  to  the  wave  function  which 
results  from  the  perturbation  of  the  dipolar  Hamiltonian.  Terms  C  and  D  are  off- 
diagonal  in  first  order,  E  and  F  in  second  order  and  contribute  therefore  only 
negligibly  to  the  eigenenergies  of  the  dipolar  system.  They  will  be  omitted  from 
here  on. 

3.2.    Continuous  Wave  Lineshape 

The  Hamiltonian  is  now  in  a  form  that  can  readily  be  evaluated  for 
orthohydrogen  with  J  =  1  and  M  ■  +1,  0,  -1.  The  wave  functions  are  |++), 
vs(|+-)  +  |-+))  and  | — )•  Bracketing  the  Hamiltonian  with  these  states 
produces  the  eigenenergies 
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E++=^Ho+M-(1-3cos20) 

^oo  =  -^(1-3cos20) 

E__=-yhH0+       (1  -  3  cos2  0) 
Ar  ' 


Transitions  occur  between  adjacent  energy  levels  due  to  the  selection  rule 
A/W  =  ±1  which  leads  to  two  different  absorption  frequencies  that  are  no  longer 
degenerate  as  they  were  in  the  pure  Zeeman  case 


2 

hv,  =  E_  -  E00  =  yhH0  -  3 cos2  0) 

Ar  ' 

hv2  =  EM  - E++  =  yhH0  -^-(1  -  3cos2  0) 


This  calculation  was  carried  out  classically.  As  shown  in  the  introduction 
the  proton-proton  intramolecular  axis  is  rapidly  rotating,  a  consequence  of  the 
quantum  rotor  property.  An  average  has  to  be  invoked  on  the  angular 
dependence.  The  frequencies  v+  and  v.  are  symmetrically  located  above  and 
below  the  undisturbed  Zeeman  frequency  v0.  The  frequencies  are 

v±=±— tf(3cos20-l\ 

1  %y2 

where  in  the  liquid  case  d  =  -    '     =  57.68  kHz.    In  solid  hydrogen  the 

5  2n{r  ) 

constant  is  slightly  different  with  a  value  of  54.2  kHz.  The  angle  0  is  measured 
between  the  external  z-direction  of  the  magnetic  field  and  the  intramolecular 
axis.  It  is  physically  more  insightful  to  transform  to  another  angle.  This  angle  (5 
defines  the  orientation  of  a  surface  or  the  crystal  field  in  the  local  environment  of 
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an  ortho  molecule  relative  to  the  magnetic  field.  The  average  over  all  space  with 
respect  to  this  direction  yields  the  2/5  factor  for  a  quantum  mechanical  object. 
The  result  is,  in  addition,  also  dependent  on  the  rotational  state  ("shape")  of  the 
molecule  characterized  by  the  normalized  quadrupolar  order  parameter 

This  leads  to  the  final  result 

v±=±|ad(3cos2/?-l) 

which  will  be  discussed  as  it  describes  the  continuous  wave  lineshape. 

For  a  fixed  order  parameter  a  at  a  fixed  angle  fi  two  sharp  lines  occur 
which  are  inhomogeneously  broadened.  This  is  the  case  for  a  single  crystal. 

If  the  angle  /?  is  distributed  from  0  to  90  degrees,  as  is  the  case  of  a 
powder,  the  lineshape  has  the  famous  Pake  doublet  form  (Figure  2).  The 
intensity  corresponding  to  the  frequency  can  be  found  from  a  space 
transformation  argument.  Since  the  dipoles  are  isotropically  distributed  the 
fraction  of  pairs  lying  in  dp  is  d(cos/7).  With  the  lineshape  function  being  g(v) 

and  the  isotropic  space  intensity       =  —  it  turns  out  that 

An 

g(v)dv  =  l{Q)dQ 


*>n  dv  18cfcrcos/7 
2nd{cosp) 
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Figure  2.  Pake  doublet  with  order  parameter  a  =  1 
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Solving  for  cos/?  in  terms  of  frequency  and  order  parameter,  using  the  frequency 
equation,  results  in  the  lineshape 


where  v0  =  .  The  complete  lineshape  is  a  superposition  of  the  positive  and 

negative  branch.  The  "+H  sign  is  valid  for  -v0  <  v+  <+2v0  and  the  "-"  sign  for 
-2v0  <  v_  <  v0. 

The  shape  can  be  understood  when  the  degree  of  space  degeneracy  is 
considered.  The  90  degree  case  is  much  more  probable  because  it  covers  the 
whole  plane  compared  to  one  parallel  alignment  in  the  0  degree  case. 

When  the  constraint  of  constant  cr  is  relaxed,  the  lineshape  becomes  a 
superposition  of  Pake  doublets,  leading  to  a  bell  shaped  curve.  This  behavior  is 
attributed  to  glassy84  85  behavior. 

A  more  detailed  equation  for  the  frequency  shifts  has  been  given86 
incorporating  the  eccentricity  r]  =  (j2x  -J2y  \  of  the  hydrogen  molecule 


V3 


3do 


vt=±-d  oi3cos2/?-1)  +  -r7Sin2/?cos2/? 


In  most  cases  the  assumption  of  zero  eccentricity  holds. 


3.3.    Relaxation  Times 


The  dipolar  Hamiltonian  is  frequently  written  in  tensorial  form 
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2 


HDD  =  £F«M<« 


where  the  FW)  are  random  functions  of  the  relative  positions  of  two  spins  in  time 
and  the  Aw  are  spin  operators.  The  products  just  rewrite  the  terms  A  to  F  from 
before.  Terms  A  and  B  are  equivalent  to  q  =  0,  C  and  D  to  q  =  ±1,  and  E  and  F 
to  q  =  ±2.  The  F{~q)  =  F[q)'  only  contain  the  spherical  coordinates  and  the 
Al~q)  =  A{q)+ ,  which  are  nonhermitian,  consist  of  all  spin  operator  components. 

The  correlation  function  g  measures  the  relationship  between  the  random 
spatial  functions  at  different  times  averaged  over  a  statistical  ensemble 


Their  Fourier  transforms  are  the  spectral  densities  J^ico)  which  describe  the 
interaction  strength  as  a  function  of  frequency 


The  solution  of  the  equation  of  motion  for  the  density  matrix  a  in  the 
interaction  picture  is 


Jqq,(a})  =  jgq,(T)e-^dT 


where  starred  quantities  indicate  the  interaction  representation,  transformed  with 
the  Zeeman  Hamiltonian.   This  is  just  another  terminology  for  the  "rotating 
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frame".  The  observable  q  in  an  experiment  is  averaged  over  the  macroscopic 
collection  of  subsystems  acted  upon  by  the  operator  Q 

q'(t)  =  (Qy=tr{o'(t)Q} 

This  is  a  relationship  that  describes  the  dynamics  due  to  the  dipolar  Hamiltonian 
decoupled  from  the  fast  and  trivial  dynamics  of  the  Zeeman  Hamiltonian. 

For  the  hydrogen  system  the  operator  is  Q  =  l\  +  /*.  Taking  the  trace  over 
the  solution  of  the  equation  of  motion  multiplied  by  Q  results  in 

where 


«,(;)  =  tr^H'DD(t-  r),[H'DD(t),(ll  +  £)]]  <r,(;)j 

which  defines  a*  and  a0  and  introduces  Tv  the  longitudinal  or  spin-lattice 
relaxation  time.  It  describes  the  time  it  takes  for  the  spin  system  to  relax  back  to 
equilibrium  along  the  external  magnetic  field.  Inserting 


into  the  system  of  equations  leads  to 
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dt 


=  -tr 


by  generating  the  spectral  density  J  from  the  random  spatial  functions. 

The  computation  of  the  commutators  leads  to  an  expression  for  the 
relaxation  time  in  the  high  temperature  regime  by  identification  of  terms 


l  =  |r^2/(/  +  1){/1V0)  +  J(2)(2^0)} 


where  co0  is  the  Larmor  frequency. 

A  similar  calculation  leads  to  the  transverse,  or  spin-spin  relaxation  time. 
It  describes  the  time  dependence  of  the  magnetization  amplitude  in  the  x-y 
plane.  The  initial  magnetization  is  assumed  to  lie  along  the  x-axis  with 
expectation  value 

(lx)'=tr{lxo}  =  tr{e->">>lxeiH<'a} 

The  time  evolution  for  the  x-component  of  the  magnetization  to  the  zero  value 
equilibrium  is 


d(ll+ll) 


dt 


=  -tr 


■  o 


The  result  of  the  evaluation  leads  to  the  transverse  relaxation  rate 
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J-=yV/(/  +  1){fJ(0)(0)  +  ^J«((yo)  +  lJ(2)(2Q;0)} 

'2 

The  crucial  difference  between  T:  and  T2  is  that  T2  has  a  contribution  at 
zero  frequency.  For  a  correlation  time  x  much  shorter  than  the  Larmor 
frequency  all  spectral  densities  become  independent  of  frequency  and  equal  to 
Jw(0).  It  is  possible  to  determine  the  ratio  between  the  different  orders  of 
spectral  densities  by  keeping  in  mind  that 


j(o).  jd).  J(2)  =  F(0)  2:F(d  2.  F(2)  2  =  6.t4 


This  implies  that  7,  and  T2  are  equal  in  the  case  of  short  correlation  times  with 
"white"  spectral  densities  for  frequencies  far  above  the  Larmor  frequency. 
3.3.1.  Spectral  Densities 
3.3.1.1.        Rotational  motion 

Functional  relationships  for  the  spectral  densities  are  needed.  They  are 
developed  starting  from  the  correlation  function 


gq(T)  =  F«\t)F«»(t+x) 


In  the  case  of  rotational  motion  the  average  extends  over  space  elements 
Q  and  Q0.  The  average  takes  place  over  a  function  p(Q,f;n0,f0)  which  defines 
the  probability  of  taking  on  the  value  Q  at  r  and  Q0  at  r0.  This  can  be 
alternatively  expressed  as 


p{Q,t;Q0,tQ)=P(n,t;nQ,t0)p(Q1t) 


56 


where  P  is  the  conditional  probability  of  taking  on  Q.  starting  out  with  Q0. 

The  rotation  of  a  molecule  can  be  compared  to  a  sphere  of  radius  b 
rotating  in  a  medium  of  viscosity  77  and  is  described  by  a  diffusion  equation 

at  b2  s 

where  As  is  the  Laplacian  on  the  sphere  and  Ds  is  Stokes'  rotational  diffusion 
coefficient 

3  Qnr]b 

The  solution  to  this  equation  is  the  probability  P(Q,t;Cl0,t0)  that  the  two  spins 
have  the  orientation  Q  at  time  f  after  having  the  orientation  Cl0  initially.  The 
initial  condition  of  the  problem  is 

¥(Q,0)  =  £(Q-Q0) 

and  the  solution  is  sought  for  in  the  form  of  a  spherical  harmonics  expansion 
^(Q,f)  =  ScJMVOw(") 

JM 

noting  that  AsVjM  =  -J(J+1)YjM.  Introducing  this  into  the  diffusion  equation  leads 
to  an  equation  in  the  expansion  coefficients 

drM  n 
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which  describes  an  exponential  relationship  with  time  constant 

b2 

An  expansion  of  the  delta  function  can  be  given  in  terms  of  spherical  harmonics 
which  determines  the  initial  coefficient  and  therefore  the  total  time  behavior  of 
the  expansion  coefficient  Cj(t)  =  YjM(Q0)e~'/T.  This  supplies  the  conditional 
probability 

P(Q,r;Q0,f0)  =  XKHW  (W"* 
needed  to  compute  the  correlation  function 

using  the  isotropic  probability  p(Q0,L)  =  — .  Since  the  random  spatial  functions 

An 

with  distance  r  between  the  two  spins  in  the  molecule  are  known  to  be 

r  r  V  1 5 
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it  implies  that  the  exponential  correlation  functions  are  also  established. 
By  Fourier  transform,  the  spectral  densities  are  given  by 


J(0V)  = 


24 

15r6 

1+o>V 

4 

r 

15/-6 

1+<uV 

16 

r 

J(>) = 


J(2)(<y)  = 


3.3.1.2.        Translational  motion 

The  motion  is  again  characterized  by  the  diffusion  equation.  Its  solution  is 
the  conditional  probability  for  two  particles  to  be  a  distance  r  away  from  each 
other.  The  integration  of  the  correlation  function  involves  spherical  Bessel 
functions  and  has  a  more  complicated  closed  form.  The  limits  are 

Jao))~ — r-r  for<2>r«1 
5  a3 

JAcq)~  -4-7      for  cm »\ 

where  a  is  the  nearest  neighbor  separation.  The  functions  are  proportional  to 
the  density  n  of  particles  in  the  material. 
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3.4.    Motional  Narrowing  and  Second  Moment 

As  already  mentioned,  motion  is  an  important  ingredient  in  understanding 
the  properties  of  hydrogen.  In  the  case  of  NMR,  the  width  and  the  shape  of  the 
absorption  line  are  affected  by  this  motion.  When  spins  are  in  motion  the  field 
experienced  by  a  neighboring  spin  is  fluctuating  on  a  time  scale  xc,  the 
correlation  time.  These  fluctuations  are  averaged  out  over  a  longer  observation 
time.  This  average  is  much  smaller  than  the  instantaneous  local  field,  thus  the 
motional  narrowing. 

There  is  an  intuitive  derivation  of  the  time  scales  involved.  The  line  width 
Acq  is  the  inverse  of  the  time  t  it  takes  for  two  spins  to  dephase  by  an  angle  of 
order  unity.  It  is  assumed  that  the  local  field  difference  Aco0  between  two  spins 
stays  the  same  with  its  sign  changing  at  an  average  frequency  of  1/  rc  »  Ao)0. 
The  dephasing  angle  after  time  rc  is 

SO  =  ±tcAq)0 

and  the  mean  square  angle  of  this  random  walk  motion  is 


AO2  =  /7<53>2  =  nAo)02Tc2 
where  n  =  tl  rc.  After  a  unit  change  of  angle  the  rate  is  given  by 
-  =  Aco  =  Aa)2rc 


60 


which  means  that  the  narrowing  is  more  distinct  for  shorter  correlation  times. 
The  effect  on  7"2  is  directly  related  to  the  correlation  time  as  it  is  the  time  it  takes 
for  spins  to  dephase  by  one  radian 

r  =JL  =  _L_ 

2    Ad)  Aco2rc 

It  also  shows  the  dependence  of  T2  on  the  line  width.  The  value  of  7"2  is  a 
function  of  motion  and  can  be  quantified  as  a  function  of  the  spectral  densities. 
In  the  case  of  rapid  motion  mc  « 1  all  spectral  densities  are  identical  while  for 
slow  motion  mc»'\  the  dominant  contribution  originates  from  J(0)(0).  This 
leads  to  a  ratio  of  (f  x  6  +  x  1  +  j  x  4)  /  (§  x  6)  =  ^  which  states  that  T2  is  larger, 
the  line  width  smaller,  by  a  factor  of  10/3  for  total  motional  narrowing. 

The  effect  of  increased  motion  on  the  lineshape  and  the  second  moment 
is  investigated  in  the  following.  For  this  purpose  the  correlation  function  of  the 
magnetization 


G(t)  =  tr{Mx(t)Mx} 


is  calculated.  The  time-dependent  magnetization  is  a  result  of  its  equation  of 
motion  and  leads  to 


G(t)  =  ^\(E0s\M0\E'0  sf  exp{iay}exp 


\jo)(t')dt' 


where 


<y(0  =  (E0s\HDD(t)\E0s)  -  (E'0  s'\HDD(t')\E'0  s') 
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For  random  motion  the  correlation  function  is  averaged  over  a  frequency  density 
f{a>)  and  collecting  the  prefactor  leads  to  an  expression 

G(f)  =  e^\P{a>)e^rwdco  =  e^le^ 


The  use  of  the  central  limit  theorem  rewrites  the  average 


e^')  =  exp(-i(jaKOrt'f) 


=  expj  -\2  j  (/  -  T)(o)(t)co(t  -  T))dr 


and  finally  leads  to 


G{t)  =  ew  exp 


-{co2)\(t-T)gMdT 


where  g  is  the  reduced  correlation  function  with  the  property  g^O)  =  1.  With  this 
result  the  two  limits  of  short  and  long  correlation  times  can  be  investigated. 

In  the  long  correlation  time  limit  (<y2)i-c2»1  the  reduced  correlation 
function  becomes  unity  and  the  integration  yields 

G(t)  =  e,a°lexp{-i(o)2)t2} 


which  is  Gaussian  and  the  rigid  lattice  form.  Its  Fourier  transform  is  the 
lineshape  and  also  a  Gaussian. 
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In  the  short  correlation  time  limit  {co2)?*  « 1  it  is  permissible  to  neglect  r 
compared  to  t.  This  leads  to  an  exponential  decay 

G(0  =  e^exp{-(^2)rrc} 

whose  Fourier  transform  indicates  a  Lorentzian  lineshape. 

The  conclusion  of  this  analysis  is  that  the  lineshape,  in  addition  to  the  line 
width,  is  an  indicator  of  the  dynamic  regime  the  system  is  in,  whether  motional 
narrowing  takes  place  or  not. 

A  discussion  of  the  moments  of  the  lineshape  follows  naturally.  The  n-th 
moment  is  defined  as 

Mn=j(o)-Q)0)nf(a))dcD 

where  f(oo)  is  the  normalized  lineshape  function  which  leads  to  vanishing  odd 
moments  if  it  is  symmetric  around  the  Larmor  frequency  a)0.  The  previously 
introduced  (<y2)  is  the  second  moment.  Its  form  can  be  Gaussian  or  Lorentzian 
with  different  results  for  the  moments. 
For  a  Gaussian 

the  even  moments  are  M2n  =  (2/7  -  1)A2".  The  half  width  at  half  height  (HWHH)  is 
defined  by 


(<»-<»o) 
2A2 


f(a>Q  +  S)=if(CD0) 
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and  leads  to  a  value  of  5  =  AV2ln2  =  1.18A. 
In  the  case  of  a  Lorentzian 


n  S2  +  (a)-coQ) 

.where  8  is  the  HWHH,  the  moments  are  diverging  if  they  are  taken  over  all 
frequencies.  A  cutoff  must  be  introduced,  limiting  the  frequency  to  an  interval 
|o>-a>0|  -  a  which,  of  course,  must  be  much  larger  than  HWHH.  The  moments, 

neglecting  terms  — ,  are  then  found  to  be 
a 

M  TtCX 

The  large  ratio  of  — iT  =  —  indicates  that  this  functional  dependence  should 

only  be  tried  under  this  condition.  In  the  Gaussian  case  the  ratio  is  much  smaller 
with  a  value  of  3. 

3.5.    Heterogeneous  Spin  System  Relaxation 

If  spins  experience  different  environments  their  relaxation  times  are  also 
likely  to  be  different.  In  the  case  of  zeolite  there  are  at  least  two  distinguishable 
environments:  the  surface  and  the  center  of  the  cavity.  A  group  of  spins  shall 
be  called  a  spin  system.  The  following  discussion  is  restricted  to  two  different 
systems  a  and  b. 

For  the  moment  it  shall  be  assumed  that  the  spin  systems,  experiencing 
two  different  environments,  are  independent  of  each  other.  Each  spin  system 
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relaxes  therefore  independently  and  has  a  well-defined  relaxation  time  or, 
expressed  differently,  inherent  relaxation  rate  Ra  and  Rb.  There  is  also  a  distinct 
fraction  of  magnetization  pa  and  pb  associated  with  each  system. 

At  this  point  the  earlier  assumption  of  independence  between  the  spin 
systems  is  dropped.  The  consequence  is  that  the  spin  systems  interact  with 
each  other  undergoing  exchange  processes.  The  exchange  process  can  be  a 
transfer  of  magnetization  which  happens  by  spin  diffusion  via  spin  flip-flop,  or 
particle  diffusion.  This  interaction  between  spin  systems  is  a  more  realistic 
assumption  for  physical  systems.  The  outcome  for  this  model  is  that  the 
measured  relaxation  rates  A+  and  A"  and  measured  magnetization  fractions  C+ 
and  C~  do  not  correspond  to  the  physical  processes  one  is  interested  in.  In 
addition  to  the  previous  relaxation  rates,  there  are  relaxation  rates  ka  and  k„ 
between  the  spin  systems  which  are  not  measured  directly.  They  can  only  be 
extracted  as  a  complex  byproduct  to  the  overall  relaxation.  A  calculation8788  will 
reveal  the  involved  quantities. 

For  simplicity  of  the  following  equations,  a  reduced  magnetization  is 
introduced  for  the  longitudinal  relaxation  (saturation  recovery) 


and  the  transverse  relaxation 


M0 


The  equations  of  motion  for  the  magnetizations  are  Bloch  equations  extended  by 
the  additional  cross  relaxation  between  spin  reservoirs 
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d(SMa) 
dt 


=  -RaSMa-ka5Ma+kb5Mb 


=  -Rb5Mb  -  kb5Mb  +  ka5Ma 


where  the  above  definitions  are  employed.  These  coupled  differential  equations 
are  solved  for  the  magnetizations  in  terms  of  the  measured  or  apparent  fractions 
and  relaxation  rates 


5M,  =C-e-x''+C*e-i'' 


with  complicated  expressions  for 


A*  =  i(Ra  +  Rb  +  ka  +  kb) ± i([Ra  - Rb  +  ka  - kbf  +  4k akb) 


with 


Cf  =  SMI0(^^)±(5MI0-5MI0)-^ 


where  subscript  0  indicates  the  value  immediately  after  the  pulse.  Another 
important  piece  of  information  is  the  equation  of  detailed  balance 
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where  pa  and  pb  are  the  spin  fractions  of  each  group.  By  these  means  the 
apparent  relaxation  rates  become  dependent  on  the  magnetization  fractions  of 
each  group  by  substituting  for  one  of  the  ks. 

For  known  pa  and  pb,  which  can  be  assumed  for  the  case  of  zeolite  from 
the  fraction  of  molecules  in  contact  with  the  surface,  the  real,  inherent  relaxation 
rates  can  be  determined.  This  is  achieved  by  inverting  the  A+  and  X~  equations 
leading  to 

[(V-A-)C-+PbA--par-ka] 

"a  =  

Pb-Pa 

which  can  be  plotted  versus  ka.  Knowing  the  functional  form  of  Ra  and  Rb,  the 
functional  form  of  the  apparent  rate  A"  can  also  be  generated.  The  ka  for  which 
the  calculated  value  of  X  corresponds  to  the  experimental  value  of  A"  is  the 
exchange  rate  between  the  two  spin  systems.  The  values  of  Ra  and  Rb  can 
then  be  determined  from  the  known  ka  value.  The  remaining  difficulty  is  to 
reliably  deconvolute  two  exponential  decays.  This  is  even  a  problem  for 
sophisticated  computer  routines  because  of  possible  local  minima  in  parameter 
space.  As  a  consequence,  initial  guesses  which  are  too  far  away  from  the  true 
solution  lead  to  a  nonconvergent  fit. 

3.6.    Methods  of  Detection 

3.6.1.  Continuous  Wave  (cw)  Method 

In  continuous  wave  mode  the  system  is  in  a  steady  state  and  the  spin 
system  therefore  in  equilibrium  with  the  lattice. 
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3.6.1.1.  Q-meter-detection 

The  Q-meter  method  relies  on  a  change  in  impedance  as  the  resonance 
is  approached.  It  is  therefore  necessary  to  study  the  impedance  of  a  L-C  circuit. 
At  first,  the  induction  must  be  determined.  The  flux  d>  through  a  coil  which  is 
induced  by  a  circulating  current  /  is  defined  as 

<D  =  BA  =  (H  +  4^M)A  =  Re{L(1  +  4^)/0e"u' } 

where  77  is  the  filling  factor,  indicating  the  ratio  between  coil  volume  and  sample 
volume.  The  magnetization  is  related  to  the  field  via  the  complex  susceptibility 
M  =  xH.  This  functional  form  rescales  the  inductance  to  an  effective  value 
incorporating  the  effect  of  the  sample  magnetization  to 

LB„  =L(U4miz) 

The  voltage  generated  from  a  constant  current  source  is  proportional  to 
the  impedance  Zof  a  parallel  circuit  consisting  of  a  sample-filled  inductor  with  a 
series  resistance  r  and  a  capacitor  C 

1  |   J  ' 

r  +  icoL{\  +  4  xtjx)  {icoC)'^ 

It  is  useful  to  introduce  the  quality  factor  Q  as  the  ratio  of  energy  stored  and  the 
dissipated  power  in  the  coil 

Q)W  =Q){LI2  =  00L 

P       \ri2  r 
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If  the  circuit  is  tuned  to  the  resonance  condition,  co2LC  =  1,  the  impedance  Zcan 

(coL)2 

be  rewritten  in  terms  of  a  parallel  shunt  resistance  R  =  - — — 


Z  =  R 


and  the  quality  factor  is  expressed  as  Q  =  —  =  This  can  be  simplified  by 
assuming  Q  »  1  which  is  a  reasonable  assumption  for  real  circuits  with  Q  - 100 


Z  -     - \Anr\Qx\  =     -  /4 -4 m\Q%x '] 

The  relative  change  in  Zis  AZ  »  -AnRQT\x"  and  is  proportional  to  the  change  in 
the  detected  signal.  The  dispersion  is  not  measured  with  this  method  which  is  a 
drawback  in  itself.  It  does,  however,  have  the  advantage  of  avoiding  the  mixing 
of  dispersion  and  absorption  signals.  Another  problem  with  the  arrangement  is 
that  the  small  signal  voltage  is  superimposed  on  top  of  the  large  carrier  voltage. 

3.6.1.2.        Bridge  method 

The  last  mentioned  problem  is  eliminated  with  a  bridge  arrangement.  The 
carrier  voltage  V0  is  partially  compensated  by  a  bridge  or  "Tee"  so  that  the 
detection  amplifier  does  not  saturate 

V  =  V00-i4xQiiz)-Vi 


It  is  preferable  to  keep 
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VQ-V\  =  aV0ex$  >  8V 


in  order  to  generate  a  dispersion  or  absorption  signal  by  having  the  correct 
phase  on  the  compensating  voltage.  This  leads  to  an  equation  for  the  voltage  V 


This  way,  an  angle  0  =  0  leads  to  the  Q-meter-like  absorption  signal  but  an 


arrangement  in  practice  is  the  adjustment  of  the  phase  angle  and  the  constancy 
over  time. 

The  appropriate  apparatus  is  a  hybrid  Tee  which  has  an  inverting  and  a 
normal  input  power  branch.  One  branch  is  connected  to  a  50  Ohm  load,  the 
other  is  the  sample  path  and  is  matched  to  50  Ohm  with  an  impedance  matching 
device.  Scanning  through  the  resonance  with  a  sweeping  magnetic  field 
achieves  measurable  detuning  proportional  to  the  response  of  the  sample. 

Both  detection  schemes  are  supplemented  with  lock-in  amplification 
because  of  signal-to-noise  (S/N)  problems.  The  lock-in  acts  as  a  narrow  band 
filter  to  improve  the  S/N.  To  achieve  this,  it  modulates  the  swept  magnetic  field 
H0  with  an  audio  frequency  Q  and  an  amplitude  much  smaller  than  the  line  width 
to  avoid  distortion.  This  modulation  leads  to  a  differentiated  signal,  oscillating  at 
the  modulation  frequency.  This  signal  is  then  mixed,  i.e.  multiplied  by,  a  sinusoid 
at  the  same  frequency  and  averaged  over  time  T 


angle  <D  =  —  reverses  the  signal  into  pure  dispersion.  The  problem  with  this 


r 

S(T)  =  J [S0  sin(Qf)  +  n(t)]s\n(Qt)dt 

0 
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where  n(t)  is  the  noise  term,  taken  to  be  the  Nyquist  noise  at  temperature  0, 
which  is  only  meaningful  as  its  mean  square  value 


With  the  approximation  of  long  average  times  (QT»  1)  the  signal  to  noise  ratio 
becomes 


implying  that  the  band  width  of  the  lock-in  is  Av  =  2/T  which  improves  the  S/N 
by  longer  integration  times  at  the  expense  of  reduced  resolution. 

3.6.2.  Transient  Method 
3.6.2.1.        Coherent  pulses 

The  coherent  pulse  method  opens  up  the  opportunity  of  exploring  line 
shapes  much  narrower  than  those  limited  by  field  inhomogeneities.  Many  pulse 
sequences  are  in  use.  Motionally  narrowed  line  shapes  are  of  special  interest  in 
this  context. 

The  Hamiltonian  for  the  system  in  the  rotating  frame,  during  the 
application  of  a  pulse  along  the  x-direction,  is  given  by 


N    ^T4k@R      °14  k@R 
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where  the  H,  field  is  much  larger  than  the  field  in  z-direction.  This  justifies  the 
approximation  H  =  -yHJx  for  the  short  time  of  the  pulse.  Applying  for  a  time  t 
means  rotating  the  magnetic  moment  by  an  angle  0  =  yH,t  described  by  the 
rotation  operator 

R(Q)  =  e-iH'  =  ew-'  =  e"'- 

Due  to  the  strength  of  the  applied  rf  field,  ideally  turns  the  spins  in  a  negligible 
time.  The  pulse  sequence  described  in  the  following  is  the  original 
90, -7-180,-r  sequence  proposed  by  Hahn89  and  is  graphically  shown  in 
figure  3. 

The  steady  field  H0  leads  to  an  excess  orientation  of  spins  in  z-direction 
(Figure  3a).  After  the  steady  state  has  been  reached,  the  Hrfield  is  switched  on 
for  a  short  time  in  order  to  tilt  the  magnetization  by  90  degrees  (Figure  3b)  from 
Mz  into  -My: 


( 

M. 


The  same  coil  producing  H,  picks  up  the  voltage  induced  by  the  magnetization 
rotating  in  the  x-y  plane.  Without  relaxation  effects,  as  expressed  by  7",  and  T2, 
the  alternating  induced  voltage  would  persist  forever.  In  reality  T2  plays  a  role 
but  the  7;  effect  is  negligible,  considering  7,  »  T2  for  a  solid.  This  gives  rise  to 
dephasing  which  destroys  the  ideally  persisting  voltage  by  phase  cancellation 
(Figure  3c).  A  free  induction  decay  (FID)  results  from  the  superposition  of  all 
individual  spin  signals  at  slightly  different  fields.  The  signal  decreases 
exponentially  in  the  case  of  a  Lorentzian  line  with  a  time  constant  T2.  This  T2  is 
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Figure  3.  Spin  echo  formation1 
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a  sum  of  all  effects,  especially  the  inhomogeneity  of  H0.  It  causes  some  spins  to 
advance,  others  to  lag  behind.  But  this  is  not  what  really  determines  the  physics. 
The  goal  is  to  eliminate  this  unphysical  influence  by  applying  a  180  degree  pulse 

The  180  degree  pulse  produces  a  mirror  image:  the  slowest  spins  are  ahead  of 
the  faster  spins  (Figure  3d).  The  /r-pulse  acts  like  a  time  reversal  operator. 
After  the  same  time  t  between  the  nl2  and  the  ^-pulse  the  fastest  spins  have 
caught  up  with  the  slowest.  Since  the  inhomogeneity  of  the  steady  magnetic 
field  has  not  changed  during  the  time  2t  the  almost  same  FID  occurs  after  this 
time  (Figure  3e).  The  difference  is  that  the  voltage  has  the  opposite  sign  of  the 
rotated  magnetization  and  a  decreased  amplitude.  For  the  90,-T-180y-r 
sequence  the  echo,  called  "solid"  echo,  has  the  same  sign  as  the  original  FID. 
The  reduction  in  signal  amplitude  is  caused  by  effects  that  give  insight  into  the 
microscopic  behavior.  Due  to  collisions  during  the  two  time  intervals  x  only  a 
diminished  spin  echo  occurs,  e.g.  as  a  result  of  diffusion  effects. 
Ty  effects  reduce  the  measured      according  to 

1      1  1 
7"an  "  Tt  +  27, 

For  T,  »  T2,  as  is  the  case  of  solids,  this  is  not  a  crucial  difference. 

It  is  important  to  notice  that  this  pulse  sequence  is  only  able  to  refocus 
spin  and  not  dipolar  components.  This  means  that  the  echo  amplitude, 
corrected  by  the  spin-spin  relaxation  time,  is  a  measure  of  the  liquid-like  fraction 
of  the  sample. 


CHAPTER  4 
EXPERIMENTAL  SETUP 


4.1.     NMR  Setup 

Many  components  of  the  NMR  spectrometer  are  used  for  pulse  as  well  as 
continuous  wave  mode.  With  the  necessary  expertise  it  is  therefore  possible  to 
switch  between  the  two  modes  within  minutes.    This  makes  the  existing 
apparatus  versatile  for  the  different  information  attainable  by  NMR. 
4.1.1.  Pulse  Apparatus 

A  diagram  with  an  overview  of  the  complete  pulse  NMR  apparatus  is 
shown  in  figure  4.  The  system  will  be  discussed  in  detail,  closely  following  the 
diagram. 

The  radio  frequency  (rf)  path  originates  at  the  rf-generator  which  supplies 
the  ultrahigh  frequency  (UHF)  of  268  MHz  continuously.  This  frequency 
corresponds  to  a  magnetic  field  H0  of  6.3  T  for  hydrogen  according  to  the 
Larmor  relationship.  This  magnetic  field  is  in  the  upper  attainable  region  of  the 
available  magnet  and  is  chosen  for  high  resolution  and  sensitivity,  i.e.  good 
signal  to  noise  ratio  (S/N),  obtainable  at  high  frequencies.  The  S/N  is  a  function 
of  the  magnetic  field  to  the  3/2  power.  It  therefore  pays  to  work  in  a  high 
frequency  regime  as  long  as  no  other  effects,  such  as  the  skin  depth  are 
counteracting.  Another  argument  is  the  short  dead-time  which  is  proportional  to 
the  reciprocal  of  the  frequency.  A  continuous  rf  signal  is  not  desired  for  the 
pulse  mode.  For  this  reason  the  pulse  generator  produces  rectangular  signals  of 
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Figure  4.  Pulse  NMR  diagram 
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variable  length  T  and  delay  time  r  between  them.  Sequences  can  be 
programmed  from  the  laboratory  computer  or  selected  manually.  Four  outlets 
connect  four  different  gates  which  allow  for  a  phase  shift  of  0,  90,  180  and  270 
degrees.  At  the  gates  the  rf  frequency  is  blocked  and  passes  only  for  the  time 
adjusted  by  the  pulse  length  T.  To  reduce  noise  from  the  rf  generator  pulse  or 
leakage  through  the  gates  in  the  off-state,  crossed  diodes  are  placed  after  the 
gates.  The  diodes  do  not  conduct  for  voltages  less  than  0.5  V  which  is  above 
the  noise  level.  For  higher  voltages  they  have  no  impact.  A  source  of 
frequencies  other  than  the  generated  268  MHz  is  the  rectangular  shape  of  the 
pulse.  A  Fourier  analysis  shows  that  higher  odd  multiples  of  this  frequency  are 
present  and  useful  power  is  lost.  This  contribution  is  also  diminished  by  the 
crossed  diodes. 

Pulse  NMR  in  solid-state  physics  requires  a  large  amplitude  of  rf  which  is 
associated  with  a  short  pulse.  The  values  in  the  present  apparatus  are  15  G  for 

the  H,  field  which  is  equivalent  to  a  -^-pulse  duration  of  less  than  20  ^isec.  In  a 

negligible  time  the  spins  have  to  be  tipped  over  by  the  desired  angle  to  achieve  a 
quasi-instantaneous  situation  where  time  effects  are  negligible.  Another  reason 
for  the  short  pulse  is  the  need  for  wide  coverage  of  the  whole  frequency 
spectrum.  These  conditions  are  met  by  the  two  high  power  amplifiers  of  40  and 
300  Watts. 

A  quarter-wave-length  line  with  attached  crossed  diodes  is  mounted  in 
parallel.  This  device  serves  as  noise  reduction  when  the  pulses  are  not  turned 
on.  For  an  understanding  of  this  effect,  transmission  line  theory90  is  necessary. 
There  are  special  characteristics  for  A/2  and  A/4  lines.  In  the  case  of  a  A/2  line 
the  impedance  on  one  end  is  simply  transformed  to  the  other  end,  except  for  the 
additional  resistance  of  the  used  cable.   For  a  A/4  line  the  input  and  output 
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impedances  Z,,  Z0  and  the  characteristic  cable  impedance  -  Z  =  50  CI  in  the 
standard  coaxial  case  -  are  related  to  each  other  by  Z,Z0  =Z2.  These  special 
cases  are  a  consequence  of  the  general  formula  for  transmission  lines 

Z,  _  Z0  +  \Z\ankl 
Z  ~  Z  +  \Z0\ankl 

2n 

where    =  —  is  the  wave  number.  For  high  voltages  (above  0.5  V)  the  diodes 

At 

conduct  and  represent  a  low  impedance.  For  the  low  impedance  at  the  end  of 
the  crossed  diodes  the  resistance  at  the  other  end  of  the  A/4  line  is  very  high. 
This  implies  that  the  high  voltage  level  of  the  frequency  corresponding  to  the 
wavelength  of  the  spectrometer  is  not  affected.  Other  frequencies  do  not  match 
and  are  dampened  proportional  to  their  frequency  deviation.  For  voltage  levels 
below  0.5  V  the  diodes  do  not  conduct  and  the  net  effect  is  an  almost  zero 
resistance  at  the  carrier  frequency.  This  mechanism  decreases  this  frequency 
noise  component  but  is  less  efficient  at  other  frequencies. 

The  center-piece  of  the  pulse  spectrometer  is  the  duplexer.  Its  task  is  to 
connect  the  transmitter,  the  sample  and  the  receiver,  but  to  simultaneously 
establish  a  separation  of  the  transmitter  from  the  receiver.  This  requirement 
becomes  clear  when  the  power  of  the  transmitter  (300  W)  is  compared  to  the  0.1 
mW  power  of  the  sample  signal:  the  highly  sensitive  receiver  would  be 
destroyed  by  the  transmitter  power.  The  duplexer  consists  of  two  impedance 
matching  circuits.  One  transforms  the  cable  impedance  of  50  Q  to  the  3  Q  of 
the  series  NMR  resonance  circuit.  The  other  circuit  matches  the  3  Q  to  the  400 
Q.  of  the  UHF  amplifier  with  an  amplification  factor  of  approximately  1000. 

The  low-level  amplifier  is  a  narrow-band,  3  stage  solid  state  amplifier  and 
is  another  key  component  of  the  spectrometer.  A  fourth  stage  is  designed  to 
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provide  variable  matching  to  the  input  impedance  of  the  incoming  signal.  The 
amplifier  gain  must  be  high  but  also  robust  against  self-oscillation.  The  narrow 
bandpass  is  achieved  with  a  resonance  circuit  for  each  stage  consisting  of  a 
fixed  coil  and  a  tunable  capacitor.  To  reduce  readiness  for  self-oscillation 
resistors  were  soldered  across  the  coils,  at  the  expense  of  increased  band  width 
and  reduced  gain.  The  used  tunable  capacitors  were  of  high  mechanical  and 
electrical  quality  in  order  to  reduce  mechanical  vibrations  and  allow  for  precise 
fine  tuning.  The  DC  power  supply  to  the  double-gate  transistors  (3N21 1,  3N213) 
was  decoupled  from  the  rf  path  by  large  inductors.  The  individual  compartments 
were  also  carefully  separated  from  each  other  to  prevent  feedback  between 
amplification  stages.  The  double-gate  transistors  were  carefully  biased  in  order 
to  operate  in  a  region  of  high  forward  admittance.  This  is  the  case  for  about  2  V 
between  the  source  and  gate  2  and  zero  voltage  between  the  source  and  gate  1 , 
which  is  the  signal  carrying  input.  Exact  tuning  of  all  resonance  circuits  is  the 
key  to  high  and  stable  amplification  levels.  The  entire  amplifier  was  placed  into  a 
brass  box  which  was  temperature  regulated.  This  seemed  necessary  due  to  the 
temperature  dependence  of  the  transistors  and  the  variable  temperatures  in  the 
laboratory.  A  temperature  regulator  with  a  temperature  probe  was  purchased.  A 
heater  was  dimensioned  to  deliver  enough  power  to  keep  the  temperature  at  32 
centigrades.  This  value  was  chosen  because  it  was  above  any  reasonable  room 
temperature. 

The  receiver  is  protected  by  two  sets  of  crossed  diodes  inside  the 
duplexer  which  shunt  the  high  level  transmitter  power  and  limit  the  input  voltage 
into  the  amplifier.  This  could  also  have  been  achieved  with  an  additional  A/4 
line  in  the  main  line.91  Its  function  is  described  as  follows.  For  the  high 
transmission  levels  the  diodes  would  turn  on  and  the  upper  point  of  the  A/4  line 
would  represent  a  very  high  resistance  from  the  transformation  of  the  short- 
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circuit  at  its  end.  This  would  decouple  the  receiver  stage.  The  small  sample 
signal  would  only  see  a  negligible  resistance  and  could  be  processed  further  on 
its  way  to  detection.  Matching  problems  prevented  this  method  from  being 
employed. 

The  NMR  circuit  is  tuned  to  resonance  (co2LC  =  1)  at  268  MHz,  where  the 
series  circuit  shows  a  minimum:  |Z|2  =  r2  +(1-  ofLC I  coC).  This  is  a  desirable 
property  of  the  series  circuit  because  the  small  resistance  of  r  =  3  Q  only 
requires  a  few  volts  of  applied  power.  An  alternative  parallel  circuit  would 
require  several  kV  for  the  same  power,  which  would  lead  to  breakdown  along  the 
1.5  m  long  coaxial  cable  to  the  coil.  The  sample  coil  must  be  dimensioned 
appropriately  in  conjunction  with  the  tunable  capacitor  for  the  applied  frequency. 
The  inductance  for  small  coils  is  determined  by  the  formula92 

L_  N2a2 
9a  +  10b 

where  2a  is  the  diameter  of  the  winding,  /Vthe  winding  number  and  b  the  length 
of  the  coil.  As  long  as  b  >  0.8a  this  is  a  good  estimate.  The  inductance  is 
chosen  relatively  large  to  obtain  a  large  quality  factor.  This  is  restricted  by  the 
corresponding  value  of  the  capacitor  which  should  not  be  below  a  few  pF.  The 
present  values  are  2.6  pF  and  150  nH.  The  quality  factor  is  about  100.  The 
used  coaxial  cable  is  important  because  it  should  not  introduce  excessive 
losses.  A  homemade  cable  was  a  combination  of  copper-nickel  tubing  and 
flexible  braid  in  the  lower  part  of  the  cryostat.  Along  the  cold  finger  a 
commercially  available,  high  quality,  rigid  coaxial  cable  was  used  in  order  to 
minimize  spurious  strain  effects.  The  pulse  system  is  adjusted  to  the  coil.  The 
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inductance  is  matched  with  the  input  capacitance  of  the  duplexer  which  can  also 
be  influenced  by  the  coax  cable  length. 

The  probe  signal  is  processed  in  a  phase  detector  after  amplification  to 
eliminate  the  carrier  frequency.  Preamplification  before  transformation  to  low 
frequencies  yields  an  S/N  advantage,  because  1/f  noise  effects  are  much 
smaller  at  high  frequencies.  The  phase  detector  or  mixer  receives  its  reference 
signal  from  the  rf-generator  and  passes  only  signals  which  are  modulated  with 
this  frequency  and  higher  harmonics  which  are,  however,  not  relevant.  By  these 
means  the  S/N  ratio  is  increased.  A  phase  shifter  allows  for  changing  the 
reference  phase.  The  possibility  of  detecting  absorption  and  dispersion  signals 
arises. 

Another  amplification  by  250  and  a  signal  averager  follow.  Signal- 
averaging  improves  the  signal  to  noise  ratio  by  repeating  the  experimental 
procedure  and  adding  up  the  digitized  amplitudes  of  the  signal.  The  signal 
improvement  increases  as  the  square  root  of  the  number  of  repetitions. 

The  broader  the  NMR  line,  the  broader  the  bandwidth  of  the  amplifiers 
must  be  in  order  to  avoid  omitting  parts  of  the  frequency  spectrum.  A  larger 
bandwidth  deteriorates  the  noise  and  gain  characteristics. 

The  elaborate  system  at  high  frequency  with  the  homogeneous  magnet 
and  the  noise  reducing  schemes  allows  for  detection  of  1016  spins  at  1  K. 

Before  cooldown  the  NMR  system  needs  to  be  checked  for  pulse 
propagation  to  the  sample  coil.  An  auxiliary  coil  is  placed  around  the  NMR  coil 
and  the  signal,  as  a  response  of  the  applied  pulse,  is  induced  in  the  pick-up  coil 
and  detected.  The  connected  oscilloscope  displays  a  voltage  of  40  V  peak  to 
peak  via  the  weak  coupling.  This  value  represented  a  standard  for  good 
operating  conditions. 
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4.1.2.  Continuous  Wave  (cw)  Apparatus 

Many  components  of  the  pulse  apparatus  are  also  used  for  cw  NMR.  For 
this  reason,  the  modifications  and  differences  compared  to  the  pulse  system  are 
described  first.  The  setup  is  then  explained  following  the  complete  cw  block 
diagram  displayed  in  figure  5. 

The  entire  pulse-generating  equipment  including  the  high  power 
amplifiers,  gates  and  pulse  generator  is  bypassed.  Instead  of  employing  a 
duplexer  there  are  two  different  methods  of  how  continuous  wave  data  are 
obtained.  Both  principles,  the  Q-meter  and  the  hybrid  T  have  already  been 
discussed  in  principle.  The  hybrid  T  supplies  the  receiver  only  with  the  sample 
response  while  the  Q-meter  is  also  subject  to  the  input  power. 

The  hybrid  T  is  an  analogue93  of  the  magic  T  at  high  frequencies.  It  works 
like  a  bridge  and  operates  with  a  dummy  load,  which  does,  however,  waste  half 
of  the  transmitted  power.  The  rf  generator  input  is  split  in  half  and  reaches  the 
dummy  load  and  the  NMR  circuit,  which  are  both  matched  to  50  Q.  As  long  as 
the  sample  is  off-resonance  the  hybrid  T  is  matched  so  that  the  two  paths  cancel 
each  other,  because  the  sample  branch  is  out  of  phase  by  180  degrees.  The 
matching  is  spoiled  at  resonance  and  only  the  NMR  change  is  detectable  by  the 
receiver.  This  is  the  ideal  case;  in  reality  a  small  deliberate  deviation  from  ideal 
phase  or  amplitude  matching  is  adjusted  as  a  signal  carrier.  The  small  error 
picks  out  the  component  of  the  NMR  signal  in  the  direction  displaying  a 
dispersion  or  absorption  signal,  respectively.94-95  The  advantage  is  an 
amplification  effect  and  a  deliberate  choice  of  signal  type  instead  of  a  mixture. 

Another  major  difference  to  the  pulse  apparatus  is  the  modulated 
sweeping  field.  The  steady  field  H0  can  be  adjusted  by  a  sweep  coil  in  order  to 
ramp  the  magnetic  field  through  resonance.    This  is  more  practical  than 
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Figure  5.  Continuous  wave  (cw)  diagram 
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sweeping  the  frequency  through  resonance  because  all  spectrometer 
components  are  carefully  frequency-tuned.  The  change  in  the  static  magnetic 
field  in  z-direction  is  not  the  only  modification.  In  addition,  the  auxiliary  sweeping 
field  is  sinusoidally  modulated  by  79.3  Hz.  This  frequency  is  deliberately  chosen 
as  an  odd  value,  such  that  no  other  frequency  matches  it  accidentally.  The 
modulation  changes  the  nature  of  the  sample  response.  Instead  of  the 
absorption  signal  itself,  its  first  derivative  is  obtained  because  the  sine  wave  is 
projected  onto  the  absorption  lineshape.  Thus,  the  absorption  lineshape  is 
probed  yielding  the  strongest  signal  for  the  steepest  slopes  and  therefore 
monitoring  the  first  derivative  of  the  line.  Another  change  is  the  desired  shift  of 
the  measurement  to  a  higher  frequency  which  reduces  1/f  noise  and  allows  for 
using  a  lock-in  amplifier.  The  lock-in  has  the  effect  of  narrowing  the  bandwidth 
and  improving  the  signal  to  noise  ratio  sharply.  Experimentally,  the  lock-in 
parameters  have  to  be  chosen  carefully.  The  time  constant  is  a  way  to 
smoothing  the  recorded  signal  by  integration.  It  is  important,  however,  not  to 
disguise  narrow  features  with  too  generous  time  constants.  The  lock-in  is 
phase-sensitive  and  a  phase  offset  by  90  degrees  totally  eliminates  any  existing 
signal.  It  is  clear  that  the  right  phase  choice  is  important  for  maximizing  the 
signal. 

The  complete  cw  system  can  now  be  understood  easily.  Again  the  rf  path 
originates  at  the  frequency  generator.  The  power  is  split  into  the  signal  and 
reference  branch.  The  signal  branch  feeds  -  via  an  adjustable  attenuator  to 
control  the  power  input  -  the  Q-meter  or  the  hybrid  T  which  leads  to  the  detection 
of  the  sample  response.  The  weak  signal  is  amplified  in  the  low-level  amplifier 
already  described  in  the  pulse  mode  section.  The  specialty  of  the  input  stage  of 
this  amplifier  in  Q-meter  operation  is  that  it  performs  two  tasks  at  once:  it 
matches  the  impedance  and  acts  as  the  resonance  capacitor  of  the  tank  circuit. 
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This  double  task,  which  is  performed  by  the  appropriate  tapping  off  the 
resonance  coil  of  the  first  stage,  is  critical  to  the  observed  signal.  In  hybrid  T 
operation  the  impedance  matching  is  already  provided  by  two  capacitors,  one  to 
ground  and  the  other  to  feed  rf.  This  exterior  impedance  matching  has  also 
been  successfully  applied  for  the  Q-meter  alternative.  The  signal  is  mixed  with 
the  phase  adjustable  reference  signal  at  the  phase  detector  or  mixer  and  then 
fed  into  another  amplifier.  The  final  stages  are  the  lock-in  and  the  chart 
recorder. 

4.2.    Dilution  Refrigerator 

In  order  to  cool  down  to  temperatures  below  4.2  K  a  dilution  refrigerator  is 
used.  Its  components  and  the  cooling  process  are  described. 

A  dilution  refrigerator  is  operated  with  a  3He/4He  mixture  circulating  in  a 
closed  gas  handling  system.  The  4He  essentially  represents  an  inert 
background  while  3He  is  cycling  through  the  refrigerator.  Two  pumps  drive  the 
gas  flow  and  determine  the  cooling  power  with  their  pumping  speed.  The 
gaseous  3He  at  room  temperature  returns  from  the  pumping  system,  enters  the 
cryostat  and  condenses  at  a  small  bath  of  liquid  helium  held  at  1  K,  the  "1-K  pot". 
Thermal  contact  is  established  by  passing  the  gas  through  a  cupronickel  tube 
which  is  wound  around  and  then  immersed  in  the  bath. 

The  pressure  drop  of  liquid  3He  from  the  condensation  condition  (100 
Torr)  to  values  corresponding  to  the  still  temperature  (3  Torr)  makes  it  necessary 
to  introduce  a  flow  limiting  impedance.  The  pressure  drop  for  a  circulation  rate 
N  is  determined  by  Poiseuille's  equation 


Ap  =  77Z/vV3 
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where  ri  is  the  viscosity,  Z  the  flow  limiting  impedance,  depending  on  tube 
parameters  /  and  d,  and  V3  the  molar  volume  of  3He.  Without  the  impedance, 
the  3He  is  in  danger  of  re-evaporating  which  gives  rise  to  a  subsequent  heat 
load  by  recondensation  and  turbulences.  On  its  way  to  the  mixing  chamber  the 
3He  must  be  precooled  because  the  dilution  cooling  mechanism  only  starts  at 
temperatures  below  the  tricritical  point  of  0.86  K.  Below  this  temperature  the 
3He/4He  mixture  spontaneously  separates  into  a  3He  rich  phase  floating  on  top 
of  a  4He  rich  phase  due  to  the  mass  difference.  The  3He  tube  is  therefore  led 
through  the  still  and  various  heat  exchangers  for  precooling  purposes. 
Especially  the  heat  exchangers  represent  an  essential  means  of  precooling  and 
increase  the  efficiency  of  the  refrigerator.  Two  different  kinds  are  in  use: 
continuous  and  discrete  heat  exchangers.96  The  continuous  kind  consists  of  an 
outer  tube  which  guides  the  inner  warmer  tube.  Instead  of  truly  discrete  heat 
exchangers  with  single  copper  blocks,  "quasi-continuous"  heat  exchangers  are 
employed.  The  heat  transfer  is  ensured  by  a  large  surface  area  of  sintered 
silver.  Despite  the  large  contact  area  the  flow  impedance  is  not  excessive.  The 
simple  continuous  counterflow  heat  exchangers  are  limited  by  the  Kapitza 
resistance  and  viscous  heating.  The  Kapitza  resistance  is  due  to  the  acoustic 
mismatch  between  the  liquids  and  the  separating  metal  tube 


where  A  is  the  effective  surface  area,  AT  the  temperature  gradient  across  the 
interface  and  Q  the  thermal  flux.  As  a  consequence  the  simple  heat  exchangers 
fail  at  temperatures  below  100  mK  and  the  more  sophisticated  designs  are 
needed  for  lower  temperatures.  More  efficient  heat  exchangers  allow  for  higher 
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circulation  rates  without  increasing  the  temperature.97  When  the  3He  leaves  the 
last  exchanger  it  is  colder  than  the  tricritical  temperature,  once  the  refrigerator  is 
working,  and  fills  up  the  concentrated  3He  phase.  At  this  point  the  dilution 
cooling  process  sets  in  as  3He  quasi-evaporates  from  the  concentrated  phase 
into  the  "He  phase,  or  dilute  3He  phase.  The  cooling  power  can  be  quantified 

Qc=N*[H3D(Tu)-H3C(TM)] 

where  TM  is  the  temperature  of  the  mixing  chamber  and  Hthe  enthalpy  of  dilute 
and  concentrated  3He.  The  available  cooling  power  QM  at  the  mixer  is 
equivalent  to  the  external  heat  leak  in  equilibrium  and  is  reduced  by  the  heat 
load  QL  of  the  incoming,  warmer  3He  which  has  the  heat  loaded  temperature  TL 

qm=qc-Ql  =  n3[h2D(tm)-h,c(tm)-h,c(tl)+h3C(tm)] 

=  N,[H3D(TU)-H3C(TL)] 

The  calculation  is  simplified  by  the  assumptions  that  only  3He  is  circulated  and 
no  frictional  heating  arises.  It  turns  out  to  be  important  for  the  performance  of 
the  refrigerator  to  establish  a  high  circulation  rate.  For  further  results  the 
involved  enthalpies  are  needed.  The  dilute  phase  can  be  interpreted  as  an  ideal 
Fermi  gas  with  the  specific  heat 


where  TF  =  0.38  K  at  xD  =  0.064,  the  maximum  3He  concentration  in  the  4He 
background  due  to  the  Pauli  principle.  For  the  concentrated  phase 
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C„  =  247 


moIeK 


which  is  found  experimentally.  The  chemical  potential  is  constant  between  the 
concentrated  and  dilute  phase  (/i3C  =  ^3D)  at  the  coexistence  curve 


The  entropy  is  calculated  to  be 


and  the  enthalpy 


H3C  =  H3  (0)  +  TjC3CdT  =  H3  (0)  + 1 272 


Using  this  information  leads  to  a  numerical  result  for  the  cooling  power 
6M=/N/3(95rM2-12rL2) 


Two  limits  can  be  investigated.  In  the  case  of  a  single  shot  experiment  or  the 
use  of  ideal  heat  exchangers,  i.e.  when  TL=TU,  the  cooling  power  is  maximal 
and  determined  by 
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In  the  other  limit  the  cooling  power  approaches  zero  if  -r-  =  2.83.  This  again 

stresses  the  need  for  efficient  heat  exchangers. 

The  mixing  chamber,  where  the  cooling  power  is  generated,  is  thermally 
linked  to  the  sample  cell  (Figure  6)  by  a  cold  finger  consisting  of  annealed 
copper  bars.  The  heat  contact  between  the  mixing  chamber  and  the  cold  finger 
is  established  by  a  cake  of  sintered  silver.  The  cold  sintering  is  performed  at  200 
atm  using  a  hydraulic  press  and  a  specifically  machined  piston  to  exert  even 
force  on  the  silver  powder.  Before  sintering,  the  copper  surface  must  be  cleaned 
with  striking  solution,  applying  voltage  in  reversed  polarity,  hence  reducing  any 
disturbing  cations.  The  so  prepared  surface  is  then  galvanized  with  a  solution  of 
silver  salt  as  a  base  for  the  sintering. 

After  crossing  the  phase  boundary  inside  the  mixer,  the  3He  atoms  return 
to  the  still  through  the  heat  exchangers,  driven  by  an  osmotic  pressure  gradient. 
The  temperature  in  the  still  must  be  high  enough  to  establish  a  sufficient  3He 
vapor  pressure  in  order  to  maintain  a  certain  circulation  rate  through  the  system. 
External  power  is  supplied  to  the  still 

Qs=N,[H3(vaP)-H3(cond)} 

This  implies  that  the  dilution  refrigerator  provides  more  cooling  power  when  the 
still  is  moderately  heated!  Almost  exclusively  3He  gas  leaves  the  still  because 
the  vapor  pressure  of  4He  is  only  about  5%  of  the  value  for  3He  at  these 
temperatures.  Impurities  such  as  air,  hydrogen  and  oil  derivates  are  separated 
from  the  3He  by  passing  the  gas  through  a  2-stage  liquid  nitrogen  trap  and  then 
pumping  3He  back  into  the  system.  By  circulating  3He,  the  process  is 
maintained  continuously.  Without  large  heat  leaks,  temperatures  of  50  mK  can 


89 


From 
Mixing 
Chamber 


Plastic 


Figure  6.  NMR  cell 


90 


be  sustained  over  a  long  period  of  time  by  refilling  the  4  K  He  bath.  The  1  K  bath 
is  filled  continuously  by  an  inlet  through  a  suitable  impedance  from  the  4  K  bath. 

In  a  dilution  process  the  number  of  3He  is  independent  of  temperature  but 
the  cooling  power  decreases  with  T£.  The  lower  the  mixer  temperature,  the  less 
cooling  power  is  available.  This  is  consistent  with  a  version  of  the  3rd  law  of 
thermodynamics,  stating  that  absolute  zero  temperature  cannot  be  attained. 

In  an  evaporation  process  the  number  of  molecules  removed  from  the 
vapor  is  proportional  to  the  vapor  pressure  which  decreases  exponentially.  The 
cooling  capacity  per  molecule,  however,  remains  approximately  constant.  For 
lower  temperatures  the  dilution  process  is  more  efficient  because 
T2  >Cexp(-1/7")  which  makes  the  dilution  refrigerator  an  efficient  way  of 
cooling  to  mK. 

The  presence  of  impurities  poses  a  problem  to  the  small  diameters  of  the 
impedances  in  the  closed  cycle  of  the  refrigerator .  If  the  system  has  a  leak  to 
air  the  gas  path  for  helium  will  start  plugging  up  at  nitrogen  temperatures.  In  the 
present  system  several  leaks  occurred  during  the  course  of  operation.  In  the 
external  gas  handling  system  three  different  leaks  were  discovered  and  fixed. 
Another  leak  at  a  quasi-continuos  heat  exchanger  inside  the  vacuum  can 
opened  up.  Such  a  leak  can  cause  losing  valuable  3He  and  prevents  from 
obtaining  low  temperatures  because  the  3He  acts  as  exchange  gas  to  the 
surrounding  4.2  K  helium  bath.  Also  the  shaft  seals  of  the  mechanical  booster 
and  the  sealed  oil  pump  had  to  be  exchanged.  Another  problem  are  pump  oil 
derivates  which  can  also  block  the  gas  pipes.  With  the  problem  of  contaminated 
mixture  gas  was  dealt  by  cycling  the  gas  through  a  trap  immersed  in  liquid 
helium.  This  is  the  most  efficient  way  of  cleaning  helium  gas  from  impurities. 
Another  way  to  prevent  clogging  in  the  pipes  was  to  run  the  refrigerator  with  a 
small  amount  of  helium  gas  upon  initial  cooldown  from  room  temperature.  As  a 


91 


consequence,  the  mixture  and  any  impurities  are  in  motion  and  cycled  through 
the  nitrogen  trap  so  that  the  probability  of  blocking  is  minimized. 

4.3.  Substrates 

4.3.1.  Zeolite 

The  name  "zeolite"  is  Greek  and  means  "boiling  stone".  It  is  a 
consequence  of  its  property  to  boil  and  melt  to  a  glass  as  it  is  heated  in  a  flame. 
Before  reaching  high  temperatures,  above  500  centigrades,  the  zeolite  releases 
water  but  does  not  disintegrate.  This  is  in  contrast  to  most  other  water-bearing 
crystals  and  results  from  the  porous  structure  of  zeolite.  There  are  other  names 
associated  with  zeolites  alluding  to  their  unique  properties:  Molecular  Sieve, 
Solid  Solvent  or  Ion  Exchanger. 

There  exists  a  great  variety98  of  naturally  occurring  as  well  as  synthetic 
zeolites.  All  of  them  are  defined  "as  aluminosilicates  with  a  frame  work  structure 
enclosing  cavities  occupied  by  large  ions  and  water  molecules,  both  of  which 
have  considerable  freedom  of  movement,  permitting  ion  exchange  and 
reversible  dehydration."99  They  are  grouped  into  Sodalites  (Faujasites), 
Chabazites,  Philipsites,  Analcimes  and  Mordenites.100  The  most  commonly  used 
zeolites  are  the  synthetically  prepared  and  commercially  available  (Union 
Carbide  Corporation,  Box  372,  South  Plainfield,  NJ  070880)  zeolites  X  and  A  of 
the  sodalite-faujasite  group.  Only  their  structure  will  be  discussed. 

The  structure101-102  of  zeolite  is  complicated  but  completely  regular  and 
resembles  a  "jungle-gym."  Figure  7  provides  a  simplified  view  of  the  structure  of 
zeolite  A.103  The  fundamental  building  block  of  any  zeolite  is  a  tetrahedron  of 
four  oxygen  ions  with  a  centered  aluminium  or  silicon  ion.  Each  oxygen  carries 
two  negative  charges,  one  of  them  being  compensated  by  Si.  The  remaining 
single  charge  on  the  oxygen  enables  it  to  combine  with  other  Si  ions  and  to 
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Figure  7.  Structure103  of  zeolite  A 
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extend  the  crystal.  If  Si  is  substituted  by  Al  one  less  charge  is  neutralized  and 
an  additional  cation,  such  as  sodium  or  calcium,  is  needed.  The  exchangeable 
cations  attach  themselves  loosely  to  the  oxygens  and  act  as  barriers  into  the 
channels  connecting  the  pores.  These  additional,  loosely  bound  cations  are  the 
reason  for  the  facilitated  exchange  of  ions,  such  as  Na  by  Ca,  in  zeolites. 

The  structure  formed  by  the  tetrahedra  is  similar  to  the  structural  unit  in 
sodalite,  the  "sodalite  unit."  It  contains  24  (Si,  Al)  and  36  oxygens  and  forms  a 
truncated  octahedron  with  one  tetrahedron  at  each  corner.  To  understand  this 
one  needs  to  know  that  a  truncated  octahedron  consists  of  eight  hexagonal 
faces,  six  square  faces,  24  vertices  and  36  edges. 

At  this  point  the  difference  between  zeolite  type  A  and  X  comes  into  play. 
For  zeolite  A  the  octahedrons  assemble  along  the  square  faces  (Figure  7)  while 
for  zeolite  X  the  hexagonal  faces  connect.  This  produces  a  central,  truncated 
cube  octahedron  with  an  internal  cavity  diameter  of  9  A  for  type  5A,  consisting  of 
eight  sodalite  units  on  a  simple  cubic  lattice.  For  type  X  a  tetrahedral 
arrangement  of  10  sodalite  units,  as  in  diamond,  is  formed  with  an  internal  cavity 
diameter  of  13  A  and  a  quasi  spherical  shape  with  four  caps  cut  off.  These 
cages  are  termed  a  cages  or  super-cages,  in  contrast  to  the  /?  cages  within  the 
initial  octahedron  of  6.6  A  diameter.  The  p  cages  are  too  small  to  be  occupied, 
except  for  crystal  water,  because  their  connecting  channels  are  only  2.2  A.  The 

s 

connecting  channels  between  the  a  cages  in  type  5A  are  about  5  A  wide, 
formed  by  eight  bridge  oxygens.  In  type  X  the  connecting  channels  are  8  A 
wide,  formed  by  a  ring  of  twelve  oxygen  ions.  Both  types  of  super-cages  allow 
for  occupation  by  small  molecules. 

The  chemical  formulas  for  the  aluminosilicates  zeolite  A  and  X  are 

Na56[(AI02)JSi02),x]x264  H20 
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and 

A/a86[(>4/O2)86(S/O2)106]x264  H2O 

respectively.  This  formula  expresses  the  similarity  between  the  two  types  with 
the  difference  of  an  altered  Al-Si  ratio,  which  leads  to  the  differences  in  the  two 
lattice  structures.  The  density  of  zeolite  X  decreases  from  1.97  g/cc  to  1.31  g/cc 
upon  dehydration.  This  is  even  more  prominent  for  type  A  with  a  reduction  from 
1 .99  g/cc  to  1 .27  g/cc.  The  void  volume  or  porosity  is  50%  for  X  and  47%  for  A. 
The  relative  surface  area  is  large  with  800  m2/g.  The  symmetry  is  cubic  in  both 
cases. 

Zeolites  can  be  synthesized.  It  is  thereby  possible  to  create  variations 
that  do  not  exist  in  nature.  The  procedure  is  to  prepare  a  highly  reactive 
aluminosilicate  gel  which  is  an  aqueous  solution  of  aluminate,  sodiumsilicate  and 
sodium  hydroxide.  The  gel  crystallizes  at  temperatures  between  room 
temperature  and  150  centigrades  under  atmospheric  or  elevated  pressures. 
Large  numbers  of  crystallite  nuclei  are  formed  which  grow  from  the 
supersaturated  gel. 

In  addition  to  the  fascinating  microscopic  structure  which  is  worthwhile 
investigating  in  itself,  zeolite  is  variably  usable104  in  every-day  applications. 

It  is  used  as  a  water  softener  exploiting  its  capability  of  exchanging  Na  for 
Ca.  The  possibility  of  custom-designing  the  width  of  connecting  channels  by 
varying  the  adsorbed  neutralizing  cations  opens  a  wide  variety  of  selective 
filtering.  A  good  example  is  the  process  of  upgrading  gasoline  by  separating 
chain  from  branching  hydrocarbons.  The  chains  slip  through  the  channels  which 
is  not  the  case  for  the  less  ignitable  branching  structure. 
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The  ionic  character  of  the  aluminosilicates  also  allows  for  selective 
separation  of  differently  polarizable  molecules.  Polar  molecules  and  saturated 
hydrocarbons  are  more  readily  adsorbed. 

Zeolites  may  even  be  used  as  carriers  of  volatile  catalysts.  The  catalysts 
are  trapped  in  the  pores  and  facilitate  the  chemical  reaction  but  are  not  lost  in 
the  process. 

Important  to  keep  in  mind  is  that  zeolites  are  very  regular  porous 
materials  with  a  wide  spectrum105  of  monodiperse  pore  sizes  connected  by 

i 

channels  of  diameters  between  3.8  and  10  A  according  to  the  chosen  type  of 
zeolite.  This  allows  the  experimenter  to  select  an  appropriate  pore  dimension 
with  the  advantage  of  monodispersity.    The  outcome  and  results  of  an 
experiment  are  therefore  more  easily  relatable  to  a  specific  pore  radius. 
4.3.2.  Vvcor 

Vycor  is  an  extensively  used  porous  substrate.  It  is  a  silica  glass 
produced  by  phase  separation  of  boron  in  boron  glass.  The  phase-separated 
boron  forms  a  connected  network  which  leaves  pores  behind  when  leached 
out.106  Vycor  is  commercially  available  (Corning  Glass  Works,  Houghton  Park, 
Corning,  N.Y.  14830)  with  a  wide  range  of  average  pore  radii.  This  makes  vycor 
a  convenient  substrate  to  use.  Interest  in  vycor,  as  a  material,  spurred  as 
experiments  revealed  a  fractal  dimension107  of  1.74  in  the  range  of  30  to  300  A. 

Average  vycor  has  a  mean  radius  of  30  A  and  a  distribution  between  20 

0 

and  55  A.  Its  porosity  is  about  28%  with  an  effective  surface  area  of  200  m2/g 
which  is  much  lower  than  for  zeolite.  The  density  is  1 .5  g/cc. 

A  drawback  compared  to  zeolite  is  that  the  pore  size  distribution  is  broad 
and  continuous.  This  makes  definite  statements  about  radial  dependence 
difficult.  Other  problems  are  reduced  interconnectivity  due  to  dead  ended  pore 
paths  which  can  act  as  dead  volume  in  special,  dynamic  experiments. 
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4.3.3.  Exfoliated  Graphite 

Exfoliated  graphite  is  the  classical  two-dimensional  substrate  with 
hexagonal  lattice  structure.  There  are  several  different  kinds  of  exfoliated 
graphites  available,  depending  on  relative  surface  area  and  z-axis  orientation  of 
the  two-dimensional  graphite  flakes.  From  grafoil  (Union  Carbide)  to  papyex 
(Carbone  Lorraine,  37-41  Rue  Jean-Jaures,  92231,  Gennevilliers)  to  UCAR  ZYX 
the  z-axis  distribution  is  narrowing  at  the  expense  of  relative  surface  area. 

Graphite  crystals  are  intercalated  with  a  strong  oxidizing  agent  such  as 
sulfuric  acid,  carefully  rinsed  and  rapidly  heated.  The  two-dimensional  character 
is  imposed  on  the  substrate  by  rolling  it  into  sheets.  The  density  is  1  g/cc  which 
represents  a  reduction  by  a  factor  of  two  by  this  process  .  The  surface  area  for 
grafoil  is  strongly  enhanced  over  graphite  with  20-30  m2/g  but  decreases  to  0.3 
m2/g  for  UCAR  ZYX.  Impurities  like  Fe,  Al,  Mn  and  Si  are  abundant  at  levels  of 
10  to  100  ppm.  Outgassing  the  sheets  by  heating  to  600  centigrades  and 
vacuum  pumping  can  reduce  these  levels  substantially,  but  impurities  still  remain 
a  problem  for  experiments  that  rely  on  clean  surfaces. 

A  severe  problem  for  high-frequency  NMR  measurements  is  the  electrical 
conductivity  in  the  graphite  planes.  For  low  temperatures  the  resistivity 
decreases  to  1.42x10"5Qm  at  low  frequencies.108  This  conductivity  would  still 
allow  for  radio-frequency  penetration  into  0.2  mm  sheets  regarding  the  skin 
depth  of  0.15  mm  at  270  MHz.  ESR  and  NMR  measurements  have  indicated, 
however,  that  the  conductivity  is  considerably  enhanced  at  frequencies  as  high 
as  270  MHz  such  that  skin  effects  expel  the  rf  from  the  substrate  and  prevent 
from  obtaining  resonance  signals  for  hydrogen  on  grafoil.109 


CHAPTER  5 
EXPERIMENTAL  DATA 

5.1     Sample  Preparation 

The  goal  is  to  prepare  a  solid  hydrogen  sample  confined  to  the  cages  of 
zeolite.  Several  steps  have  to  be  undertaken  to  achieve  this. 

Zeolite  13X  is  commercially  available  in  pellets  and  was  crushed  into  a 
coarse  powder  to  fit  into  the  sample  cell  with  dimensions  of  3  mm  in  length  and 
diameter.  This  powder  had  to  be  outgassed  in  order  to  eliminate  any  adsorbed 
water,  oxygen  and  other  impurities.  The  crystal  water  which  is  part  of  the  zeolite 
complex  is  not  affected  by  the  outgassing  procedure,  otherwise  the  zeolite  would 
"boil"  and  disintegrate.  The  outgassing  procedure  was  performed  for  10  hours  at 
about  500  centigrades  while  pumping  on  the  powder,  establishing  a  vacuum  of 
10"5Torr. 

Loading  the  sample  cell  without  readsorption  of  water  and  oxygen  was 
succeeded  at  by  performing  the  powder  transfer  in  a  dry  nitrogen  atmosphere. 
The  powder  was  taken  up  in  one  end  of  a  3  mm  hollow  pipe  and  pushed  into  the 
sample  cell  with  a  pressure  of  nitrogen  at  its  other  end  and  then  compressed 
with  the  endcap  piston  that  provided  the  vacuum  seal. 

Once  the  sample  cell  was  loaded  with  zeolite  the  vacuum  can  and  the 
dewar  were  mounted.  The  sample  cell  is  therefore  not  directly  accessible  as  is 
characteristic  for  a  low-temperature  system.   The  only  access  is  a  capillary 
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leading  from  the  external  hydrogen  gas-handling  system  through  the  cryostat  to 
the  sample  cell  at  the  bottom  of  the  dewar. 

Hydrogen  is  gaseous  at  room  temperature  and  contracts  substantially 
upon  cooling,  thereby  liquifying  at  20.4  K  and  solidifying  at  13.8  K  for  the  bulk 
system.  This  requires  a  special  technique  for  the  preparation  of  hydrogen 
samples.  Once  the  cryostat  is  cooled  to  4.2  K  with  liquid  helium,  the  sample  cell 
must  be  heated  to  above  the  melting  point  of  bulk  hydrogen,  to  about  16  K,  in 
order  to  avoid  blocking  the  hydrogen  filling  capillary.  This  is  done  by  applying  a 
voltage  to  the  hydrogen  line  heater,  which  is  wrapped  around  the  filling  capillary, 
and  adding  a  small  amount  of  helium  exchange  gas  in  the  vacuum  can  in  order 
to  establish  a  stable  equilibrium  temperature.  The  heating  power  on  the  line  and 
the  cooling  power  of  the  exchange  gas  must  balance  at  the  desired  temperature. 
Once  the  system  was  prepared  to  these  conditions  a  calibrated  amount  of 
ImMole  of  hydrogen  gas  was  introduced  into  the  capillary  by  opening  a  valve 
that  connects  the  calibration  volume  with  the  sample  volume.  The  capillary 
volume  is  hereby  neglected.  This  is  acceptable  because  most  of  the  capillary  is 
at  a  much  higher  temperature  than  the  sample  cell  and  represents  only  a  small 
volume  fraction  due  to  the  small  capillary  dimensions.  The  condensed  hydrogen 
sample  was  then  annealed  for  about  30  minutes  at  a  temperature  above  the  bulk 
solidification  point  in  order  to  reduce  crystal  imperfections.  Depending  on  the 
experimental  requirements  the  hydrogen  line  heater  was  then  reduced  at 
different  rates  and  exchange  gas  added  to  lower  the  temperature  to  equilibrium 
with  the  surrounding  liquid  helium  bath. 

Two  different  NMR  methods  at  268  MHz  were  employed  to  investigate  the 
hydrogen-zeolite  system.  Continuous  wave  spectroscopy  was  used  to  observe 
the  evolution  of  the  quasistatic  orientational  order  parameter  which  determines 
the  lineshapes  at  low  temperatures.  The  molecular  dynamics  were  probed  using 
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coherent  pulse  NMR  techniques.  Nuclear  spin  echo  amplitudes,  spin-spin  or 
transverse  relaxation  rates  T2  and  nuclear  spin-lattice  or  longitudinal  relaxation 
times  7",  were  measured  to  determine  the  spectral  densities  of  the  motion  at  low 
frequencies  and  high  frequencies,  respectively.  It  is  important  to  keep  in  mind 
that  the  total  nuclear  spin  quantum  number  for  parahydrogen  is  /  =  0  and  is 
therefore  not  detected  with  NMR. 

5.2.    Pulse  Work 

"Ordinary"  90,  -r-180,  -  r  spin  echo  techniques  were  used  to  monitor 
the  liquid  fraction  of  the  hydrogen  sample  directly.  The  dipolar  broadened  lines 
from  the  solid  component  of  the  sample  do  not  contribute  to  the  ordinary  echo, 
and  the  echo  amplitude,  after  correction  for  the  T2  decay,  provides  a  reliable 
measure  of  the  liquid  fraction,  rf  amplitudes  of  15  gauss  were  used  at  268  MHz 
for  high  sensitivity  and  the  spectrometer  recovery  time  was  35  //sec.  The 
results  of  the  pulse  data  have  been  published  in  Rail110  et  al.  mainly  focussing  on 
the  supercooling  effect. 
5.2.1.  Transverse  Relaxation  Time 
5.2.1.1.       Temperature  dependence 

The  sample  was  cooled  to  4.2  K  within  1-15  hours.  During  cooling,  the 
spin-spin  relaxation  time  T2  was  measured  by  recording  the  spin  echo  decay  as 
a  function  of  the  time  r  between  the  first  and  second  pulse  at  discrete 
temperatures.  The  relaxation  can  be  interpreted  as  two  exponential  decays  with 
a  faster  rate  initially.  This  can  be  seen  from  a  fit  of  the  echo  amplitudes  for  later 
r  whose-fit  undercuts  the  echo  amplitudes  at  early  r.  The  attempt  was  made  to 
extract  the  amplitudes  and  the  decay  rates  of  the  two  involved  functions.  The 
quality  of  the  fit  depends  on  the  initial  guesses  used.    For  a  few  distinct 
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Figure  9.  Intrinsic  transverse  relaxation  times 
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temperatures  the  deconvolution  was  performed  and  the  two  relaxation  times  are 
shown  in  figure  8.  Both  relaxation  times  display  almost  identical  dependence  on 
the  temperature.  This  is  the  reason  for  using  the  single,  average  relaxation  time 
for  further  discussion.  The  measured,  deconvoluted  values  were  used  to 
determine  the  physical  or  intrinsic  rates  from  the  heterogeneous  spin  system 
relaxation  theory  described  in  chapter  3.5.  Choosing  pbulk  =  20  %  for  the  zeolite 
system  allows  for  determination  of  T2wa",  r2bulk  and  /cb"u1lk,  the  interbath  relaxation 
time,  in  the  graphical  fashion  already  described.  The  outcome  is  shown  in  figure 
9.  One  of  the  main  results  is  that  T2b u,k  displays  a  stronger  peak  than  72wa"  which 
suggests  that  the  peak  feature  is  mainly  a  result  of  the  bulk  component.  The 
very  long  k^lk  peak  value  is  an  indication  for  negligible  interaction  between  spin 
systems  at  the  point  of  critical  behavior. 

In  the  following  discussion  average  relaxation  times  are  used  for  reasons 
mentioned  above.  The  data  can  reasonably  well  be  fit  by  a  single  exponential 
which  has  the  advantage  of  stable  fitting  results.  The  relaxation  times  show 
(Figure  10a)  a  strong  temperature  dependence  at  around  10  K.  They  peak 
sharply  at  this  temperature,  termed  the  supercooling  transition  temperature  7"^, 
displaying  similar  values  above  and  below  the  peak.  No  effect  was  observed  at 
13.8  K,  the  bulk  solidification  point.  This  behavior  has  been  interpreted  as 
supercooling,  the  suppressed  solidification  of  hydrogen  in  porous  material.  The 
observation  of  a  finite  width  peak  in  T2  instead  of  a  mere  drop  from  liquid  to  solid 
values  is  characteristic  for  the  investigated  system.  It  is  explained  in  terms  of  the 
inhomogeneity  of  the  microscopic  environment  for  different  hydrogen  molecules. 
The  width  of  the  peak  represents  therefore  a  superposition  of  critical  regions 
rather  than  one  distinct  collective  transition.  The  interpretation  of  this  data  in 
terms  of  supercooling  will  be  supported  by  further  data  in  pulse  and  continuous 
wave  mode. 
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Figure  10.  Transverse  relaxation  time  transitions 
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5.2.1.2.        Ortho-para  dependence 

The  same  measurements  have  been  performed  for  different 
orthohydrogen  concentrations  (Figure  10).  This  was  possible  by  waiting  for  the 
ortho-para  conversion  to  proceed  until  it  reached  the  desired  concentration  level. 
The  different  shapes  of  the  peaks  are  related  to  the  orthohydrogen  contents  and 
the  sample  preparation  time.  The  important  new  result  of  these  studies  is  the 
ortho  concentration  dependence  of  the  supercooling  transition  temperature  T^. 
The  supercooling  temperatures  obtained  are  between  10.8  and  8.5  K, 
suppressed  from  13.8  K  for  bulk.  The  rsc  could  be  fit  to  a  linear  function  of  the 
parahydrogen  concentration  (Figure  11).  This  had  not  been  observed  in 
previous  studies63  and  is  interesting  because  it  indicates  that  supercooling  is 
easier  to  achieve  for  parahydrogen,  the  spherically  symmetric  species.  This  is 
not  unexpected  because  of  the  additional  orbital  degeneracy  and  the  anisotropic 
interactions  of  the  orthohydrogen  molecules.  The  7"^  inferred  from  the  T2 
studies  extrapolates  to  13.4±0.5K  at  zero  parahydrogen  concentration,  and 
3.3+0.8  K  for  zero  ortho  content  with  a  slope  of  -0.10 ±0.01  K/%.  The  values 
of  for  low  ortho  contents  are  especially  interesting  because  they  are  within 
the  range  of  the  predicted  superfluid  transition  of  6  K.  If  the  entire  sample 
consisted  of  parahydrogen  and  the  relation  between  supercooling  and 
parahydrogen  contents  remained  linear  down  to  that  region,  superfluid  hydrogen 
was  a  realistic  expectation.  This  region,  however,  could  not  be  probed  using 
NMR  because  of  the  weakness  of  the  NMR  signals  for  very  low  ortho 
concentrations.  The  corresponding  ortho-para  ratios  in  the  sample  were 
obtained  from  continuous  wave  studies  discussed  below. 
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Figure  1 1 .  Supercooling  transition  as  a  function  of  para-H2  concentration 
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5.2.1.3.  Hysteresis 

It  has  been  established  by  thermally  cycling  through  the  peaks  that  the  T2 
transitions  are  non-hysteretic.  The  T2  peak  in  figure  10a  was  measured  by 
warming  and  cooling  about  the  transition  temperature  and  is  extremely  sharp.  It 
is  therefore  deduced  that  this  supercooling  is  purely  thermodynamic,  and  is 
different  from  the  dynamic  supercooling  observed  by  other  groups  for  droplets. 
For  specific  heat  measurements  of  hydrogen  in  porous  materials  a  word  of 
caution  is  in  place.  In  general,  hysteresis  is  a  phenomenon  that  describes  the 
different  behavior  of  a  system  when  it  is  alternatively  probed  upon  cooling  and 
heating.  Using  NMR  as  a  tool  of  investigation  is  unproblematic  in  this  concern 
as  it  is  not  based  on  the  application  of  heat  itself.  This  is,  however,  not  the  case 
for  a  specific  heat  measurement.  Here,  no  matter  whether  the  critical  region  is 
approached  from  above  or  below,  the  sample's  response  to  a  heating  pulse  is 
detected.  It  is  therefore  not  quite  justified  to  speak  of  probing  the  transition  upon 
cooling.  Instead,  heat  is  always  applied  and  the  quality  of  such  data  concerning 
the  delicate,  entirely  heat-dependent  property  of  hysteresis  is  not 
unquestionable.  In  the  present  case  of  NMR  data  there  is  sample  dependent 
sharpness  of  the  T2  peak  and  the  values  of  T2,  but  this  is  not  a  result  of  thermal 
cycling. 

Another  strong  argument  against  hysteresis  will  be  given  using  the  7*2 
data  of  hydrogen-deuteride  in  zeolite. 
5.2.2.  Spin  Echo  Amplitude 

Not  only  the  spin-spin  relaxation  times  T2  display  a  characteristic  behavior 
around  10  K  but  also  the  recorded  NMR  echo  amplitudes.  The  general  behavior 
is  a  steep  increase  around  11  K  forming  a  distinct  peak  at  about  10  K  and  then 
dropping  sharply  by  a  factor  of  20.  The  latter  change  is  also  attributed  to  the 
liquid/solid  transition.    This  behavior  has  been  observed  for  many  different 
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samples  (Figure  12)  which  were  prepared  under  slightly  different,  but  controlled 
conditions,  to  be  discussed  below.  It  was  also  the  initial  indication  for 
supercooling  in  the  hydrogen-zeolite  system. 

The  recorded  echo  amplitudes  must  be  corrected  for  the  Curie  law 
according  to  which  the  magnetization  in  the  high  temperature  limit  is  proportional 
to  the  inverse  temperature.  The  echoes  are  therefore  multiplied  by  the 
temperature  at  which  they  were  taken.  It  should  be  remarked  that  the  initial 
change  in  measured  echo  amplitudes  is  due  to  the  occurrence  of  the  dramatic 
peaking  of  the  nuclear  spin-spin  relaxation  time  T2.  The  final  drop  of  the  echo 
amplitudes  upon  cooling  is  then  a  consequence  of  the  diminished  liquid  fraction 
of  the  hydrogen  sample.  The  average  relaxation  times  were  shown  in  figure  10 
for  the  two  samples,  but  two  components  could  be  distinguished  which  are  also 
attributed  to  the  narrow  and  broad  components  of  the  cw  NMR  line  shape  to  be 
discussed  later. 

Although  the  general  features  are  the  same  for  different  samples  the 
shape  of  this  transition  varies.  The  supercooling  temperature  Tsc  for  the  echo 
amplitude  falls  in  the  range  of  9.5  to  12.4  K.  The  width  and  relative  heights  of 
the  transition  are  equally  sample  dependent.  The  varying  behavior  depends  on 
several  conditions:  the  sample  preparation  which  includes  the  temperature  at 
which  the  hydrogen  gas  was  condensed  and  the  cooling  rate  after  condensation, 
the  orthohydrogen  concentration  and  whether  the  sample  was  warmed  or  cooled 
while  the  echo  amplitudes  were  measured.  Two  effects  are  observed.  Firstly, 
the  sharpness  of  the  peak  decreases  as  the  cooling  time  increases  as  seen  in 
figures  12a,  12b,  12d  and  12e.  The  peak  becomes  broader  and  the  relative 
height  decreases.  The  discontinuity  upon  approaching  the  peak  from  high 
temperatures  also  disappears  in  this  sequence.  Secondly,  a  weak  connection 
seems  to   exist   between   a   longer  cooling   time   and   higher  as 
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demonstrated  by  figures  1 2a  at  9.5  K,  1 2c  at  1 0  K,  1 2d  at  1 0.7  K  and  1 2e  at  1 2.4 
K.  Slower  cooling  allows  the  sample  to  form  more  homogeneously  which  leads 
to  a  less  suppressed  solidification.  The  warming  transitions  in  figure  12f  and  12g 
show  a  less  prominent  peak.  This  is  assumed  to  be  the  case  because  parts  of 
the  hydrogen  sample  leave  the  cell  upon  warming.  An  additional  reason  for  the 
different  shapes  is  that  the  sample  preparation,  the  temperature  and  the 
pressure  of  condensation,  and  the  filling  of  the  pores  are  never  precisely 
reproducible. 

Two  exceptional  features  could  be  observed.  Figure  12b  shows  the  bulk 
liquid/solid  transition  which  is  measured  at  13.8  K  by  the  drop  in  echo  amplitude. 
The  explanation  for  this  effect  is  that  bulk  hydrogen  was  still  outside  the  zeolite. 
This  was  the  only  time  the  pure  bulk  feature  occurred.  The  other  feature  is 
shown  in  figure  12e.  Instead  of  a  single  peak  the  feature  rather  consists  of  two 
or  even  three  peaks.  A  speculation  is  that  the  pores  (12  A)  exhibit  a  different 
transition  temperature  than  the  channels  (9  A)  in  zeolite.  The  transition  occurs 
slowly  enough  to  distinguish  between  the  two  different  diameters.  No  other 
experimental  data,  however,  lead  to  such  a  conclusion  which  makes  this  a 
tentative  explanation. 

A  more  sophisticated  analysis  removes  the  T2  effect  on  the  echoes  by 
extrapolating  the  echo  amplitudes  back  to  "zero"  time.  The  fully  corrected  echo 
amplitudes  are  proportional  to  the  liquid  component.  Figure  13  shows  7*2  and 
the  fully  corrected  echoes  taken  at  an  ortho  concentration  of  73  %.  It  is  obvious 
how  the  T2  peak  coincides  with  the  drop  of  the  corrected  echo  amplitude.  This 
leads  to  the  conclusion  that  at  this  temperature  the  solidification  takes  place. 
5.2.3.  Longitudinal  Relaxation 

For  more  information  about  the  dynamics  of  the  molecules  the  spin-lattice 
relaxation  time  7,  was  also  measured.  The  chosen  method  to  extract  7,  was  to 
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vary  the  repetition  time  of  the  set  of  double  pulses.  The  idea  is  that  the  spins 
need  time  to  recover  back  to  the  z-direction  and  less  magnetization  will  be 
measured  when  the  repetition  time  is  short.  This  can  again  be  monitored  by  the 
resulting  spin  echo  amplitudes  and  fit  to  a  relaxation  time  assuming  exponential 
behavior.    The  problem  with  this  method  are  the  varying  amounts  of  heat 
introduced  into  the  spin  system  by  the  changing  repetition  times  as  the  echoes 
are  scanned  for  one  relaxation  time.  The  data  (Figure  14)  are  not  of  high  quality 
due  to  these  problems  with  the  data  analysis  but  do  not  show  real  evidence  for  a 
transition  over  the  temperature  range  3  <  T  <  15  K.  The  7j  data  for  the  cooling 
process  (Figure  14a)  seem  to  show  a  peak  at  13  K  and  a  dip  at  10.5  K.  This 
could,  however,  not  be  confirmed  upon  warming  (Figure  14b).  The  result  is  that 
the  relaxation  time  decreases  overall  with  lower  temperature  with  a  substantial 
scatter,  not  clearly  exhibiting  any  sharp  feature  at  10  or  13.8  K,  the  temperatures 
in  question.  This  leads  to  the  conclusion  that  the  dynamics  of  the  system  is 
mostly  governed  by  the  zero-frequency  regime  where  T2  has  contributions  from 
the  spectral  densities  but  not  7j.  As  the  sample  is  cooled  the  spectral  density 
changes  as  a  function  of  frequency.  The  area  under  the  spectral  density  is 
conserved  so  that  an  additional  contribution  at  lower  frequencies  must  be  taken 
from  higher  frequencies  as  the  temperature  decreases.  The  peak  in  T2  is  a 
consequence  of  the  zero  frequency  contribution.  Although  the  spectral  density 
at  low  frequencies  rises,  it  must,  however,  dip  at  zero  frequency  in  order  to 
explain  the  observed  behavior.  While  the  motion  of  larger  clusters  close  to  the 
critical  region  shifts  towards  lower  frequencies  the  zero-frequency  contribution 
even  diminishes  because  no  motion  at  all  is  only  favorable  in  the  solid  state. 
This  is  the  reason  why  T2  drops  again  below  the  critical  region  in  the  solid.  The 
observation  that  the  7,  data  does  not  show  evidence  of  a  phase  transition 
supports  the  interpretation  of  the  peak  in  Tz  and  the  spin-echo  amplitudes  in 
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terms  of  the  solidification  of  the  liquid  component  rather  than  an  abrupt  change 
in  the  spectral  density  of  the  system. 

5.3.    Continuous  Wave  Work 

cw  lineshapes  were  recorded  to  evaluate  the  quasistatic  orientational 
order  parameter  a.  The  orientational  ordering  of  a  small  number  of  quantum 
rotors  with  highly  frustrated  interactions  raises  a  number  of  questions  concerning 
the  nature  of  the  phase  transitions  and  type  of  order  achieved.  Answers  are 
expected  from  the  shape  and  the  time-  and  temperature  dependence  of  the 
continuous  wave  NMR  line.  The  results  are  published  in  Rail111  et  al.  with 
emphasis  on  the  anomalous  ortho  para  conversion. 
5.3.1.  Ortho-Para  Conversion 

It  was  initially  tried  to  track  the  orthohydrogen  concentration  by  measuring 
the  spin  echo  amplitude  overtime.  This,  however,  turned  out  to  be  an  erroneous 
method  under  the  given  circumstances  where  the  line  was  too  wide  for  the 
pulses  to  cover  the  entire  spectrum  satisfactorily.  Instead  of  measuring  the 
changing  ortho  contents  an  additional  change  in  the  lineshape  constituents  was 
folded  in.  This  is  verified  by  the  changing  spin-spin  relaxation  time  over  the 
course  of  time  (Figure  15).  The  lineshape  shifts  its  relative  weights,  where  the 
wings  gain  more  importance  relative  to  the  center  component  as  orthohydrogen 
decays  away.  The  echo  amplitude  is  a  slightly  oscillatory  function  of  the  time 
between  the  two  applied  pulses  with  a  frequency  of  3  kHz  (Figure  16).  This 
nonexponential  behavior  is  possibly  a  consequence  of  the  different  degrees  of 
orientational  ordering  in  the  system. 

Continuous  wave  NMR  was  therefore  used  to  track  the  orthohydrogen 
concentration.  The  derivative  of  the  NMR  absorption  lineshapes  is  recorded 
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using  standard  magnetic  field  modulation  and  lock-in  techniques.  The  recorded 
information  with  such  a  setup  is  the  differential  lineshape.  It  is  integrated  in 
order  to  obtain  the  lineshape  itself.  As  the  area  of  the  lineshape  is  a  direct 
measure  of  the  susceptibility  of  the  sample  and  therefore  also  of  the  number  of 
orthohydrogen  molecules  it  is  integrated  a  second  time.  The  differential 
lineshape  is  therefore  doubly  integrated  to  determine  the  total  nuclear 
magnetization  and  thus  the  orthohydrogen  concentration.  The  way  this 
integration  is  performed  is  crucial  for  the  outcome  of  the  quality  of  the  data 
analysis.  In  the  ideal  case,  the  differential  lineshape  is  perfectly  symmetric  about 
its  inflection  point  in  terms  of  amplitude  and  width.  Any  deviation  in  this 
symmetry  affects  the  integration.  It  is  therefore  important  to  incorporate  a 
measure  of  symmetry  in  the  integration  routine.  This  was  solved  by  writing  a 
routine  that  integrated  each  half  of  the  line  separately  and  compared  the  results 
for  the  two  sides.  Another  advantage  of  this  method  is  that  the  baseline 
correction  is  performed  individually  for  both  sides  and  also  compared  with  each 
other  as  a  further  measure  of  data  quality.  It  is  thereby  guaranteed  that  no 
spurious  baseline  contributions  enter  into  the  integration. 

The  orthohydrogen  concentration  was  measured  over  time  for  three 
different  samples.  The  data  were  taken  at  4.2  K  up  to  1550  hours.  After 
carefully  stabilizing  the  NMR  spectrometer  by  regulating  the  temperature  of  the 
low-level  amplifier  to  ±0.1°C,  a  1550  hour  continuous  experiment  was 
performed.  The  result  is  shown  in  figure  17a.  The  data  could  be  fit  with  the 
solution  of  a  bimolecular  process 
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where  S  is  the  signal,  x(0)  =  75%,  which  is  the  initial  orthohydrogen  contents  at 
room  temperature,  and  t  the  time.  For  the  first  500  hours  the  conversion  rate 
constant  was  determined  to  be  k,  =0.425 ±0.006% /h  and  for  the  remaining 
time  k2  =2.21  ±0.075% /h.  Figure  17b  shows  the  inverse  conversion  data 
which  are  fit  with  two  straight  lines.  The  relative  error  increases  for  small  x 
because  the  signal  decreases  while  the  absolute  error  remains  the  same. 

The  initial  conversion  rate  is  anomalously  slow  compared  to  the  bulk 
value  of  1 .9%/h  at  4.2  K  by  a  factor  of  approximately  4.  Equally  remarkable  is 
the  increase  to  an  even  larger  value  for  long  times. 

In  order  to  confirm  these  results,  two  further  samples  were  investigated  to 
exclude  amplifier  drifts  and  other  destabilizing  influences.  The  conversion  rate 
constants  were  kg  =0.35 ±0.08% /h  and  kc  =0.54±0.13%/h.  The  errors  in 
these  latter  experiments  were  larger  due  to  the  short  measurement  time  but 
entirely  consistent  with  the  short-term  conversion  of  the  extensively  studied 
sample.  This  test  was  important  to  establish  the  stability  of  the  NMR  apparatus 
which  was  thermally  regulated  and  the  gain  monitored  carefully  over  the  duration 
of  the  experiments.  The  use  of  an  independent  standard  proton  NMR  sample  in 
order  to  monitor  the  overall  gain  of  the  spectrometer  was  impractical  at  the  UHF 
frequencies  used  for  high  sensitivity  studies,  and  proved  to  be  unnecessary. 

Assuming  a  binomial  distribution  Wn  for  orthohydrogen  in  an  equilibrium 
ortho-para  mixture  the  mean  number  of  orthohydrogen  nearest  neighbors 
surrounding  an  ortho  molecule  is  given  by 

M  =  f^nWn=Nx 

where  x  is  the  orthohydrogen  concentration  and  N  the  number  of  nearest 
neighbors.  As  long  as  the  self-diffusion  in  the  sample  is  substantial  it  is  capable 


122 


of  homogenizing  the  distribution  which  otherwise  deviates  from  equilibrium  as  a 
result  of  conversion.  An  explanation  for  the  slow  initial  conversion  rate  is  the 
reduced  number  of  nearest  neighbors  for  the  molecules  on  the  wall  which 
constitute  about  80%  of  the  sample.  A  simple  calculation  yields  a  reduced 
average  number  of  nearest  neighbors.  It  is  assumed  that  about  20%  of 
molecules  are  in  the  center  of  the  pores  having  the  standard  12  nearest 
neighbors  and  the  remaining  80%  having  only  5  nearest  neighbors  due  to 
contact  with  the  curved  wall  which  represents  less  than  the  half-space.  This 
leads  to  an  average  number  of  nearest  neighbors  of  6.4.  This  is  more  than  half 
the  value  of  bulk.  This  would  explain  a  drop  of  the  conversion  rate  by  a  factor  of 
3.5,  since  the  rate  varies  as  the  square  of  the  number  of  nearest  neighbors. 

The  conversion  rate  is  also  expected  to  depend  on  the  hydrogen-wall 
adsorption  energy  which  will  be  shown  to  be  about  270  ±50  K  using  pressure 
versus  temperature  data  when  the  sample  was  taken  out.  This  value  for  the 
absorption  energy  is  smaller  but  comparable  to  the  value  of  400  ±100  K  inferred 
by  Monod  et  al.112  The  adsorption  potential  inhibits  the  mobility  and  the 
fluctuations  of  the  molecules  and  reduces  the  catalytic  influence  at  the  wall. 

The  conservation  of  energy  during  the  conversion  requires  that  the 
excess  rotational  kinetic  energy  (2B,  =  171  K)  be  carried  away  by  phonons.  The 
Debye  energy  for  bulk  hydrogen  at  zero  pressure  is  120  K,  and  this  requires  a 
two-phonon  process  for  conversion  in  the  bulk.  One  expects  that  the  surface 
binding  reduces  the  density  of  surface-phonon  states  and  thus  slow  the 
conversion  rate  involving  a  molecule  near  the  wall.  Theoretical  estimates  of  the 
phonon  states  have  been  carried  out,  and  show  that  the  mechanism  is  still  a  two- 
phonon  process.  The  density  of  phonon  states  is  necessary  for  a  quantitative 
comparison  with  the  observations.  The  simple  dimensional  consideration  of  2 
versus  3  phonon  modes  on  the  wall  reduces  the  density  of  states  by  this  ratio.  It 
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is  intuitively  clear  that  the  tortuosity  of  the  porous  material  also  restricts  the  long- 
wavelength  phonons  reducing  their  effect  of  carrying  away  energy.  These 
additional  influences  are  difficult  to  express  quantitatively  but  lead  to  a  further 
reduction  in  the  conversion  rate  coming  close  to  the  observed  factor  of  4. 

While  a  bimolecular  process  ("B"  in  Figure  17c,d)  appears  to  provide  the 
best  fit  at  later  times,  the  relative  experimental  errors  are  larger  at  lower 
concentrations  and  an  exponential  conversion  could  fit  the  data.  Such  a  fit,  as 
shown  in  figure  1 7c,  is,  however,  inconsistent  with  the  behavior  at  early  times 
because  it  would  numerically  dominate  the  bimolecular  rate  at  all  times  because 
any  impurity-induced  conversion  would  be  constant  in  time  and  not  be  "switched- 
on"  at  t  =  500  hours.  The  decay  represented  by  the  dashed  line  A1  of  figure  17c 
corresponds  to  an  exponential  decay 

x(0  =  x(0)exp(-r/r1) 

with  a  decay  constant  z,  =340±10  h  for  T  <  500  hours.  As  shown  by  the  solid 
line  A2  of  figure  17d,  the  decay  for  late  times  can  be  fit  by  an  exponential  decay 
with  a  time  constant  r2  =  573±15  h. 

The  behavior  observed  at  late  times  is  similar  to  the  Schmidt  result  for 
bulk  hydrogen  for  lower  temperatures  (1.57  K).  He  found  that  diffusion  is  large 
enough  at  4  <  T  <12  K  to  sustain  the  same  conversion  constant  while  for  1.57  K 
the  quadrupole  effect  overcompensates  the  reduced  diffusion  and  leads  to 
clustering,  thus  accelerating  the  conversion.  The  enhanced  conversion  after  500 
hours  must  therefore  be  explained  by  the  fact  that  the  mean  number  of  nearest 
neighbors  is  larger  than  M0=Nx.  This  may  be  due  to  the  formation  of  small 
clusters  whose  existence  has  been  inferred  from  the  structure  of  the  cw  NMR 
lineshape  and  its  corresponding  quasistatic  order  parameter.    The  quantum 
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diffusion  of  ortho  molecules  in  the  interior  of  the  zeolite  pores  will,  however,  be 
determined  by  the  dipolar  interactions  as  it  is  in  bulk  solid  hydrogen.  The  overall 
influence  of  the  wall  potential  is  therefore  expected  to  raise  the  temperature  for 
which  clustering  in  pairs  or  larger  units  occurs  with  respect  to  the  bulk  because 
of  the  stronger  anisotropic  potentials  at  the  surface  and  the  preferential 
orientation  of  the  ortho  molecules  at  the  surface  sites.  This  has  been  observed 
for  the  temperature  dependence  of  the  line  width  upon  cooling.  It  can  explain 
why  clustering  effects  are  observed  at  temperatures  as  high  as  4  K  in  these 
geometries  compared  to  an  upper  limit  of  1.5  K  for  bulk  hydrogen.  For  lower 
orthohydrogen  concentration,  ortho  molecules  can  diffuse  more  easily  as  the 
condition  for  energy  conservation  for  the  process  of  tunneling  is  satisfied  more 
probably.3  This  allows  for  ortho  cluster  formation  at  lower  ortho  concentrations 
and  later  times.  Clusters  of  pairs  and  triplets  correspond  to  a  mean  number  of 
nearest  neighbors  of  2  and  3.  This  is  of  negligible  influence  for  an  original  mean 
value  of  about  6  orthohydrogen  molecules.  As  soon  as  the  mean  number  drops 
to  2  the  clusters  start  to  contribute  an  above-average  value.  From  then  on  the 
conversion  rate  constant  increases  although  the  absolute  number  of  converting 
molecules  is  still  smaller  than  for  earlier  times. 
5.3.1.1.  Impurities 

The  attempt  to  fit  the  orthohydrogen  concentration  (Figures  17c  and  17d) 
using  a  single  exponential  fit  did  not  yield  satisfactory  results.  A  concentration  of 
approximately  100  ppm  of  paramagnetic  impurities  in  the  zeolite  would  be 
needed  to  generate  surface  conversion  rates  comparable  to  the  overall  rates 
observed.  This  estimate  was  obtained  from  a  direct  comparison  of  the  studies  of 
Cunningham  and  Johnston80  of  the  conversion  of  hydrogen  on  paramagnetic 
impurities  adsorbed  on  porous  material.  For  Fe  impurity  concentrations 
xFa  =5*10"2  on  known  surface  areas,  they  observed  decays  of  kFg  =1.7%  per 
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sec,  while  the  present  value  is  /f  =  0.4%  per  hour.  This  leads  therefore  to  an 
estimate  of  the  required  paramagnetic  impurity  concentrations 


^  para  ~  *Fe 


^Fe    =130  ppm 


where  iiFe  =5.92/is  is  the  magnetic  moment  of  Fe  and  fj.g  the  Bohr  magneton. 
This  indicates  that  magnetic  impurities  in  zeolite  do  not  play  an  important  role  in 
the  conversion  process,  and  the  sample  preparation  methods  were  successful  in 
keeping  the  surface  free  from  contamination. 
5.3.2.  Lineshape 

As  discussed  in  chapter  3  about  NMR,  diatomic  molecules  give  rise  to  a 
special  line  shape  pattern.  For  axial  quadrupoles  the  order  parameter  is  given 
by  g  =  -±(3J]  -  J2).  A  single  crystal  sample  would  result  in  two  NMR 
absorption  lines  at  q)l±Aq)(P,g)  with  /?  the  angle  between  magnetic  field  and 
the  local  molecular  symmetry  axis,  coL  being  the  Larmor  frequency.  For  a 
powder  sample  this  leads  to  a  Pake  doublet  with  sharp  features.  If  in  addition  a 
does  not  have  a  fixed  value,  the  line  shape  becomes  a  smooth  bell  shape.  For  a 
frustrated  system  like  the  one  in  question  it  is  highly  likely  to  observe  lineshapes 
that  display  some  type  of  glassy  behavior.  The  width  and  the  specific  shape 
lead  to  knowledge  about  the  degree  of  ordering  and  frustration.  The  recorded 
lineshapes  are  characteristic  for  glassiness  and  have  a  surprisingly  large  width 
of  about  130  kHz  at  4.2  K  (Figure  18).  There  is  a  dominant  center  component 
with  a  width  at  half  height  of  25  kHz  and  a  wing  component  extending  out  to  the 
full  line  width.  The  time  and  temperature  dependence  of  the  lineshape  give 
further  information  about  the  nature  and  origin  of  the  lineshape. 
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5.3.3.  Time  Dependence 

The  lineshape  was  followed  as  the  sample  aged  (Figure  18),  and  it  was 
observed  that  the  narrow,  central  part  of  the  resonance  line,  associated  with  a 
disordered  contribution,  with  orientationai  order  parameter  <j  =  ±(3J*  -  J2)  =»0, 
decays  faster  than  the  wings.  This  can  be  seen  by  a  comparison  of  the  decaying 
lineshapes.  The  relative  weight  of  the  narrow  center  component  is  fading 
because  of  the  faster  conversion  in  the  bulk  environment  with  more  nearest 
neighbors.  This  is  entirely  consistent  with  the  previous  interpretation  for  the 
conversion  mechanisms. 
5.3.3.1         Second  Moment 

The  total  second  moment  M2  saturates  after  500  hours,  indicating  that 
there  is  no  further  increase  in  Mfn  associated  with  increased  orientationai 
ordering,  i.e.  higher  values  of  o.  This  confirms  the  explanation  given  above  in 
terms  of  the  formation  of  clusters  at  late  times,  for  which  a  is  fixed  with 
cravg  =0.45.  This  was  determined  using  a  computer  simulation  (Chapter  6) 
superimposing  Pake  doublets,  corresponding  to  different  order  parameters. 

Figure  19  shows  the  behavior  and  the  large  values  of  the  second  moment 
over  time.  The  values  are  scattered  due  to  the  uncertainty  introduced  by  the 
asymmetry  of  the  differential  line  shape  and  the  numerical  double  integration. 

5.3.4.  Temperature  Dependence 

The  cw  line  shapes  evolve  smoothly  from  0.94  K  to  15.8  K  (Figure  20). 
There  is  no  sharp  feature  at  the  supercooling  transition  temperature  associated 
with  a  change  in  the  quasistatic  order  parameter.  The  initially  dominant  center 
component  decreases  in  weight  with  increasing  temperature.  For  high 
temperatures  only  the  broad  contribution  remains.  Interesting  in  this  context  is 
the  remaining  width  at  high  temperatures  of  15.8  K  (Figure  20j).  The  width  is  still 
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Figure  19.  Second  moment  as  a  function  of  time 
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governed  by  the  wing  component  while  the  center  component  has  disappeared. 
From  the  pulse  data  it  has  been  concluded  that  hydrogen  is  liquid  at  these  high 
temperatures,  which  are  even  above  the  bulk  melting  temperature.  This  implies 
that  the  liquid  is  strongly  constrained  by  the  wall  potential  and  is  inhibited  in  its 
motion  as  is  expressed  by  the  wide  line  which  is  not  motionally  narrowed  at  such 
high  thermal  excitation. 

The  nonexistence  of  a  sharp  transition  around  10  K  (Figures  20e,  20f, 
20g,  20h)  and  the  change  in  the  nature  of  the  line  shape  implies  that  the  narrow 
lineshape  component  of  the  hydrogen  sample  leaves  the  cell  first  upon  warming. 
This  is  related  to  the  contribution  in  the  center  of  the  pores  which  is  almost 
unaffected  by  the  wall  potential.  The  remaining  hydrogen  is  a  restricted  liquid 
bound  to  the  surface.  This  again  points  to  a  two-component  system  of  two  types 
of  environments  for  the  adsorbed  hydrogen. 

Using  a  dilution  refrigerator,  the  sample  was  cooled  to  lower 
temperatures.  The  question  arose  whether  the  system  would  orientationally 
order.  This  represents  a  crucial  test  whether  the  observed  transition  in  the  spin- 
spin  relaxation  time  is  supercooling  or  a  severely  raised  onset  temperature  of 
orientational  ordering.  For  the  bulk  system  the  onset  for  ordering  is  related  to 
the  quadrupolar  interaction  constant  of  r  =  0.8  K  and  is  naturally  dependent  on 
the  orthohydrogen  contents.  There  is  a  possibility  that  this  temperature  shifts 
up,  aided  by  the  surface  potential  that  locks  the  molecules  in  place  relative  to 
each  other  and  the  rough  substrate  wall.  This  potential  has  an  anisotropic 
component  which  preferentially  binds  ortho  molecules.  Clustering  and  ordering 
is  expected  to  occur  predominantly  at  the  zeolite  surface  because  of  its  strong 
binding  potential.  The  overall  influence  of  the  wall  potential  is  therefore  expected 
to  raise  the  temperature  for  which  clustering  in  pairs  or  larger  units  occurs  with 
respect  to  the  bulk  because  of  the  stronger  anisotropic  potential  at  the  surface 
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Figure  20.  --  continued 
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and  the  preferential  orientation  of  the  ortho  molecules  at  the  surface  sites. 

The  occurrence  of  two  conversion  constants  confirms  this  expectation.  In 
analogy  with  the  bulk  data  at  1.57  K  and  the  explanation  of  the  fast  rate  at  later 
times  in  terms  of  clustering  at  4.2  K  some  orientational  ordering  is  present  at 
higher  temperatures  than  for  bulk. 

This  effect  is,  however,  unlikely  to  be  the  reason  for  the  transition  at  10  K 
because  the  cw  lineshapes  do  not  exhibit  any  change  in  width  over  this 
temperature  range.  Indeed,  the  NMR  lines  broaden  from  130  kHz  to  200  kHz 
upon  cooling  from  4  Kto  0.94  K  (Figure  20a,  20b),  which  indicates  that  additional 
ordering  takes  place  at  these  characteristic  bulk  ordering  temperatures.  This 
implies  that  the  surface  supports  the  ordering,  but  not  substantially,  and  that  the 
transitions  observed  in  pulse  mode  are  not  due  to  orientational  ordering.  This  is 
a  strong  indication  for  supercooling  at  10  K. 

The  overall  temperature  behavior  leads  to  the  interpretation  of  a  two- 
component  system.  The  lineshapes  observed  at  4  K  are  found  to  consist  of  two 
components:  (1)  a  relatively  narrow,  central  (liquid-like)  component 
corresponding  to  approximately  20%  of  the  spins  which  is  the  bulk  disordered 
contribution  with  a  =  0  in  the  center  of  the  pores;  (2)  a  broad  line,  130  kHz  in 
width,  originating  from  small  clusters  in  contact  with  the  surface  walls.  This  is 
associated  with  a  broad  distribution  of  Pake  doublets.  It  is  deduced  from  a 
computer  simulation  for  the  Pake  doublet  superposition,  discussed  in  the  next 
chapter,  that  the  order  parameters  a  have  a  distribution  around  a  mean  of  about 
0.45  with  a  width  of  0.3.  In  this  analysis  it  is  found  that  the  outer  wings  of  the 
Pake  doublet  superposition  were  below  noise  levels  and  not  observable.  The 
interpretation  of  this  cw  NMR  result  is  that  at  least  a  certain  percentage  of  ortho 
molecules  are  partially  ordered  as  small  clusters.  The  order  parameters  for  an 
orthohydrogen  pair  is  a  =  0.25  and  for  triplets  cr  =  0.38.    Complete  ordering 
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occurs  for  a  value  of  unity.  The  roughness  and  inhomogeneity  of  the  zeolite 
surface  are  expected  to  lead  to  such  a  distribution  of  a  for  partially  ordered 
states. 

5.4.  Hydrogen-Deuteride 

The  idea  for  using  hydrogen-deuteride  (HD)  is  that  it  is  very  similar  to 
parahydrogen.  In  HD,  one  hydrogen  atom  is  substituted  by  one  neavy  hydrogen 
atom,  the  deuteron.  HD  has  slightly  more  mass  than  hydrogen  and  therefore 
slightly  less  zero-point  motion.  It  otherwise  behaves  very  similar  to 
parahydrogen  because  its  groundstate  is  also  J  =  0  and  spherically  symmetric. 
There  is  only  one  type  of  HD  because  it  is  heteronuclear,  unlike  in  the  case  of 
hydrogen.  This  is  also  the  reason  for  an  exceptional  property  of  HD:  its  center 
of  mass  and  center  of  rotation  do  not  coincide. 

The  replacement  of  hydrogen  by  HD  has  the  advantage  that 
parahydrogen  can  be  simulated.  This  is  of  special  interest  in  the  case  where  it  is 
necessary  to  distinguish  between  orientational  and  translational  properties.  HD 
does  not  have  orientational  degrees  of  freedom  because  it  is  spherically 
symmetric.  Orientational  ordering  can  therefore  not  take  place.  This  is  another 
valuable  test  to  distinguish  between  orientational  ordering  and  supercooling  at 
the  observed  1 0  K. 

Four  sets  of  data  were  taken  for  HD  in  zeolite  which  overlap  with  the  work 
on  hydrogen  in  zeolite. 

The  longitudinal  relaxation  time  7,  increases  for  lower  temperatures  when 
the  sample  is  cooled  (Figure  21)  or  warmed  and  is  much  longer  than  for 
hydrogen  because  of  the  missing  dipolar  interaction.  Several  data  points  were 
taken  in  the  critical  region  but  no  discontinuity  or  peak  were  observed.  Similar 
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Figure  21 .  Longitudinal  relaxation  time  of  hydrogen-deuteride 
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Figure  22.  Echo  amplitudes  of  hydrogen-deuteride 
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Figure  23.  Transverse  relaxation  times  of  hydrogen-deuteride 
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caution  as  in  the  case  of  hydrogen  applies  to  this  data  . 

The  corrected  echo  amplitudes  decrease  weakly  upon  cooling  up  to  10  K 
and  collapse  to  small  values  thereafter  (Figure  22).  The  identical  behavior  was 
observed  upon  warming. 

The  T2  data  recorded  over  the  interesting  temperature  range  show  similar 
behavior  as  for  hydrogen  (Figure  23).  Again  a  peak  is  observable  at  10.5  K  for 
cooling  (Figure  23a)  and  9.8  K  for  warming  (Figure  23b).  The  peak  is  only  steep 
on  the  low  temperature  side  and  less  distinct  for  temperatures  above  the  peak. 
The  difference  in  transition  temperatures  for  cooling  and  heating  is  small  and 
can  hardly  be  made  responsible  for  hysteresis.  This  is  especially  so  because 
the  hysteresis  should  lead  to  a  lower  transition  temperature  for  the  cooling 
procedure  rather  than  a  higher  temperature.  This  is  another  powerful  argument 
against  hysteresis. 

The  continuous  wave  line  is  of  similar  shape  as  for  hydrogen  with  an 
increased  width  of  170  kHz. 

The  convincing  outcome  of  this  experiment  is  that  the  observed  effect  is 
indeed  supercooling,  for  HD  and  hydrogen.  Worth  mentioning  is  the  fact  that  the 
suppression  of  the  melting  point  is  quite  substantial  in  the  case  of  HD  because 
the  bulk  melting  point  is  18  K  which  represents  supercooling  by  an  entire  8  K. 
This  is  an  extension  of  figure  11  -  the  supercooling  transition  temperature  as  a 
function  of  parahydrogen  -  because  HD  simulates  100  %  parahydrogen. 

5.5     Isochoric  Pressure  versus  Temperature  Data 

Valuable  information  can  be  extracted  from  heating  the  sample  from  the 
solid  to  the  gaseous  state  while  keeping  the  volume  constant.  The  way  the 
pressure  rises  is  characteristic  for  the  system.   This  procedure  is  performed 
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when  all  measurements  on  the  sample  are  completed  and  it  is  the  last  step 
towards  taking  out  the  sample. 

The  hydrogen  line  heater  is  gradually  increased  so  that  the  temperature 
and  the  corresponding  pressure,  measured  at  the  filling  line  on  top  of  the 
cryostat,  are  always  in  equilibrium.  In  practice  the  isochoric  condition  is  not 
perfectly  satisfied.  When  the  sample  is  heated  the  sample  cell  volume  becomes 
connected  with  the  capillary,  which  had  been  disconnected  by  a  plug  of  solid 
hydrogen.  The  capillary  volume  is,  however,  small  compared  to  the  sample 
volume.  There  are  four  distinct  regions  of  pressure  observable,  slightly  shifted 
for  different  samples,  as  the  temperature  rises  (Figure  24).  Initially  there  is  no 
response  in  pressure  until  about  8  K  where  the  first  bulk  material  in  the  hotter 
capillary  and  the  hydrogen  in  the  center  of  the  pores  liquifies  and  raises  the 
pressure.  This  continues  until  about  12  K  from  where  on  a  mere  expansion 
follows.  At  this  point  no  more  molecules  contribute  to  the  pressure  but  according 
to  the  ideal  gas  law  the  increasing  temperature  raises  the  pressure. 

The  third  region,  starting  at  about  18  K,  is  again  characterized  by  a  steep 
rise  in  pressure.  This  is  attributed  to  the  desorption  of  hydrogen  molecules  off 
the  substrate  walls.  It  was  attempted  to  fit  this  increase  by  an  exponential 

(  e  \ 
p  =  p0exp  — js. 

V    Kl  J 


where  Eads  is  the  adsorption  energy  of  hydrogen  on  the  zeolite  walls  which  is 
extracted  from  the  rise  of  the  pressure  curve.  There  are  problems  with  this 
analysis  in  addition  to  the  possibility  of  nonequilibrium  data  points.  Firstly,  it  is 
not  clear  where  to  set  the  baseline  of  the  exponential  including  the  problem  of 
the  linear,  ideal  gas  law  increase  and  the  choice  of  the  origin  for  the  pressure 
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Figure  24.  Isochoric  pressure  versus  temperature  data 
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and  temperature.  Secondly,  as  the  pressure  rises  the  number  of  molecules 
changes  continuously.  This  has  an  impact  on  the  pressure  constant  p0,  which  is 
neglected.  The  results  for  the  three  different  samples  are  204,  289  and  329  K 
with  errors  of  several  Kelvin.  This  results  in  a  value  of  the  adsorption  energy 
E  ^  =270  ±50  K  and  is  lower  than  an  earlier  estimate  of  400  K.  Compared  to 

3QS 

graphite  with  a  strong  adsorption  potential  of  500  K  the  higher  value  for  zeolite 
seems  especially  high. 

The  last  region  is  again  a  linear  thermal  expansion  region  where  all 
molecules  are  gaseous,  setting  in  at  22  K.  The  expansion  was  followed  up  to  50 
K  with  pressures  of  1200  Torr  and  demonstrated  the  high  quality  of  the  heat 
shrinking  vacuum  seal  of  Kel-F  on  polished  brass. 

Besides  the  similarities,  the  samples  exhibit  some  differences.  Sample  5 
had  almost  completely  converted  out  to  only  5%  orthohydrogen.  Sample  7  was 
at  x  =  50%  and  sample  8  at  x  =  70%.  This  may  be  the  reason  for  the  different 
onset  temperatures  of  the  second  steep  pressure  increase.  The  absolute 
pressures  are  a  function  of  the  hydrogen  filling,  the  amount  of  exchange  gas  and 
the  preparation  conditions. 

The  result  from  the  data,  gained  from  heating  the  hydrogen  sample 
isochorically,  is  the  existence  of  a  two-component  system.  One  contribution  is 
the  more  volatile  bulk-like  component  in  the  pore  center  which  becomes  gaseous 
and  leaves  the  sample  first,  leading  to  an  increase  in  pressure.  The  other 
contribution  originates  from  the  more  strongly  bound  component  on  the 
substrate  wall  which  stays  longer  in  the  zeolite  and  leads  to  an  estimate  for  the 
adsorption  energy.  This  is  in  impressive  agreement  with  the  indications  gained 
from  the  pulse  and  continuous  wave  work. 
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5.6.    Zeolite  5A  versus  13X 

Additional  studies  were  carried  out  using  zeolite  with  smaller  dimensions 
in  order  to  determine  the  relative  significance  of  the  wall  behavior.  Zeolite  5A 
was  chosen  because  of  its  smaller  dimensions  with  5  A  channels  and  9  A  pores. 
This  would  allow  for  comparison  with  the  previous  data  since  the  chemical 
composition  and  the  microscopic  structure  is  closely  related.  Remarkable  is  that 
virtually  all  hydrogen  molecules  have  direct  contact  with  the  substrate  walls. 

The  supercooling  transitions  observed  for  hydrogen  in  zeolite  5A  are 
much  sharper.  The  measured  (Figure  25a)  and  corrected  (Figure  25b)  echo 
amplitudes  are  almost  constant  until  about  13  K  and  drop  to  zero  at  10  K.  The 
irregularities  at  1 1  K  are  interpreted  as  fluctuations  of  the  occurring  liquid/solid 
phase  transition.  T2  drops  from  2  msec  above  1 1  K  to  values  <10//sec  which  is 
shorter  than  the  recovery  time  of  the  spectrometer  and  echoes  cannot  be 
observed  anymore.  The  cw  NMR  line  shapes  differ  in  a  controlled  fashion.  In 
contrast  to  zeolite  13X,  a  narrower  wing  contribution  is  observed  (Figure  26). 
The  width  at  half  height  has  changed  little  (20  kHz)  but  the  line  is  only  50  kHz 
wide  instead. of  130  kHz  below  4  K  (Figure  26  b).  This  confirms  the  view  that 
orientational  ordering  is  reduced  in  very  small  clusters.  The  line  shape  broadens 
from  40  kHz  at  5.95  K  to  110  kHz  at  1.75  K  upon  cooling  (Figure  26)  also 
indicating  a  higher  degree  of  ordering.  The  picture  arising  from  the  width  of  the 
resonance  line  is  that  the  cluster  size  decreases  with  the  reduction  in  available 
space. 

The  adsorption  potential  is  again  extracted  from  isochoric  pressure  versus 
temperature  data  and  is  found  to  be  smaller  than  for  13X,  being  only  100  K. 

This  complementary  data  on  another  type  of  zeolite  (5A)  confirms  the 
findings  made  for  zeolite  13  X  with  the  tendency  towards  more  extreme  results. 
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Figure  25.  Echo  amplitudes  for  hydrogen  in  zeolite  5A 
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Figure  26.  Temperature  dependence  of  cw  lineshape  in  zeolite  5A 
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Figure  26.  --  continued 


CHAPTER  6 
COMPUTER  SIMULATION 


From  the  data  described  in  chapter  5  the  system  has  been  deduced  to  be 
a  two-component  system.  It  is  therefore  desirable  to  quantify  the  contributions 
and  the  exact  order  parameters  involved  in  the  quasistatic  lineshape.  The 
glassiness  of  the  line  suggests  a  superposition  of  Pake  doublets  which  where 
discussed  in  chapter  3.  The  question  is  how  wide  the  order  parameter 
distribution  is  and  whether  the  observed  lineshape  is  reproducible  with  a 
combination  of  a  Pake  doublet  superposition  and  a  Gaussian  line.  The  Pake 
doublets  represent  the  partially  ordered  molecules  on  the  substrate  walls 
whereas  the  Gaussian  is  the  result  of  the  bulk-like  molecules  in  the  center  of  the 
pores. 

Starting  from  this  hypothesis  a  computer  program  was  written  to  simulate 
the  observed  lineshape  theoretically.  The  program  was  written  on  a  personal 
computer  with  processor  80386  at  20  MHz  with  a  math  coprocessor.  The 
computation  time  is  still  about  one  hour  per  spectrum.  The  language  used  was 
QUICKBASIC  in  the  framework  of  LabWindows,  a  powerful  data  analysis 
software  program.  The  task  is  split  into  three  parts.  Program  1,  incorporating 
most  of  the  work,  performs  the  generation  and  superposition  of  the  Pake 
doublets.  Program  2  generates  the  Gaussian  lineshape  while  program  3  adds 
the  two  components  for  a  plot  of  the  simulated  lineshape. 

Going  through  program  1  will  reveal  the  problems  and  solutions  to  the 
variable,  unbiased  and  symmetric  generation  of  Pake  doublets  and  their 
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summation.  Therefore,  a  printout  of  the  program,  as  it  was  run  for  the  final 
version  of  the  simulated  lineshape,  is  presented: 

PROGRAM  1 

PAKE  DOUBLET  SUPERPOSITION 

DIM  shared  1(3001) 
DIM  shared  X(3001) 
DIM  shared  v(3001) 

Dim  shared  funct.filename  as  string  *  30 
Dim  shared  freq.filename  as  string  *  30 
Dim  shared  graph.filename  as  string  *  30 
Dim  shared  f1  as  string  *  30 

c$="c:\LW\3dpake\" 

print  "This  program  assumes  data  is  written  to  ";c$ 
input  "Name  for  Function  file*;f1$ 

a1  %=fmt(funct.filename$,"%s<%s%s",c$,f  1  $) 

b1%=fmt(freq.filename$,"%s<%s%s%s",c$,f1$,"fr")  'name  frequency  file  with  "fr" 
c1%=fmt(graph.filename$,"%s<%s%s%s",c$,f1$,"g")  'name  graph  file  with  "g" 

oif%  =  openfile(funct.filename$,0,0,1 )  'save  array 
ovf%  =  openfile(freq.filename$,0,0,1 )  'save  array 

d  =  54.2  'solid  state  constant;  for  liquid  state  the  constant  is  57.68  kHz 
p  =  0 
s  =  1 

CALL  Grf  Reset  (12) 

z  =  1000  '  Half  number  of  sampling  points 

FOR  p  =  0  to  3001      1  p  is  the  array  counter  for  the  function  values 

l(p)-0 

X(p)  =  0 

v(p)-0 
NEXT  p 

'  Introduce  a  small  parameter  e  to  smoothing  the  divergences 
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e  =  1e-2 

'  Scan  through  the  order  parameter 

FOR  t%  =  1 0  to  75  stepl 

s  =  t%/100 
print  "s=";  s 

'  Compute  the  function  values  of  one  wing  of  the  Pake  doublet 

'  The  number  of  function  values  is  proportional  to  the  order  parameter 

'  such  that  summation  of  lines  is  unbiased 

p  =  3/2*z*s 

FOR  l&  =  50*z  to  -100*z  step  -100/s 
j=l&/100/z 

v(p+(z-z*s))  =  3*j*d*s 

X(p+(z-z*s))  =  (3*1/2)*(1  -2*v(p+(z-z*s))/(3*d*s))*(1/2)/(1 8*d*s*(1  -2*v(p+z-z*s)/(3*d*s)+e)) 

P  =  P- 1 
NEXT  l& 

1  Compute  the  other  wing  of  the  Pake  doublet 
p  =  3/2*z*s  +  1 

FOR  m&  =  -50*z  to  100*z  step  100/s 
k#=m&/100/z 
v(p+(z-z*s))  =  3*k*d*s 

X(p+(z-z*s))  =  (3A1/2)*(1  +2*v(p+(z-z*s))/(3*d*s))A(1/2)/(1 8*d*s*(1  +2*v(p+z-z*s)/(3*d*s)+e)) 
p  =  p+  1 
NEXT  m& 

'  Summation  of  the  two  wing  components 

'  The  function  values  are  zero  up  to  where  the  line  starts  according  to  s 

FORp  =  0to  z-1-z*s 

X(p)  =  0 

v(p)  =  v(p) 
NEXT  p 

1  Only  one  wing  contributes  to  the  first  part  of  the  complete  line 
FOR  p  =  z+1-z*s  to  z-z*s/2 
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X(p)  =  X(p) 
v(p)  =  v(p) 
NEXT  p 

'  In  the  center  both  components  contribute 

FOR  p  =  (z-z*s/2)  to  z 

X(p)  =  X(p)  +  X(p-(-1+z*s) 

v#(p)  =  v#(p) 
NEXT  p 

'  Notice  the  frequency  correction  for  the  line  beyond  the  center  point 

FORp  =  (z+1)to  (z+z*s/2) 

X(p)  =  X(p)  +  X(p+1+z*s) 

v#(p)  =  v#(p  +  1  +  z*s) 
NEXTp 

'  Only  the  other  component  contributes 

FOR  p  =  z  +  1  +  z*s/2  to  z  +  z*s 

X(p)  =  X(p  +  z*s  +  1) 

v(p)  =  v(p  +  z*s  +  1) 
NEXTp 

'  The  function  values  are  zero  again  beyond  the  line  according  to  s 

FOR  p  =  z  +  1  +  z*s  to  3*z 
X(p)  =  0 

V(p)  =v(p  +  z*s  +  1) 
NEXTp 

'  Integrate  the  line  in  order  to  normalize  by  the  area 

Integra  =  0 

FOR  p  =  0  to  3000 

Integra  =  Integra  +  (X(p+1)  +  X(p))  /  2  *  (v(p+1)  -  v(p)) 
NEXTp 

print  "INTEGRAL  of  LINE  with  s=";s;"is";lntegra 

'  Normalize  by  Area  of  Line 

'  and  weight  the  lines  by  their  order  parameter  s 
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FOR  p  =  0  to  3000 
X(p)  =  X(p)/lntegra 

If  s  <  0.25  then 
s  =  0.5-s 
endif 

X(p)  =  X(p)/s  '  Stronger  contribution  for  pair  clusters 
NEXTp 

FOR  p  =  0  to  2000 
l(p)  =  X(p)  +  l(p) 
NEXTp 

'  End  the  order  parameter  loop 
NEXTt% 

'  Rle  the  function  and  frequency  values 

nif%  =  fmtfile(oif%,"%2001f[w10]\  l() ) 
nvf%  =  fmtfile(ovf%,"%2001f[w10]",  vO  ) 

cif%  =  closefile(oif%) 
cvf%  =  closefile(ovf%) 
'  Set  up  the  Graphics 

CALL  SetAxName  (0,  "FREQUENCY  (kHz)") 

CALL  SetTxAlign  (1,2) 

CALL  SetAxName  (1,  "PAKE  DOUBLETS") 

CALL  SetTrtle  ("3D  PAKE  DOUBLET  SUPERPOSITION") 

CALL  SetAxGridVis(-1,  0) 

CALL  SetGrdFrame  (1) 

CALL  SetCurv2D  (0) 

CALL  SetPointStyle  (10) 

CALL  SetXDataType  (4) 

CALL  SetYDataType  (4) 

CALL  GrfCurv2D  (v(),  l(),  2001) 

call  saveGraphfile(graph.filename$,  0)      '  save  ail  ports 
END  OF  PROGRAM  1 
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In  order  to  obtain  reasonable  accuracy  and  smoothness  of  the  line  the 
arrays  are  chosen  maximally  for  a  4  megabyte  RAM.  The  function  and 
frequency  values  were  computed  and  stored  in  files  such  that  immediate  access 
without  recalculation  was  possible.  This  was  useful  considering  the  long 
calculation  times  involved. 

Starting  from  the  Pake  doublet  equation 


where  v0  =  -^— ,  it  is  known  that  the  complete  Pake  doublet  lineshape  is  a  sum 

of  the  positive  and  negative  branches.  The  plus-sign  is  valid  for  -  v0  <  v+  <  +2  v0 
and  the  minus-sign  for  -2v0  <  v_  <  v0.  The  first  difficulty  was  to  guarantee  that 
the  frequencies  for  lineshapes  of  different  order  parameters  a,  here  called  s, 
coincide  and  have  the  same  sampling  density.  This  was  necessary  for  the 
unbiased  summation  of  different  lines.  The  problem  was  solved  by  s  dependent 
function  sampling  centered  around  the  middle  of  the  total  sampling  interval,  here 
2000  (two  times  z)  data  points.  The  order  parameter  is  scanned  through  from  a 
minimum  to  a  maximum  value  in  steps  of  0.01  to  get  reasonable  smoothness  of 
the  line.  One  complete  Pake  doublet  itself  consists  of  two  separate  wings. 
These  have  to  be  generated  individually.  Each  contribution  is  stored  in  1500 
times  s  array  points  representing  the  equal  frequency  sampling  density.  The 
reverse  order  for  the  first  contribution  is  necessary  so  that  the  line  is  symmetric 
for  all  order  parameters  s  around  array  1000.  The  main  complication  with  the 
Pake  doublets  are  the  divergences  at  the  positive  and  negative  frequency 
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Vn  = 


3da 


,  where  the  cosine  of  90  degrees  becomes  zero,  corresponding  to  the 


infinite  contribution  from  the  two-dimensional  plane.  These  divergences  must  be 
regularized  in  a  fashion  that  is  insensitive  to  an  infinitesimal  regularization 
parameter  e.  In  principal  there  are  four  different  ways  to  attempt  this  task.  The 
correct  method  is  the  one  where  the  result  is  independent  of  e  in  the  limit  of 
vanishing  e.  The  methods  can  be  separated  into  introducing  the  parameter  for 
the  frequency  itself  or  in  the  entire  denominator  and  whether  it  is  real  or 
complex.  For  a  real  regularization  parameter  for  the  frequency  the  Pake  doublet 
equation  takes  on  the  form 


a±j  v\  =      —     1       oc  pi  at  divergence 


and  for  the  whole  denominator 


0*2  ( v)  -  —  ,       —   -  at  divergence 


For  an  imaginary  regularization  parameter  only  the  real  part  of  the  expression 
must  be  considered 


3*3(10  =  Re 
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where  v  =  1 — -  which  is  zero  at  divergence.  The  final  possibility  is 


gt<(v)  =  Re- 


\Qdo 


V3 


The  interesting  property  of  all  expressions  is  the  behavior  at  divergence. 
Expressions  1  through  3  have  an  obvious  e-dependence  at  the  divergence 
frequency  which  renders  them  inept  for  the  regularization  procedure  because 
they  are  highly  sensitive  to  the  externally  introduced  regularization  parameter. 
The  last  expression,  however,  has  the  special  property  of  being  zero  at  the 
divergence  frequencies,  and  is  therefore  independent  of  e.  The  important 
difference  is  that  the  frequency  in  the  denominator  is  small  compared  to  the 
square  of  e,  which  is  small  but  finite,  and  the  numerator  itself  goes  to  zero 
independent  of  the  value  of  e  as  the  divergence  is  approached.  The  line 
converges  quickly  to  its  final  form  hardly  changing  for  e  values  smaller  than  10~2. 
This  value  was  chosen  to  incorporate  smoothing  of  the  divergences  which  in 
practice  is  a  consequence  of  the  ZPM,  the  temperature  effect  and  the 
intermolecular  interactions. 

The  summation  of  the  two  components  requires  a  careful  addition  of  the 
regions  where  both  are  zero,  only  one  contributes  and  where  both  parts  are 
nonzero.  This  care  is  necessary  for  shifting  the  second  wing  component  down  to 
the  same  frequencies  as  for  the  first  wing.  The  lines  are  then  integrated  to 
normalize  them  by  dividing  each  line  by  its  area. 

Another  issue  is  the  final  weighting  of  the  lines  with  different  order 
parameters.  It  is  conceivable  in  zeroeth  order  to  weight  them  evenly  and  sum 
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them  directly.  In  first  order  the  lines  can  be  multiplied  or  divided  by  the  order 
parameter.  The  most  physical  process  is  to  attribute  more  importance  to  small 
order  parameters.  Small  s  implies  small  clusters  which  are  very  likely  for  the 
small  pores  of  zeolite.  A  lower  cutoff  at  the  value  for  pairs  of  hydrogen 
molecules  can  be  introduced  which  is  equivalent  to  s  =  0.25.  The  Gaussian  line 
for  disordered  material  covers  the  remaining  narrow  contribution.  The  calculated 
lines  are  then  added  to  yield  the  superposition  of  Pake  doublets.  The  final  part 
of  the  program  is  concerned  with  plotting  the  lineshape.  The  result  is  shown  in 
figure  27  for  constant  weighting  of  a  between  0.1  and  0.75. 

Program  2  is  a  simple  generation  of  a  Gaussian  line  with  varying  area  and 
line  width.  The  function  values  are  also  stored  in  files  and  plotted  (Figure  28). 
This  line  contributes  about  20%  to  the  total  lineshape  area  and  has  a  Gaussian 
width  parameter  of  5  kHz. 

Program  3  takes  the  calculated  values  from  programs  1  and  2,  adds  them 
up  and  plots  them.  The  lines  are  somewhat  fuzzy  in  their  contour  because  the 
number  of  lines  superimposed  is  still  not  large  enough  to  yield  a  completely 
smooth  curve.  The  line  is,  however,  clearly  visible  and  interpretable. 

The  results  are  displayed  in  figure  29a  for  constant  weighting  of  s 
between  0.1  and  0.75.  The  shape  is  not  quite  satisfactory  for  the  wings  where 
the  experimental  line  is  strongly  convex.  An  attempt  was  made  to  remedy  this 
shortcoming.  The  large  values  of  s  from  0.75  to  0.25  were  weighted  inversely. 
This  introduces  the  discussed  cutoff  at  the  value  for  pair  clusters.  The  order 
parameters  smaller  than  0.25  were  weighted  inversely  after  correcting  them  to 
0.5  -  s.  This  reduces  their  importance  as  well.  The  reasons  are  that  inverse 
weighting  of  small  s  introduces  an  unphysical  dominance  of  clusters  that  are  on 
average  smaller  than  pairs  and  that  the  Gaussian  represents  the  smallest  order 
parameters.  The  result  is  shown  in  figure  29b. 
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Figure  27.  Pake  doublet  superposition 
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Figure  29.  Synthesized  cw  lineshapes 
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The  result  of  the  simulation  is  that  the  order  parameters  cr,  i.e.  s,  follow  a 
distribution  around  a  mean  of  about  0.45  with  a  width  of  0.3.  The  important 
conclusion  in  this  analysis  is  that  the  outer  wings  of  the  Pake  doublet 
superposition  are  below  the  noise  level  of  the  experiment  and  therefore  not 
observable  experimentally.  This  explains  why  the  line  is  by  far  not  as  wide  as  a 
fully  ordered  hydrogen  line  with  a  line  width  of  about  330  kHz.  On  the  other 
hand  it  shows  why  the  line  is  broad  with  some  components  close  to  perfect 
ordering  (s  =  1).  The  interpretation  of  this  cw  NMR  result  is  that  at  least  a  certain 
percentage  of  ortho  molecules  is  partially  ordered  in  small  clusters.  The 
roughness  and  inhomogeneity  of  the  zeolite  surface  is  expected  to  lead  to  such 
a  distribution  of  a  for  partially  ordered  states. 

The  simulated  line  is  compared  with  the  experimental  line  and  resembles 
it  closely  in  width  and  appearance. 


CHAPTER  7 
CONCLUSION 


In  chapter  5  the  measured  data  were  presented  and  explanations  given 
for  the  origin  of  the  behavior  of  hydrogen  in  zeolite.  In  contrast,  the  present 
chapter  will  reorganize  the  data  in  terms  of  physical  effects.  The  data  are  used 
to  develop  ideas  about  the  system  and  its  behavior  as  a  whole.  The  conclusion 
chapter  is  in  that  sense  paralleling  the  introductory  chapter  in  that  I  review  the 
open  questions  about  confinement  for  quantum  solids.  This  chapter  will  try  to 
give  answers  to  the  raised  questions  as  far  as  they  could  be  resolved  during  the 
course  of  the  measurements. 

The  studies  mainly  focussed  on  Nuclear  Magnetic  Resonance  data.  The 
hydrogen-zeolite  system  was  probed  in  coherent  pulse  and  continuous  wave 
(cw)  modes  at  268  MHz.  These  two  methods  complemented  each  other  with 
their  information  about  dynamical  properties,  the  quasistatic  order  parameter  a 
and  the  lineshape.  Additional  information  was  gained  from  pressure  data  as  a 
function  of  temperature. 

1)  Substrate  Potential  and  Two-Component  System 

The  substrate  potential  that  acts  on  the  adsorbed  molecules  constitutes  a 
major  difference  to  a  three  dimensional  system.  Its  impact  is  very  dominant 
concerning  the  behavior  of  the  hydrogen  molecules.  An  important  implication  of 
the  data  is  the  existence  of  a  two-component  system  which  has  been  attributed 
to  the  surface  potential.  This  manifests  itself  in  all  experiments  performed.  The 
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relaxation  process  exhibits  two  relaxation  rates  for  the  pulse  data.  The 
temperature  dependence  of  the  lineshape  shows  the  disappearence  of  the 
center  component  upon  warming  while  the  broad  contribution  remains.  Also,  as 
the  sample  ages  the  narrow  center  component  decays  faster  than  the  wings.  A 
final  indication  is  the  isochoric  pressure  versus  temperature  data.  The  fashion  in 
which  the  hydrogen  leaves  the  sample  clearly  demonstrates  two  independent 
contributions.  The  adsorption  energy  of  hydrogen  on  the  zeolite  surface  was 
determined  to  be  270 ±50  K  from  the  slope  of  the  second  sieep  increase  in 
pressure.  The  explanation  for  the  two-component  system  is  the  substrate  wall 
potential.  It  binds  the  molecules  in  contact  with  the  wall  much  more  strongly 
than  molecules  in  the  pore  center. 

2)  Zero-Point  Motion 

The  impact  of  the  zero-point  motion  was  not  measured  directly.  From  the 
behavior  of  other  parameters  it  is,  however,  possible  to  draw  conclusions  of  its 
importance.  From  the  previously  mentioned  property  of  a  two-component 
system  it  follows  that  the  ZPM  is  not  large  enough  to  compensate  for  the 
attraction  of  the  surface  potential.  If  that  was  the  case  the  molecules  on  the 
walls  would  have  enough  energy  to  be  similar  in  their  behavior  to  the  molecules 
in  the  bulk-like  environment.  This  implies  that  the  ZPM  is  relatively  weak 
compared  to  the  substrate  potential.  Another  indication  for  the  importance  of  the 
ZPM  is  the  width  of  the  transition.  The  fluctuations  related  to  the  ZPM  widen  the 
transition.  In  the  liquid/solid  transition  of  this  experiment  the  transition  width  is 
large.  This  width  is,  however,  interpreted  as  a  result  of  the  system 
heterogeneity,  instead  of  an  excessive  ZPM  contribution.  The  environments 
experienced  by  the  hydrogen  molecules  are  various  due  to  different  degrees  of 
contact  with  the  substrate  walls  as  a  consequence  of  the  nonspherical  pore 


163 


shape.  The  critical  region  is  therefore  spread  out  over  1  to  3  degrees.  The 
reason  is  that  different  sample  contributions  experience  the  critical  region  at 
different  temperatures.  The  displayed  transition  is  a  superposition  of  individual 
transitions.  This  is  in  distinct  contrast  to  the  three-dimensional  system  where  the 
critical  region  is  very  narrow. 

3)  Orientational  Ordering 

The  expectation  that  quadrupolar  ordering  takes  place  at  higher 
temperatures  than  for  bulk  originates  from  the  assistance  of  the  ordering  process 
by  the  surface  potential.  It  is  expected  that  the  surface  potential  reduces 
fluctuations  and  locks  the  molecules  in  place  before  this  happens  for  an 
unconfined  system.  The  peak  in  T2  could  be  interpreted  as  the  onset  of 
orientational  ordering.  This  was,  however,  rejected  for  several  reasons. 

The  line  width  was  unaffected  by  a  change  in  temperature  around  the 
transition  temperature.  It  was  still  broad  at  temperatures  as  high  as  15.8  K 
which  implies  the  existence  of  nonzero  order  parameters  for  these  temperatures. 
The  remarkable  line  width  of  over  100  kHz  at  16  K  must  be  attributed  to  partial 
ordering. 

Another  piece  of  information  is  the  line  broadening  from  4.2  to  0.8  K  which 
evidences  additional  ordering  as  the  bulk  ordering  temperature  is  approached.  It 
is  important  to  note,  however,  that  even  for  the  lowest  temperatures  investigated 
no  complete  ordering  is  achieved  in  contrast  to  the  bulk  samples  with  <j  =  1. 

Measurements  on  hydrogen  deuteride,  which  has  no  orientational 
degrees  of  freedom,  resulted  in  similar  findings  as  for  hydrogen  which  leads  to 
the  conclusion  that  the  peaks  at  10  K  are  not  a  result  of  orientational  ordering. 
This  experimental  evidence  is  a  very  convincing  argument. 
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Hydrogen  in  zeolite  5A  with  pores  of  only  8  A  was  also  investigated. 
Similar  but  more  extreme  results  were  obtained  in  agreement  with  the  previous 
data.  The  obtained  lineshapes  are  50  kHz  wide  at  4  K.  The  small  width  is 
related  to  a  maximum  order  parameter  of  0.3,  indicating  small  clusters.  The 
picture  that  arises  from  the  line  widths  measured  for  hydrogen  in  bulk,  zeolite 
13X  and  zeolite  5A  is  that  the  cluster  size  decreases  with  the  smaller  geometry. 
This  is  an  intuitively  obvious  statement. 

4)  Melting  Transition  and  Supercooling 

When  a  system  shows  unusual  orientational  properties  the  logical 
extension  is  to  investigate  the  translational  transition.  The  pulse  data  revealed  a 
strong  peak  in  the  transverse  relaxation  time  T2  at  10  K.  A  similar  peaking  was 
also  observed  for  the  nuclear  spin  echoes.  This  behavior  was  interpreted  in 
terms  of  supercooling  and  implies  that  hydrogen  in  zeolite  solidifies  about  5 
degrees  below  its  bulk  melting  point.  Two  features  are  especially  remarkable 
about  the  supercooling  transition.  Firstly,  the  T2  was  found  to  be  nonhysteretic 
which  emphasizes  a  purely  energetically  driven  process.  This  result  is  strongly 
backed  by  the  T2  peaks  of  HD  which  were  almost  identical  upon  cooling  and 
warming.  In  addition,  their  variation  is  such  that  the  temperature  difference  is 
reversed,  compared  to  what  would  be  consistent  with  hysteresis.  Secondly,  the 
supercooling  transition  temperatures  exhibit  a  linear  parahydrogen  concentration 
dependence  of  -O.10±0.01  K/%.  This  implies  that  parahydrogen  can  be 
supercooled  more  easily  than  orthohydrogen. 

Supercooling  can  be  explained  using  two  alternative  theories:  nucleation 
theory  or  the  raised  free  energy  for  the  solid  state  because  of  grain  boundaries. 
In  the  grain  boundary  model  a  thin  layer  wets  and  solidifies  on  the  zeolite 
surface.  The  temperature  is  depressed  because  grain  boundaries  in  the  solid 
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state  cost  additional  energy  and  the  solid  state  becomes  favorable  only  at  lower 
temperatures.  In  the  framework  of  nucleation  theory  the  surface  is  not  wetted. 
Instead,  the  solid  develops  from  a  fluctuation  in  the  center  of  the  pore  which 
grows  to  its  equilibrium  size,  the  pore  diameter  of  zeolite.  The  two-component 
system  found  in  the  present  case  suggests  a  combination  of  these  models.  A 
monolayer  of  hydrogen  is  in  direct  contact  with  the  surface,  wetting  it.  Its  state  is 
different  from  the  bulk  case.  This  can  be  seen  from  the  high  temperature 
lineshape  which  is  extremely  wide.  Although  probably  liquid,  hydrogen  is  still 
strongly  confined  to  the  substrate  wall  by  the  attraction  potential.  The  remaining 
material  in  the  center  of  the  pore  stays  liquid  until  it  solidifies  according  to 
nucleation  theory.  The  immediate  solidification  on  the  surface  layer  is  inhibited 
by  the  grain  boundary  energy. 

An  alternative  to  supercooling  is  the  dimensionality  effect.  From  a 
molecular  dynamics  study  by  Brougthon  et  al.,  it  was  shown  that  a  solid  melts  in 
layers.  The  outermost  layer  melts  at  72  %  of  the  bulk  melting  temperature.  This 
leads  to  a  melting  onset  of  just  below  10  K.  This  is  in  close  agreement  with  the 
experiment  and  would  imply  that  the  melting  in  the  pore  is  taking  place  in  a  two- 
dimensional  fashion  from  the  layer  in  contact  with  the  wall.  The  bulk  contribution 
is  small  and  therefore  negligible.  Instead  of  supercooling  it  would  represent  a 
two-dimensional  melting  point.  This  seems  unlikely,  however,  because  of  the 
surface  interaction  which  does  not  allow  for  viewing  the  molecules  on  the  surface 
as  a  free,  two-dimensional  layer. 

The  investigation  of  the  T2  transit'on  in  terms  of  supercooling  was  a  major 
part  of  the  performed  experiment.  Several  indications,  however,  confirm  this 
assumption.  One  indirect  argument  is  that  an  orientational  ordering  transition 
can  be  excluded,  as  already  discussed  in  the  previous  section  about 
orientational  ordering.  The  direct  argument  is  an  experiment  performed  with 
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hydrogen  deuteride  (HD)  instead  of  hydrogen.  HD  is  in  a  pararotational 
groundstate  and  does  not  have  orientational  degrees  of  freedom,  in  this  sense 
imitating  parahydrogen.  HD  is,  however,  detectable  with  NMR.  Echo 
amplitudes,  T2  and  7",  data  are  in  beautiful  agreement  with  the  hydrogen  data 
and  therefore  completely  rule  out  an  interpretation  of  the  10  K  peak  in  terms  of 
orientational  ordering.  The  lack  of  irregularities  at  13.8  K,  which  is  the  bulk 
melting  point,  also  points  toward  suppressed  melting  for  the  confined  system. 
For  one  echo  amplitude  transition  a  discontinuity  was  observed  which  is  a  further 
confirmation  because  it  represents  the  bulk  effect  in  addition  to  the  suppressed 
transition  at  10  K.  The  ostentatious  coincidence  of  the  T2  peak  with  the  drop  in 
corrected  echo  amplitude  is  a  final  confirmation  of  the  liquid/solid  transition  at 
this  temperature.  The  measured  transition  at  8.5  K  is  the  largest  value  of 
supercooling  in  the  literature  achieved  for  hydrogen. 

5)  Orientational  Glass 

The  delayed  and  inhibited  solidification  has  another  implication.  When 
the  solid  finally  prevails,  it  incorporates  a  high  density  of  grain  boundaries  and 
strain.  Under  these  conditions  it  is  easily  conceivable  that  the  order  parameter  is 
widely  distributed  as  a  consequence  of  the  rough  substrate  surface.  This  is 
exactly  what  the  continuous  wave  data  show:  a  glass-type  lineshape  with  a 
width  of  130  kHz.  This  shape  represents  a  wide  line  for  disordered  hydrogen 
and  a  narrow  line  for  completely  ordered  hydrogen.  It  is  a  combination  of  broad 
wings  with  a  narrow  component  in  the  center.  This  becomes  obvious  upon 
warming  from  4.2  K  when  the  narrow  line  disappears  but  the  broad  contribution 
remains  without  a  sharp  transition  at  10  K.  The  understanding  of  the  lineshape 
found  its  manifestation  in  a  successful  simulation  which  is  a  result  of  the  sum  of 
a  superposition  of  Pake  doublets  -  with  a  distribution  of  the  order  parameter 
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between  0.1  and  0.75  -  and  a  narrow  Gaussian  line.  A  distribution  of  the  order 
parameter  is  a  classical  indication  for  an  orientational  glass.  Its  frustration  is  the 
incomplete  orientational  ordering  which  represents  the  failure  of  adopting  the 
minimum  energy. 

6)  Ortho-Para  Conversion 

It  is  necessary  to  know  the  ortho-para  hydrogen  conversion  rate  in  order 
to  characterize  the  constantly  changing  ortho-para  ratio  of  the  sample.  Any  data 
must  be  labelled  by  its  concentration  because  of  its  strong  dependence  on  the 
ortho-para  contents.  The  confined  environment  is  expected  to  cause  changes  in 
the  rate.  The  conversion  parameters  were  extracted  from  the  time  dependence 
of  the  lineshape  area.  The  conversion  presents  itself  as  a  bimolecular  process 
with  two  different  constants  for  early  and  late  times.  For  the  first  500  hours  the 
conversion  is  characterized  by  k,  =0.425 ±0.006% /h  which  is  only  a  fourth  of 
the  bulk  value.  This  behavior  is  explained  by  the  reduced  number  of  nearest 
neighbors.  The  molecules  on  the  wall  have  less  than  half  the  number  of  nearest 
neighbors  because  the  curved  wall  represents  more  than  half  of  a  molecules' 
space  on  the  wall.  The  surface  potential  reduces  the  fluctuating  magnetic  stray 
field  originating  from  a  molecule  by  restricting  its  motion.  Another  effect  on  the 
conversion  rate  is  the  reduced  density  of  phonon  states  with  the  dimension  due 
to  the  decrease  from  3  to  2  phonon  modes.  For  later  times  the  conversion 
accelerates  to  k2  =  2.21  ±  0.075% /h  and  is  explained  in  terms  of  orthohydrogen 
clustering  effects  which  make  conversion  more  efficient  due  to  more  neighbors 
and  the  decreased  separation  between  molecules.  The  cw  lineshape  changes 
over  time  in  the  sense  that  the  narrow  component  decays  •  faster.  This  is 
consistent  with  the  explanation  of  nearest  neighbors  and  a  bulk-like  environment 
for  the  molecules  in  the  center  of  the  pores. 
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7)  Super-fluid  Hydrogen 

The  ambitious  perspective  of  the  project  was  to  search  for  superfluid 
molecular  hydrogen.  The  idea  was  to  supercool  hydrogen  below  the  region  of 
superfluidity  which  is  predicted  to  be  below  6  K.  The  additional  supercooling 
from  the  small  pore  diameter  of  zeolite  gave  rise  to  the  hope  of  achieving  such 
extreme  values.  The  results  express  some  encouragement.  The  supercooling 
temperature  is  fit  by  a  linear  function  of  the  parahydrogen  concentration  with  a 
dependence  of  -0.10 ±0.01  K/%.  This  implies  more  efficient  supercooling  for 
para-  than  orthohydrogen.  Although  the  lowest  measured  liquid/solid  transition 
at  8.5  K  is  still  several  degrees  short  of  the  predicted  superfluid  temperature  this 
method  has  some  potential.  The  reason  is  that  if  the  linear  interpolation  of  the 
supercooling  temperature  to  parahydrogen  concentration  relationship  held  up  to 
100  %  parahydrogen,  the  predicted  liquid/solid  transition  temperature  was  at 
3.3 ±0.8  K.  This  could  be  within  the  reach  of  superfluidity.  The  NMR  method  is 
unfortunately  inept  to  make  a  statement  about  such  concentrations  because 
parahydrogen  by  itself  is  not  detectable. 

The  surface  tension  aLS  was  mentioned  to  be  a  crucial  parameter 
for  successful  supercooling  in  nucleation  theory.  This  value  was  found  to  be 
greater  than  1  erg/cm  by  Tell63  et  al.  for  hydrogen  in  vycor  with  a  radius  of  27  A. 
This  aLS  would  theoretically  be  large  enough  for  successful  supercooling  down 
to  the  superfluid  transition.  This  value  for  aLS  must  be  corrected  downward, 
however,  by  the  present  zeolite  experiment  with  the  result  of  one  fourth  the 
radius  but  only  twice  the  temperature  suppression.  This  renders  supercooling 
more  difficult  to  achieve  than  expected  by  the  earlier  measurements. 

The  question  is,  however,  if  it  makes  sense  at  all  to  expect  superfluidity 
for  a  system  that  only  comprises  about  20  to  30  molecules  per  pore.    It  is 
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doubtful  whether  collective  behavior  such  as  Bose  condensation  is  possible  for 
such  a  system. 

The  performed  experiment  on  hydrogen  confined  to  the  pores  of  zeolite 
explains  a  variety  of  questions  and  makes  a  contribution  to  the  general 
understanding  of  properties  that  are  related  to  a  system  characterized  by 
confinement.  Further  investigations  are  necessary  to  shed  light  on  the 
mechanism  of  the  suppressed  solidification,  the  understanding  tor  properties  of 
small  numbers  of  particles,  where  statistical  mechanics  do  not  apply  in  the 
traditional  sense,  and  the  existence  of  superfluid  hydrogen. 


BIBLIOGRAPHY 


1  Slichter,  C.  P.,  Principles  of  Magnetic  Resonance  3rd  edition,  Springer-Verlag, 
Heidelberg,  1990 

2  Silvera,  I.  F.,  Goldman,  V.V.,  J.  Chem.  Phys.  69,  4209  (1978) 

3  Rail,  M.,  Master's  Thesis,  ,  unpublished,  University  of  Florida  (1988) 

4  Brewer,  D.  F.,  Rajendra,  J.,  Sharm,  N.,  Thomson,  A.  L,  Xin,  J.,  Physica  B 
165&166,  577  (1990) 

5  Vidali,  G.,  Cole,  M.  W.,  Schwartz,  C,  Surf.Sci.  87,  L273  (1979) 

6  Novaco,  A.  D.,  Phys.  Rev.  Lett.  60,  2058  (1988) 

7  Novaco,  A.  D.,  Wroblewski,  J.  P.,  Phys.  Rev.  B  39,  1 1364  (1989) 

8  Novaco,  A.  D.,  Phys.  Rev.  B  19,  6493  (1979) 

9  Nielsen,  M.,  Ellensen,  W.,  in  Procedings  of  the  14th  International  Conference  on 
Low  Temperature  Physics,  Otaniemi,  Finland,  1975,  edited  by  M.  Krusius  and 
M.  Vuorio,  North-Holland,  Amsterdam,  Vol.  4,  437  (1975) 

10  Berthold,  J.  E.,  Bishop,  D.  J.,  Reppy,  J.  D.,  Phys.  Rev.  Lett.  39,  348  (1977) 

11  Silvera,  I.  F.,  Rev.  Mod.  Phys.  52,  393  (1980) 

12  Nielsen,  U.,  Halstead,  D.,  Holloway,  S.,  Norskov,  J.  K.,  J.  Chem.  Phys  93, 
2879(1990) 

13  Huber,  T.  E.,  Huber,  C.  A.,  Physica  B,  LT19  (1990) 

14  Brewer,  D.  F.,  Evenson,  A.,  Thomson,  A.  L,  J.  Low  Temp.  Phys.  3,  603  (1970) 

15  Ambegaokar,  V.,  Halperin,  B.  I.,  Nelson,  D.  R.,  Siggia,  E.  D.,  Phys  Rev  B  21, 
1806  (1980) 

16  Chester,  M.,  Yang,  L.  C,  Phys.  Rev.  Lett.  31,  1377  (1973) 

17  Huber,  T.  E.,  Huber,  C.  A.,  Phys.  Rev.  Lett.  59,  1 120  (1987) 

18  Bishop,  D.  J.,  Reppy,  J.  D.,  preprint  (1980) 

19  Matsuda,  H.,  Van  Den  Meijdenberg,  J.  N.,  Physica  26,  939  (1960) 


170 


171 

2°  Berker,  A.  N.,  preprint  (1980) 

21  Cohen  de  Lara,  E.,  Kahn,  R.,  J.  Physique  42,  66  (1981) 

22  Huber,  T.  E.,  Huber,  C.  A.,  Phys.  Rev.  B  40,  12527  (1989) 

23  Dozier,  W.  D.,  Drake,  J.  M.,  Klafter,  J.,  Phys.  Rev.  B  56,  197  (1986) 

24  Kaerger,  J.,  Lenzner,  J.,  Pfeifer,  H.,  Schwabe,  H.,  Heyer,  W.,  Janowski,  F., 
Wolf,  F.,  Zdanov,  S.  P.,  J.  Am.  Ceram.  Soc.  66,  69  (1983) 

25  Smith,  J.  M.,  Chemical  Engineering  Kinetics,  McGraw-Hill,  New  York,  1981 

26  Lindemann,  F.  A., Phys.  Zeit.  11,  609  (1910) 

27  Born.M.,  Chem.  Phys.  7,  591  (1939) 

28  Kuhlmann-Wilsdorf,  D.,Phys.  Rev.  A  140,  1599  (1964) 

29  Cotterill,  R.  M.  J.,  Jensen,  E.  J.,  Kristensen,  W.  D.,  Paetsch,  R.,  Esbjorn,  P.  O., 
J.  Physique  36  C2,  35  (1975) 

30  Broughton  J.  Q.,  Woodcock  L  V.,  J.  Phys.  C  11,  2743  (1978) 

31  Tallon,  J.  L,  Robinson,  W.  H.,  Smedley,  S.  I.,  Philos.Mag.  36,  741  (1979) 

32  Mermin,  N.  D.,  Phys.  Rev.  176,  176  (1968) 

33  Peierls,  R.  E.,  Helv.  Phys.  Acta  7,  81  (1923) 

34  Kosterlitz,  J.  M.,  Thouless,  D.  J.,  J.  Phys  C  5,  L124  (1972) 

35  Kosterlitz,  J.  M.,  Thouless,  D.  J.,  J.  Phys.  C  6,118  (1973) 

36  Nabarro,  F.  R.  N.,  Theory  of  Crystal  Dislocation,  Oxford,  Clarendon,  1967 

37  Halperin,  B.  I.,  Nelson,  D.  R.,  Phys.  Rev.  Lett.  41,  121  (1978) 

38  Tallon,  J.  L,  Phys.  Rev.  B  22,  453  (1980) 

39  Tallon,  J.  L,  Phys.  Lett.  76A,  139  (1980) 

40  Lasocka,  M.,Phys.  Lett.  51  A,  137  (1974) 

41  Chui,  S.-T.,  Phys.  Rev.  Lett.  48,  933  (1982) 


172 

42  Landau,  L  D.,  Lifshitz,  E.  M.,  Statistical  Physics,  Addison-Wesley,  Reading, 
Massachusetts,  1969 

*5  Ginzburg,  V.  L,  Sobyanin,  A.  A.,  Sov.  Phys.  JETP  Lett.  38,  329  (1974) 

44  Maris,  H.  J.,  Seidel,  G.  M.,  Huber,  T.  E.,  J.  Low  Temp.  51,  471  (1983) 

45  Brewer,  D.  F.,  Rajendra,  J.,  Sharma,  N.,  Thomson,  A.  L,  Physica  B  165&166 
569  (1990) 

46  Widom,  A.,  Phys.  Rev.  176,  254  (1968) 

4?  Hohenberg,  P.  C,  Phys.  Rev.  158,  383  (1967) 

48  Shevchenko,  S.  I.,  Sov.  J.  Low  Temp.  Phys.  11,  363  (1985) 

49  Alekseeva,  L  A.,  Krupski,  1. 1.,  Sov.  J.  Low  Temp.  Phys.  10,  170  (1984) 
so  Geilikman,  B.  T.,  Sov.  Phys.  Solid  State  15,  2194  (1974) 

51  Bogolyubov,  N.  N.,  Sov.  Phys.  JETP  7,  289  (1958) 

52  Smith,  E.  N.,  Brewer,  D.  F.,  Liezhao,  C,  Reppy,  J.  D.,  Physica  107B,  585 
(1981) 

53  Deguchi,  H.,  Moriyasu,  Y.,  Amaya,  K.,  Kitaoka,  Y.,  Kohori,  Y.,  Jpn.  J.  Appl. 
Phys.  Suppl.  26,  311  (1987) 

^Turnbull,  D.,  Fischer,  J.  F.,  J.  Chem.  Phys.  17,  71  (1949) 

55  Davies,  J.  T.,  Rideal,  E.  K.,  Interfacial  Phenomena,  Academic  Press,  New 
York  and  London,  1963 

56  Ahmed,  K.(  Khwaja,  F.A.,  Rahman,  T.-U.,  Int.  J.  Mod.  Phys.  3,1661  (1989) 

57  Seidel,  G.  M.,  Maris,  H.  J.,  Williams,  F.  I.  B.  Cardon,  G.,  Phys.  Rev.  Lett.  56, 
2380  (1986) 

58  Maris,  H.  J.,  Seidel,  G.,  M.,  Williams,  F.  I.  B.,  Phys.  Rev.  B  36,  6799  (1987) 
Lie-zhao,  C,  Brewer,  D.  F.,  Girit,  C,  Phys.  Rev.  B  33,  106  (1986) 

60  Shimoda,  M.,  Mizusaki,  T.,  Suzuki,  T.,  Hirai,  A.,  Phys.  Lett.  102A,  426  (1984) 
si  Adams,  E.D.,Tang,Y.  H.,  Uhlig,  K.,  Haas,  G.  E.,  J.  Low  Temp.  66,  85  (1987) 


173 

62  Bretz,  M.,  Thomson,  A.  L,  Phys.  Rev.  B  24,  467  (1981) 
«  Tell,  J.  L,  Maris,  H.  J.,  Phys  Rev.  B  28,  5122  (1983) 

64  Dash,  J.  G.,  Phys.  Rev.  B  25,  508  (1982) 

65  Balibar,  S.,  Edwards,  D.  O.,  Laroche,  C,  Phys.  Rev.  Lett.  42,  782  (1979) 

66  Landau,  J.,  Saam,  W.  F.,  Phys.  Rev.  Lett.  38,  23  (1973) 

67  Landau,  J.,  Lipson,  S.  G.,  Maeaettaenen,  L.  M.,  Balfour,  L.  S.,  Edwards,  D.  O., 
Phys.  Rev.  Lett.  45,  31  (1980) 

68  Motizuki,  K.,  Nagamiya,  T.,  J.  Phys.  Soc.  Jpn.  11,  93  (1956). 

69  Berlinsky,  A.  J.,  Hardy,  W.  N.,Phys.  Rev.  B  8,  5013  (1973) 
7°  Berlinsky,  A.  J.,  Phys.  Rev.  B  12, 1482  (1975). 

71  Rose,  M.  E.,  Elementary  Theory  of  Angular  Momentum,  John  Wiley  and  Sons, 
New  York,  1957 

72  Roberts,  R.  J.,  Daunt,  J.  G.,  J.  Low  Temp.  Phys.  6,  97  (1972) 

73  Nielsen,  M.,  Phys.  Rev.  B  7,  1626  (1963) 

74  Yucel,  S.,  Phys.  Rev.  B  39,  3104  (1989) 

75Cremer,  E.,  Polanyi,  M.,  Trans.  Farad.  Soc.  28,  435  (1932);  Z.  Phys.  Chem.B 
21,459(1933) 

76  Cremer,  E.,  Z.  Phys.  Chem.  B  28,  199  (1935);  Z.  Phys.  Chem.B  39,  445 
(1938) 

77  Ahlers,  G.,  J.  Chem.  Phys.  40,  3123  (1964). 

78  Nielsen,  S.E.,  Dahler,  J.S.,  J.  Chem.  Phys.  46,  732  (1967). 

79  Motizuki,  K.,  Nagamiya,  T.,  J.  Phys.  Soc.  Japan  14,  1639  (1959) 

80  Cunningham,  C.  M.,  Johnston,  H.  J.,  J.  Amer.  Chem.  Soc.  80,  2377  (1958). 

81  Motizuki,  K.,  Nagamiya,  T.,  J.  Phys.  Soc.  Japan  11,  654  (1956) 

82  Schmidt,  F.,  Phys.  Rev.  B  10,  4480  (1974). 


174 

83  Abragam,  A.,  Principles  of  Nuclear  Magnetism,  Clarendon  Press,  Oxford,  1961 

84  Esteve,  D.,  Devoret,  M.,  Sullivan,  N.S.,  J.  Phys.  C  15,  5455  (1982) 

85  Sullivan,  N.  S.,  Devoret,  M.,  Esteve,  D.,  Phys.  Rev.  B  30,  4935  (1984) 

86  Sullivan,  N.S.,  Devoret,  M.,  Cowan,  B.  P.,  Urbina,  C,  Phys.  Rev  B  17,  5016 
(1978) 

87  Zimmerman,  J.  R.,  Brittin,  W.  E.,  J.  Phys.  Chem.  61,  1328  (1957) 

88  Edzes,  H.  T.,  Samulski,  E.  T.,  J.  Magn.  Reson.  31,  207  (1978) 

89  Hahn,  E.  L,  Phys.  Rev.  80,  580  (1950) 

^Chipman,  R.  A.,  Schaum's  Outline  Series,  Transmission  Lines,  McGraw-Hill 
Book  Company,  New  York,  St.  Louis,  San  Francisco,  Toronto,  Sydney,  1968 

91  Fukushima,  E.,  Roeder,  S.  B.  W.,  Experimental  Pulse  NMR,  Addison- 
Wesley  Publishing  Company,  Reading,  Massachusetts,  1981 

92  Orr,  W.  I.,  The  Radio  Handbook,  Editors  and  Engineers,  Summerland, 
California,  14th  edition,  1956 

93  Sullivan,  N.  S.,  Pound,  R.  V.,  Cryogenics  13,  28  (1973) 

94  Andrew,  E.  R.,  Nuclear  Magnetic  Resonance,  Cambridge  University  Press, 
London,  1969 

95  Pound,  R.  V.,  Prog,  in  Nuclear  Phys.  2,  21  (1952) 

96  Frossati,  G.,  J.  de  Phys.  C6,  1578  (1978) 

97  Huiskamp,  W.  J.,  Lounasmaa,  O.,  V.,  Rep.  Prog.  Phys.  36,  423  (1973) 

98  Breck,  D.  W.,  Zeolite  Molecular  Sieves,  John  Wiley  &  Sons,  London,  Sydney, 
Toronto,  1974 

99  Smith,  J.  V.,  Mineralogical  Society  of  America,  Special  Paper  No.  1 ,  1 963 

100  Breck,  D.  W.,  J.  Chem.  Ed.  41,  678  (1964) 

101  Broussard,  L,  Shoemaker,  D.  P.,  J.  Am.  Chem.  Soc.  82,  1041  (1960) 

102  Barrer,  R.  M.,  Bultitude,  F.  W.,  Sutherland,  J.  W.,  Trans.  Farad.  Soc.  53, 1111 
(1957) 


175 

103  Greenwood,  N.  N.,  Earnshaw,  A.,  Chemistry  of  the  Elements,  Pergamon 
Press,  Oxford,  1984 

104  Breck,  D.  W.,  Smith,  J.  V.,  Scientific  American,  Jan.  1959 

105  Barrer,  R.  M.,  Zeolites  and  Clay  Minerals  as  Solvents  and  Molecular  Sieves, 
Academic  Press,  London,  New  York,  San  Francisco,  1978 

106  Schwertz,  F.  A.,  J.  Am.  Ceram.  Soc.  32,  390  (1944) 

107  Even,  U.,  Rademann,  K.,  Jortner,  J.,  Manor,  N.,  Reisfeld,  R.,  Phys.  Rev.  Lett. 
52,2164(1984) 

108  Rapp,  R.  E.,  Dillon,  L.  D.,  Godfrin,  H.,  Cryogenics  25,  152  (1985) 

109  Rail,  M.,  Evans,  M.  D.,  Sullivan,  N.  S.,  to  be  published  (1991) 

11°  Rail,  M.,  Brison,  J.  P.,  Sullivan,  N.  S.,  Phys.  Rev.  B  44,  9639  (1991) 

111  Rail,  M.,  Brison,  J.  P.,  Sullivan,  N.  S.,  Phys.  Rev.  B  44,  9932  (1991) 

112  Monod,  P.,  Cowen,  J.,  A.,  Hardy,  W.  N.,  J.  Phys.  Chem.  Solids  27,  727  (1966) 


BIOGRAPHICAL  SKETCH 


Markus  Rail  was  born  on  the  20th  of  July  1964  in  Sindelfingen,  Germany, 
as  the  son  of  Albrecht  and  Doris  Rail.  After  the  regular  13  years  of  Gymnasium 
(high  school),  he  received  the  Abitur  (Baccalaureat)  as  valedictorian  in  1983 
from  the  Albert  Einstein  Gymnasium,  Boeblingen.  Two  years  of  undergraduate 
studies  at  the  University  of  Stuttgart  were  completed  with  the  Physics  Vordiplom 
(Bachelor  of  Science)  in  October  of  1985. 

An  intermediate  15  months  of  mandatory  military  service  followed. 

He  continued  his  studies  at  the  University  of  Stuttgart  for  two  more 
semesters  before  transferring  to  the  Graduate  School  at  the  University  of  Florida 
in  the  fall  of  1987.  There  he  received  a  Master  of  Science  in  physics  in  the  fall  of 
1988  with  Prof.  N.  S.  Sullivan  and  extended  his  stay  to  pursue  the  degree  of 
Doctor  of  Philosophy. 


176 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope 
and  quality,  as  a  thesis  for  the  degree  of  Doctor  of  Philosophy. 


Neil  S.  Sullivan,  Chairman 
Professor  of  Physics 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  h  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope 
and  quality,  as  a  thesis  for  the  degree  of  Doctor  of  Philosophy. 


(I    >  - 


E.  Raymond  Andrew 

Graduate  Research  Professor  of 

Physics 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope 
and  quality,  as  a  thesis  for  the  degree  of  Doctor  of  Philosophy. 


Charles"  F.  wooper 
Professor  of  Physics 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope 
and  quality,  as  a  thesis  for  the  degree  of  Doctor  of  Philosophy. 


David  B.  Tanner 
Professor  of  Physics 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope 
and  quality,  as  a  thesis  for  the  degree  of  Doctor  of  Philosophy. 


John  Kflauder  -J 
P/rofe^sor  of  Physics  and 
Mathematics 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope 
and  quality,  as  a  thesis  for  the  degree  of  Doctor  of  Philosophy. 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope 
and  quality,  as  a  thesis  for  the  degree  of  Doctor  of  Philosophy. 


This  dissertation  was  presented  to  the  Graduate  Faculty  of  the 
Department  of  Physics  in  the  College  of  Liberal  Arts  and  Sciences  and  to  the 
Graduate  School  and  was  accepted  as  partial  fulfillment  of  the  requirements  for 
the  degree  of  Doctor  of  Philosophy. 


December  1991 


xfqsepn  H.  Simmons 
(Professor  of  Material  Science 
and  Engineering 


Dean,  Graduate  School 


